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1.0  INTRODUCTION  AND  SUMMARY 


This  report  is  concerned  with  the  narrow  V-shaped  wakes  seen 
behind  surface  ships  by  the  SEASAT  synthetic  aperture  radar.  An 
example  is  shown  in  Figure  1-1. 

The  following  characteristics  appear  to  be  typical.  The  wake 
(full)  angles  vary  greatly,  ranging  from  as  little  as  3°  to  (rarely) 
the  full  39°  of  the  Kelvin  wake  caustics.  There  seems  to  be  a 
correlation  between  wake  angle  and  ship  size  or  speed:  smaller 
(slower?)  ships  have  larger  wake  angles  than  larger  (faster?)  ships. 
The  overall  wake  lengths  can  be  large,  up  to  10  or  15  km  in  some 
cases.  There  is  apparently  no  correlation  between  the  wake  angle 
and  the  aspect  angle  of  the  ship's  track — that  is,  the  angle  between 
the  ship's  track  and  the  satellite  track.  The  complete  wake  con¬ 
sists  of  an  outer  V  which  appears  white,  reflecting  an  increased 
radar  cross  section,  together  with  an  axial  dark  strip,  reflecting  a 
decrease  in  radar  cross  section.  Often,  however,  only  part  of  this 
structure  is  visible;  one  may  see  only  one  side  of  the  white  V, 
together  with  or  without  the  dark  centerline,  or  sometimes  Just  the 
dark  centerline  itself.  No  statistics  seem  to  exist  on  how  often  no 
wake  at  all  is  visible.  No  correlation  has  been  made  between 


specific  SEASAT  images  and  specific  ships,  so  the  size  and  speed  of 


the  ships  Is  unknown — apart  from  the  qualitative  Implications  of  the 
strength  of  the  radar  return  from  the  ship  itself . 


To  understand  what  is  being  seen  by  the  SEASAT  radar,  it  is 
first  necessary  to  recall  the  characteristics  of  this  radar.  These 
are  given  in  detail  in  Section  2.0.  It  is  next  necessary  to  under- 
stand  the  surface  signatures  left  by  surface  ships;  this  occupies 
the  balance  of  the  report. 

Surface  ships  produce  three  different  wakes:  the  turbulent 
wake,  the  internal  wake,  and  the  Kelvin  wake. 

It  is  very  likely  that  the  dark  centerline  observed  by  SEASAT 
is  the  turbulent  wake.  The  turbulent  wake  is  known  to  smooth  out 
small  ambient  waves  on  the  sea  surface;  this  will  reduce  the  radar 
return  leading  to  a  dark  strip  in  the  radar  image.  The  turbulent 
wake  grows  slowly  in  size  (like  x*^  where  x  is  distance  behind  the 
ship)  and  can  persist  for  many  kilometers.  A  more  detailed  discus¬ 
sion  of  the  turbulent  wake  is  given  in  Section  5.0. 

The  internal  wave  wake  of  a  ship  should  form  a  narrow  V,  the 
outer  angle  being  determined  by  the  relative  velocity  of  the  ship 
and  the  fastest  internal  wave,  the  speed  of  which  is  determined  by 


the  local  thermocllne.  For  typical  parameters  one  would  expect  a 
full  wake  angle  of  a  few  degrees.  Qualitatively,  slower  ships 
should  have  larger  wake  angles.  In  principle,  there  should  be 
several  successive  crests.  More  details  appear  in  Section  3.1. 

Natural  Internal  waves  can  be  easily  seen  on  SEASAT  images. 
Nevertheless,  we  do  not  believe  the  V-shaped  ship  wakes  can  be 
Internal  waves,  because  ship  generated  internal  waves  are  too  weak 
by  approximately  three  or  four  orders  of  magnitude  to  explain  the 
observed  radar  returns.  The  mechanism  by  which  internal  waves 
generate  a  radar  signal  is  as  follows.  Internal  waves  produce  a 
horizontal  current  on  the  ocean  surface;  they  do  not  displace  the 
surface  vertically  (except  by  an  exceedingly  small  amount).  The 
surface  current  strains  the  ambient  surface  waves  and  changes  their 
mean  slope;  this  changes  the  radar  cross  section.  If  there  are  no 
ambient  waves,  there  is  no  radar  return  and  no  way  to  see  an 
Internal  wave.  Thus  the  change  in  radar  cross  section  produced  by 
an  internal  wave  is  proportional  to  the  ambient  cross  section 
already  present:  A o  «  eo  .  The  proportionality  constant  e  is 
(essentially)  the  surface  strain. 

The  ambient  cross  section  is  largely  Bragg  scattering  from  a 
rough  surface  at,  for  SEASAT,  a  Bragg  wavelength  of  30  cm.  Since 


the  surface  is  rough,  the  phase  change  between  scatterers  is  large. 

The  variance  in  the  cross  section  is  therefore  equal  to  the  mean: 

2  2 

<(o  -  <o»  >  *  <o>  .  Hence  to  extract  a  change  in  cross  section 

2 

of  eo  from  background  noise  requires  1/e  independent  samples, 
for  a  signal-to  noise  of  one. 


For  SEASAT,  the  number  of  Independent  samples  is  the  number  of 
pixels  in  the  wake,  not  more  than  10^.  The  values  of  e  can  be 
.1  to  .3  for  natural  internal  waves,  but  only  10-^  or  so  (depending 
on  ship  size  and  speed  and  on  the  local  thermocline)  for  a  ship. 
Hence,  the  natural  internal  waves  can  be  seen,  but  not  ship 
generated  ones,  except  in  conditions  of  exceptionally  strong  thermo- 
clines  with  correspondingly  large  values  of  £  .  Such  situations  do 
occasionally  arise.  We  have  seen  one  SEASAT  image  taken  in  the 
Strait  of  Georgia  showing  a  ship  wake  with  several  parallel  crests; 
this  is  likely  to  be  an  internal  wave  wake.  However,  the  Strait  of 
Georgia  is  well  known  as  a  region  of  huge  thermoclines;  internal 
wave  wakes  can  even  be  seen  visually  there. 


The  argument  briefly  summarized  here  is  given  in  detail  in 
Section  3.2. 
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If  the  V-shaped  wakes  are  not  produced  by  internal  waves,  then 
by  default  they  must  be  manifestations  of  the  Kelvin  wake. 

The  kinematics  of  the  Kelvin  wake  are  well  known  (and 
described  in  detail  in  Section  4.1);  the  dominant  features  are  two 
sets  of  waves,  divergent  and  transverse,  meeting  at  a  caustic  at  a 
(full)  wake  angle  of  39°  independent  of  ship  speed  or  size. 
Obviously,  however,  this  caustic  is  not  what  is  (generally)  seen  by 
SEASAT. 


The  Kelvin  wake,  of  course,  changes  the  ocean  surface  verti¬ 
cally,  so  the  radar  is  responding  to  a  quite  different  effect  than 
it  is  for  internal  waves.  When  the  slope  of  the  ocean  surface  is 
changed  by  the  Kelvin  wake,  the  Bragg  waves  that  are  responsible  for 
the  major  part  of  the  radar  cross  section  are  tilted;  this  produces 
a  shift  in  the  apparent  aspect  angle  of  the  radar  and  therefore 
changes  the  effective  Bragg  wavelength.  Because  of  the  k~^  spectrum 
of  ambient  surface  waves,  the  tilting  therefore  changes  the  radar 
cross  section.  Evidently,  then,  the  most  important  part  of  the 
Kelvin  wake  to  the  radar  is  the  region  of  maximum  slope.  We  begin 
our  Kelvin  wake  discussion  with  a  more  detailed  description  of  this 
effect  (Section  4.2). 
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The  dynamics  of  the  Kelvin  wake  are  not  well  known,  and  in 
particular,  it  is  difficult  to  calculate  quantitatively  the 
amplitudes  and  slopes  of  waves  in  the  wake  as  a  function  of  ship 
characteristics;  this  is  especially  true  in  the  near  axis  region  of 
the  wake  in  which,  it  turns  out,  is  what  we  are  most  interested. 

The  standard  treatment  is  known  as  "thin  ship  theory;”  this  is 
described  in  Section  4.1  as  a  baseline  calculation.  "Thin  ship 
theory"  is  only  an  approximation  to  the  full  linearized  surface  wave 
problem,  which  seems  to  be  known  in  the  trade  as  the  "Kelvin-Neumann 
problem."  This  is  addressed  in  Section  4.3.  Our  conclusions, 
however,  are  that  this  treatment  is  probably  also  not  adequate  to 
describe  the  near  axis  wake.  At  small  wake  angles,  it  is  quite 
likely  that  non-linear  effects  will  be  Important.  One  of  these 
which  may  be  particularly  significant  is  the  convection  by  the 
potential  flow  around  the  hull  of  short  surface  waves  produced  near 
the  bow. 

Other  sources  of  Kelvin  wakes  besides  the  hull  itself  may  also 
be  important.  In  Section  4.5  we  discuss  the  Kelvin  wake  produced  by 
the  ship's  propeller.  It  turns  out  that,  kinematically,  the 
caustics  in  propeller  wake  are  at  quite  small  angles,  a  few  degrees 
instead  of  the  39°  of  the  conventional  Kelvin  wake.  Also,  there  is 


bigger  angles.  However,  It  seems  unlikely  that  the  propeller  source 


strength  Is  adequate  to  explain  the  SEASAT  Images. 


i 

I 

In  conclusion,  we  are  left  with  a  somewhat  unsatisfying  situa¬ 
tion.  We  believe  that  narrow  wakes  cannot  be  produced  by  ship 
j  generated  internal  waves,  (except  under  extreme  conditions)  because 

these  are  too  weak  by  three  or  four  orders  of  magnitude.  Therefore, 
they  must  be  due  to  the  Kelvin  wake.  But  we  have  not  yet  been  able 
|  to  analyze  the  Kelvin  wake  sufficiently  well  theoretically  to 

identify  just  which  feature  of  the  Kelvin  wake  SEASAT  is  seeing. 

|  It  may  well  be,  then,  that  the  situation  will  only  be  cleared 

up  by  an  experiment.  Some  thoughts  as  to  what  this  should  include 
appear  in  Section  6.0. 


2.0  RADAR  ASPECTS 


2.1  Radar  Aspects  of  SEASAT-SAR  Wake  Images 


2.1.1  Introduction 


The  SEASAT  satellite  shown  in  Figure  2-1  operated  from  July  4 

to  October  10,  1978,  gathering  (along  with  other  data)  synthetic 

aperture  radar  (SAR)  images  equivalent  to  about  20%  of  the  Earth's 

8  2 

surface,  i.e.,  an  area  of  some  10  km  .  The  image  resolution 
varied  from  about  6  m  to  40  m,  depending  upon  direction  and  process¬ 
ing  as  discussed  below.  Since  such  a  radar  can  operate  night  and 
day  and  through  clouds,  it  is  clearly  a  formidable  tool  for  surveil¬ 
lance  and  remote  sensing  of  environmental  conditions  over  the  ocean. 
Table  2-1  summarizes  some  basic  characteristics  of  the  SEASAT  SAR 
system.  These  parameters  are  to  some  extent  optimized  to  observe 
ocean  surface  phenomena.  In  Figure  2-2  we  illustrate  the  SEASAT  SAR 
observational  geometry. 

As  with  many  remote  sensing  techniques,  SEASAT  SAR  images 
enable  large  scale  features  to  be  observed.  Ship  wakes  are  salient 
features  in  SEASAT  SAR  images  in  part  because  they  are  often 
observed  to  have  lengths  of  over  10  to  20  km.  Further,  the 


Scatterometer 

Antennas 


.  Link  Antenna 

Figure  2-1.  The  SEASAT  satellite  showing  the  2.2  x  10.7  m  synthetic  aperture  radar  antenna 
The  iong  axis  of  the  satellite  points  toward  the  satellite's  nadir  and  the  SAR 
antenna  observes  downward  and  to  the  right  of  the  satellite  velocity  vector. 


Figure  2-2.  Observational  geometry  for  SEASAT  SAR.  The  mean  ocean  surface  is  the  (x  \y ') 
plane.  Wave  phenomena  are  illustrated  with  the  wave  travelling  along  TT  at  an 
angle  <fr  to  the  y '  direction.  In  SAR  parlance  the  x  '-direction  is  referred  to  as  the 
range  or  cross-track  direction,  while  the  y  '-direction  is  referred  to  as  the  azimuth 
or  along-track  direction. 


observational  parameters  for  SEASAT  SAR  remain  reasonably  constant 


over  the  100  km  wide  swath,  e.g.,  the  radar  incidence  angle  9' 
varies  only  about  6°  (23+3°).  In  contrast,  aircraft  SAR  images, 
even  from  high  altitudes,  cover  much  smaller  swath  widths  and  the 
observational  parameters  vary  much  more  widely  over  the  swath.  For 
example,  an  aircraft  SAR  flying  at  25,000  ft.  (7.6  km)  altitude 
covers  a  swath  only  10  km  wide  for  incidence  angles  from  0  to  55°. 
Further,  an  aircraft  flying  at  400  knots  requires  about  8  minutes  to 
fly  100  km  while  the  SEASAT  satellite  requires  only  about  15  seconds 
to  traverse  the  same  distance.  The  point  here  is  that  a  satellite 
SAR  image  is  significantly  different  from  an  aircraft  SAR  image  and 
thus  large  scale  features  so  evident  in  satellite  SAR  images  are  not 
necessarily  so  obvious  in  aircraft  SAR  images  of  the  same  area. 

Here  in  Section  2  we  focus  on  radar  aspects  of  ship  wake 
features  in  SEASAT  SAR  images.  Review  papers  and  reports  by  Beal 
et  al.  (1981),  Fu  and  Holt  (1982)  and  Vesecky  and  Stewart  (1982) 
give  more  comprehensive  and  detailed  information  on  SEASAT  SAR 
images  of  the  ocean. 


2-5 


Formation  of  SAR  Images  From  Raw  Radar  Echo  Data  And  Display 


Of  SAR  Images 

2.2.1  Optical  and  Digital  Imaging  of  SAR  Scenes 

The  radar  echo  signal  received  by  a  synthetic  aperture  radar 
requires  extensive  manipulation  before  an  image  emerges.  We  shall 
refer  to  this  image  formation  process  as  "imaging”  and  to  manipula¬ 
tion  of  an  SAR  image  after  formation  as  image  "processing."  Both 
imaging  and  image  processing  may  be  done  optically,  (i.e.,  by  an 
optical  train,  which  is  an  analog  system)  or  digitally  (i.e.,  by 
digital  computer  manipulation  of  digital  data).  Optical  imaging  and 
processing  of  SAR  data  is  relatively  less  expensive  and  faster,  but 
also  less  accurate  and  less  flexible.  Digital  imaging  and  process¬ 
ing  of  SAR  data  is  relatively  expensive  and  slower,  but  can  produce 
the  highest  quality  Images  with  sophisticated  corrections  and  accu¬ 
rate  geographical  registration.  In  general,  for  a  40  x  40  km  SEASAT 
SAR  scene  optical  imaging  would  cost  tens  to  hundreds  of  dollars  in 
a  bulk  run  and  take  seconds  of  processor  time.  A  comparable, 
digitally  imaged  scene  would  cost  a  few  thousand  dollars  and  take 
hours  of  processor  time  on  a  medium-sized  computer,  e.g.,  a  Prime 
750  with  array  processor.  Digitally  imaged  scenes  are  generally 
required  for  detailed  analysis  of  features  such  as  ship  wakes.  Wake 


features  clearly  evident  in  digital  SAR  images  are  often  completely 
invisible  in  optical  SAR  images  of  the  same  ship  wake.  The  key 
point  here  is  that  images  of  varying  quality  result  from  the  same 
SEASAT  SAR  raw  data  depending  on  how  the  imaging  process  was  done. 

2.2.2  SEASAT  SAR  Images  -  Standard  Data  Products  Optically 
Imaged  at  Jet  Propulsion  Laboratory 

The  standard  data  product  image  for  SEASAT  SAR  was  necessarily 
optically  imaged  because  of  the  large  throughput  rate  required  to 

Q  O 

produce  images  for  all  the  SEASAT  SAR  data  collected,  some  10  km  . 
These  images  were  delivered  to  NOAA  for  public  distribution  through 
the  Environmental  Data  Information  Service.  The  resolution  in  these 
images  runs  from  25  to  40  m.  Portions  of  the  SEASAT  SAR  data  have 
been  optically  imaged  elsewhere.  For  example,  European  SEASAT  SAR 
passes  were  optically  imaged  at  the  Environmental  Research  Institute 
of  Michigan  (ERIM).  These  small  special  imaging  runs  done  at  ERIM, 
JPL,  or  elsewhere  are  usually  superior  in  quality  to  the  standard 


data  product  which  was  done  first  and  under  pressure  to  get  a  quick 


look  at  the  data. 


2.2.3  Digitally  Imaged  SEASAT  SAR  Data 


Digitally  imaged  SAR  scenes  from  SEASAT  were  first  produced  by 
Bennett  and  Cumming  (1979)  at  McDonald  Detweiler  and  Associates  in 
Vancouver,  B.C.  This  company  can  provide  high  quality,  digitally 
Imaged  SEASAT  SAR  scenes.  Their  computer  algorithm,  referred  to  as 
the  "MDA"  algorithm,  has  been  installed  on  other  computers,  e.g.  at 
DFVLR,  Oberpfaffenhofen,  W.  Germany.  Other  digital  Imaging 
algorithms  exist  or  are  under  development  at  the  Royal  Aircraft 
Establishment  (RAE),  Farnborough,  England;  The  Norwegian  Defense 
Research  Establishment,  Kjeller,  Norway;  Jet  Propulsion  Laboratory; 
Mitsubishi  Electric,  Japan;  Johns  Hopkins  Applied  Physics  Labora¬ 
tory,  and  the  Naval  Research  Laboratory.  These  various  digital 
imaging  algorithms  produce  high  quality  images,  but  a  particular 
algorithm  may  be  tuned  to  emphasize  speed  at  some  sacrifice  in 
quality.  The  RAE  algorithm  puts  emphasis  on  quality  and  many  of  the 
images  used  as  illustrations  in  this  report  were  produced  at  RAE. 

2.2.4  Properties  of  Digitally  Imaged  SAR  Scenes 


the  SEASAT  SAR  system  gathered  data  such  that  images  could  be 
produced  with  a  best  practical  resolution  of  about  25  m  along  the 
range  direction  and  6  m  along  the  azimuth  direction.  Generally  4  of 
the  6  m  resolution  cells  along  the  azimuth  direction  were  incoher¬ 
ently  averaged  to  produce  a  square  resolution  cell  of  25  x  25  m 
resolution.  Such  an  averaging  process  is  referred  to  as  '4  look' 
imaging  along  the  azimuth  direction.  Some  ship  wake  features  may  be 
entranced  by  using  the  6  m  resolution  (see  Fig.  2-4)  or  other 
imaging  options  available  with  a  given  digital  processing  algorithm. 

Many  subtleties  of  SAR  image  formation  are  discussed  in  a 
review  by  Tomiyasu  (1978).  An  example  of  these  subtleties  which  is 
Important  for  satellite  SAR  Images  of  ship  wakes  concerns  "range 
cell  migration."  This  term  refers  to  the  fact  that  for  SEASAT  SAR 
images  Earth  surface  curvature  and  other  effects  cause  a  point 
target  at  given  surface  range  (x'-  coordinate  in  Figure  2)  to 
migrate  in  echo  time  delay  over  the  SAR's  integration  period  (~  2.3s 
for  SEASAT).  In  order  to  calculate  the  doppler  spectrum  for  a  given 
range  cell,  the  echo  signal  in  different  time  delay  cells  must  be 
used.  Use  of  a  crude  single-step  algorithm  generates  'ghost' 
features  in  the  SAR  image  as  shown  in  Figure  2-3a.  These  ghosts  can 
be  removed  by  a  more  sophisticated  algorithm  to  handle  range  cell 


Note  (a)  Shows  a  4  x  4  km  arsa  near  Gotdstone,  CA,  where  an  array  of  nine  comer  reflectors  were  deployed  in  flat  dry  lake 
beds.  The  stronger  echos  from  the  two  reflectors  at  the  ends  ot  the  array  show  the  side  lobes  in  both  the  along  track 
(y)  and  cross  track  directions  For  the  digital  imaging  algorithm  used  by  RAE  to  produce  these  images  the  first  side  lobe 
is  13.3  dB  down.  However,  the  photographic  brightness  in  this  figure  is  proportional  to  echo  amplitude  and  thus  in  the 
figure  the  first  sidelobe  is  only  a  factor  4.6  less  bright  than  the  peak.  Note  that  for  the  two  end  reflectors  there  are  two 
sets  of  sidelobes  one  of  which  is  an  artifact  or  ghost.  The  wide,  nearly  vertical  stripe  near  the  right  side  ot  2-3(a)  is  the 
enormous  sidelobe  pattern  created  by  the  enormous  echo  from  NASA's  210'Goldstone  tracking  antenna,  located 


Note:  (b)  We  see  the  SAR  image  of  the  lower  left,  comer  reflector  processed  using  a  more  sophisticated  algorithm  which 
suppressed  the  ghost  sidelobe  pattern.  Use  of  the  superior  algorithm  is  said  to  "deg host"  the  image.  RAE  images  used 
in  this  report  (aside  from  2-3(a)j  are  deg  hosted  Images  kindly  provided  by  Brian  Barber  at  the  Remote  Sensing  Centre, 
Space  Department,  RAE  Famborough. 


Figure  2-3.  SEASAT  SAR  images  of  comer  reflectors  near  Goldstone.  CA.  In  both  images 
the  SAR  velocity  vector  runs  parallel  to  the  ordinate  (i.e.  y-axis  of  Figure  2-3  is 
vertical). 


$ 


migration.  Image  quality  is  Improved  as  shown  in  Figure  2-3b,  but 
at  the  expense  of  increased  computation  time.  'Deghosted'  images 
are  clearly  advantageous  in  looking  for  subtle  features  in  the  ocean 
surface  phenomena  such  as  ship  wakes .  Deghosted  images  are  used  in 
this  report. 

Dynamic  range  is  an  important  issue  in  examining  the  physics 
behind  features  in  SAR  images.  This  is  principally  because  photo¬ 
graphic  displays  of  SAR  images  can  be  misleading  in  terms  of  the 
relative  radar  backscatter  cross  section  (o°)  implied  by  dark  and 
light  portions  of  a  photograph.  For  the  SEASAT  SAR  the  dynamic 
range  of  the  receiver  front  end  amplifier  was  about  60  dB  (factor  of 
106).  A  digitized  image  using  16  bit  word  size  has  a  dynamic  range 
of  about  45  dB.  To  overcome  dynamic  range  problems  some  images  are 
printed  photographically  with  photographic  tone  proportional  to  the 
square  root  of  SAR  image  Intensity.  These  images  are  called 
amplitude  images.  Most  of  the  images  used  here  are  amplitude 
images.  Photographic  film  has  a  dynamic  range  of  at  best  about 
30  dB  and  a  photographic  print  probably  less  than  20  dB.  Hence, 
when  one  looks  at  a  photographic  print  only  a  portion  of  the  total 
dynamic  range  of  the  system  can  be  displayed.  It  is  usually  a 
mystery  just  what  the  relative  radar  echo  intensity  is  between  the 
blackest  black  and  the  whitest  white  in  a  photographic  print. 


This  brings  up  one  of  our  principal  recommendations ,  namely 
that  examination  of  SAR  data  regarding  ship  wakes  should  be  done 
using  digitally  imaged  data  in  a  digital  format.  By  viewing  such  an 
image  from  a  digital  format  via  CRT  display,  an  image  can  be  mani¬ 
pulated  in  a  known  manner  and  the  relative  intensities  of  image 
features  can  be  accurately  known.  Similarly  photographic  images 
generated  from  a  digital  format  image  using  an  appropriate  algorithm 
can  also  be  useful.  By  magnifying  small  areas  of  a  digital  format 
image,  a  photographic  (or  CRT)  image  can  be  generated  which  shows 
image  pixels  explicitly  and  their  relative  intensities  clearly.  The 
images  presented  in  this  report  are  photographs  generated  from 
digital  format  images  of  digitally  imaged  data.  However,  these 
images  were  printed  as  described  below  and  are  not  necessarily 
optimized  to  show  important  wake  features. 

2.2.5  Contrast  Stretched  Images 

In  order  to  fit  the  larger  dynamic  range  of  a  SAR  image  into 
the  smaller  dynamic  range  of  a  photograph,  a  data  compression 
process  called  linear  contrast  stretching  is  often  used.  Depending 
on  the  particular  image  the  stretching  may  be  an  expansion  or  a 
compression.  The  process  described  here  is  that  used  at  RAE  and 


corresponds  to  some  of  the  figures  used  in  this  report.  First,  a 


statistical  distribution  of  pixel  brightness  is  calculated.  Here  we 
use  brightness  B(x,y)  to  refer  to  the  brightness  of  the  photograph 
which  may  correspond  to  SAR  image  intensity  I(x,y),  i.e.  to  radar 
echo  power,  or  to  /T,  i.e.  to  radar  echo  amplitude.  Knowing  the 
histogram  or  probability  density  function  for  B(x,y)  over  a  16  bit 
(0-65,535)  count  range  consider  a  range  from  0  to  some  percentile, 
e.g.  0-98. 5%  which  contains  98.5%  of  the  image  pixels.  This  0-98.5% 
range,  which  usually  has  a  probability  density  peak  at  about  1000  to 
2000  counts,  is  divided  into  256  bins  spaced  linearly  and  these  bins 
constitute  the  24  dB  dynamic  range  of  the  photograph.  The  lowest 
photo  brightness  bin  incorporates  the  lowest  intensity  SAR  image 
pixels  as  a  black  tone.  However,  the  highest  Intensity  SAR  image 
pixels  are  excluded  from  the  photograph.  Many  of  the  Images 
displayed  in  this  report  use  the  0-98.5%  linear  contrast  stretch 
discussed  here.  Most  other  images  use  a  linear  contrast  stretch, 
but  the  parameters  and  specifications  vary. 

Although  a  specific  case  is  discussed  here,  it  is  clear  that 
the  photographic  display  of  a  SAR  image  can  compress  or  expand  the 
dynamic  range  of  the  pixels  in  the  original  image.  Thus,  until  one 
has  complete  specifications  of  just  how  the  contrast  compression  or 
expansion  was  done  for  a  particular  image,  it  is  not  possible  to 
deduce  the  relative  radar  cross  section  variations  displayed  in  a 


photographic  print  of  a  SAR  image.  Why  such  specifications  are  not 
printed  as  part  of  the  caption  on  all  SAR  image  photographs  is  hard 
to  understand.  In  any  case,  such  information  is  readily  accessible 
when  a  digitized  image  is  read  from  magnetic  tape  and  displayed  on 
an  image  analysis  system  such  as  manufactured  by  Comptal  or  Sun. 

2.2.6  Importance  of  Examining  Digital  SAR  Image  Data 

Having  discussed  the  risks  of  trying  to  draw  conclusions  from 
photographic  prints  of  SAR  images  alone,  we  make  the  point  again 
very  explicitly  that  these  pitfalls  can  be  avoided  by  using  the  full 
dynamic  range  output  of  a  digital  SAR  imaging  system.  This  full 
dynamic  range  can  be  faithfully  recorded  on  digital  media,  such  as 
magnetic  tape  or  disc,  and  subsequently  displayed  on  a  video  display 
terminal.  We  emphatically  recommend  that  such  a  system  be  used  in 
further  investigation  of  SAR  images  of  ship  wakes. 

2.3  Features  of  SEASAT  Ship  Wake  Images 

While  an  exhaustive  study  of  ship  wakes  in  SEASAT  SAR  images 


was  not  within  the  scope  of  this  study,  we  have  examined  some  ten's 
of  ship  wake  images  from  various  sources.  Below  we  describe 
features  of  these  ship  wake  images  which  we  believe  are  Important  in 


discovering  the  physical  mechanisms  which  produce  the  wake  images . 
For  the  most  part,  these  wake  features  are  illustrated  using  images 
generated  digitally  at  RAE.  These  Images  are  of  exceptionally  high 
quality. 

2.3.1  Classic  Kelvin  Wakes 

The  image  of  Figure  2-4,  illustrating  the  classic  Kelvin  wake, 
was  obtained  on  SEASAT  orbit  834  in  the  English  Channel  southwest  of 
the  Isle  of  Wight.  The  dark  linear  feature  parallel  to  the  long 
axis  of  the  ship  is  thought  to  be  associated  with  the  'turbulent 
wake'  of  the  ship.  This  makes  sense  since  a  linear  stripe  of  rela¬ 
tively  smooth  water  is  often  seen  behind  ships  and  such  relatively 
smooth  water  would  appear  dark  in  a  SAR  image  since  it  backscatters 
radar  waves  less  well  (at  oblique  Incidence)  than  the  rougher 
surrounding  water.  The  ship  (white  oblong  shape  at  right  center)  is 
displaced  from  its  expected  position  at  the  head  of  the  wake.  This 
is  a  well  known  phenomena  in  SAR  images  since  the  ship  is  moving 
relative  to  the  mean  ocean  surface.  Taking  the  dark  turbulent  wake 
feature  as  the  ship's  track,  we  can  calculate  the  ship's  velocity 
from  the  relation 

d  -  v  x'/V  (2-1) 

x 
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where  d  is  the  displacement  of  the  ship  from  its  true  position  (at 
the  head  of  the  wake),  is  the  component  of  ship  velocity  along 

the  x'-  direction  in  Figure  2-1,  x'  is  the  x'-  coordinate  of  the 
ship  and  V  is  the  satellite  velocity.  Knowing  d  *  564m  (scaling 
from  Figure  2-4),  x'  *  353,000m  &  V  •  7,000  ms  *  ,  we  have 
v »  11.2  ms  *.  Correcting  for  wake  trajectory  relative  to 
V,v-11.5ms^  (or  about  23  kts).  For  Kelvin  wakes,  the  waves 
transverse  to  the  ship  track  have  a  phase  velocity  equal  to  the  ship 
velocity.  Ue  can  calculate  the  phase  velocity  of  these  transverse 
waves  from  their  observed  wavelength  (scaled  from  Figure  2-4), 
namely  *  93m  .  From  the  deep  water  gravity  wave  dispersion 
relation  the  phase  velocity  c^  -  (gA  /  2w)^2  *  12  ms  1  .  So 
observationally  the  ship  displacement  and  transverse  wake  wavelength 
are  consistent  to  within  about  5%. 

Several  other  features  of  this  Kelvin  wake  are  of  Interest. 
First,  the  wake  is  a  weak  feature  barely  visible  above  the  speckle 
noise  of  the  image.  Second,  the  line  connecting  the  cusps  of  the 
wake  (see  Section  2.5  and  Figure  2-12  for  a  wake  description)  in  the 
upper  half  of  the  picture  is  dark  relative  to  the  background,  while 
in  the  lower  half,  this  line  along  the  wave  cusps  is  light  on  the 
outside  and  dark  on  the  inside,  i.e.,  this  signature  rises  above, 
then  falls  below  the  background.  The  angles  of  both  cusp  lines  are 


(within  the  errors  of  measurement)  at  the  classical  angle  of  near 
20°.  The  transverse  waves  of  the  wake  are  more  evident  in  the  upper 
than  the  lower  half  of  the  wake.  Finally,  (from  Figure  2-4  and 
other  Images)  the  signature  of  the  turbulent  wake  directly  behind 
the  ship  becomes  wider  and  more  pronounced  as  one  observes  further 
astern  of  the  ship,  finally  fading  away  some  7.5  km  astern  of  the 
ship.  Nothing  definite  about  the  ship  aside  from  its  observed 
speed,  heading  and  length  of  some  260  m  is  known  to  the  authors  of 
this  report. 


2.3.2  Dark  Turbulent  Wakes 


Perhaps  the  most  common  feature  in  SEASAT  SAR  images  of  ship 


wakes  is  the  dark  linear  feature  directly  astern  of  the  ship.  As 
discussed  in  2.3.1  above,  it  is  reasonable  to  associate  this  image 
feature  with  the  relatively  smooth  water  of  a  ship's  turbulent  wake. 
As  one  can  observe  from  the  images  shown  here,  it  is  evidently  the 
most  stable  feature  of  ship  wakes  in  SAR  images.  However,  the  tur¬ 
bulent  wake  signature  is  not  always  present,  e.g.,  see  Figure  2-5. 


2.3.3  Wake 8  With  Very  Narrow  Opening  Angles 

A  very  curious  feature  shown  occassionally  in  SEASAT  SAR 
images  is  illustrated  in  Figure  2-5.  Here  two  bright  lines  form  an 
opening  angle  of  only  some  6°  and  no  turbulent  wake  is  evident. 
Further  there  is  structure  in  the  upper  arm  of  the  wake  showing 
several,  roughly  parallel,  linear  features.  Perhaps  most  remarkable 
is  the  length  of  20  to  25  km.  From  the  size  of  the  ship,  it  is 
likely  to  be  a  bulk  cargo  carrier.  Based  on  a  rather  small  number 
of  observed  cases ,  it  seems  that  narrow  angle  wakes  occur  when  the 
velocity  vector  of  the  ship  is  fairly  near  the  direction  of  the  SAR 
platform  movement  (  y'-  direction  in  Figure  2-1).  To  the  author's 
knowledge  no  narrow  angle  wakes  have  been  observed  for  ships 
travelling  nearly  perpendicular  to  the  SAR  platform  movement 
(  x'-  direction  in  Figure  2-1). 

There  has  been  some  confusion  regarding  the  water  depth  below 
the  ship  in  this  figure.  Estimates  range  from  20  to  100  or  200 
fathoms.  Using  nautical  charts  and  the  full  JPL  digital  image 
(M  100  x  100  km)  we  estimate  the  depth  at  20  to  25  fathoms 
consistent  with  the  estimate  of  McCandless  (1984).  There  is  some 
uncertainty  in  the  estimate  because  only  a  small  portion  of  the 
coast  line  appears  in  the  image.  A  more  sophisticated  estimate  of 
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In  addition  to  the  very  narrow  wakes  of  Figure  2-5,  and  the 
classical  Kelvin  wake  of  Figure  2-4,  there  are  wake  images  from 
SEASAT  SAR  with  intermediate  opening  angles.  Figure  2-6  illustrates 
these  sorts  of  wakes  with  total  opening  angles  of  about  22  to  25°. 

It  is  interesting  that  this  fleet  of  small  fishing  boats  produces 
very  similar  wake  images  in  each  case.  The  dark  turbulent  wake 
features  are  absent.  In  this  image,  the  SAR  platform  moves  verti¬ 
cally  from  bottom  to  top  and  looks  from  the  left  of  the  image. 

Since  the  boats  are  moving  nearly  along  the  direction  of  the  SAR 
platform  (presumably  at  -  10  kts.  ),  we  observe  little  of  the 
displacement  of  the  ship  from  the  wake  which  is  so  pronounced  in 
Figure  2-4.  We  would  expect  maximum  displacement  of  boat  from  wake 
for  the  boat  at  lower  left  and  indeed  it  appears  that  the  strong 
boat  echo  is  displaced  vertically  upward  nearly  along  the  left  arm 
of  the  wake.  Compare  the  lower  left  boat  and  wake  with  the  boat  and 
wake  at  upper  right  when  the  strong  boat  echo  (including  sidelobe 
structure  as  in  Figure  2-2)  clearly  stands  near  the  vertex  of  the 
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Note:  This  image  was  collected  on  10  August  1978,  orbit  633  and  digitally  Imaged  at  RAE  Famborough.  The  SAR 

platform  is  moving  from  bottom  to  top,  parallel  to  the  vertical  axis  of  the  figure,  and  views  the  scene  from  left 
to  the  image.  The  approximate  processing  resolution  is  Ay- 6m  in  azimuth  and  Ax'- 18m  in  range.  The  image 
measures  3.1  km  along  the  abscissa  and  2.3  km  along  the  ordinate.  The  Image  was  kindly  provided  by  Brian 
Barber  at  the  Space  Department,  FAE  Farnborough,  U.K. 


Figure  2-6.  SEASAT  SAR  image  of  fishing  fleet  in  the  Irish  Sea  off  Balfyquintin  Pt.,  east  coast 
of  Northern  Ireland. 


wake.  The  wakes  in  the  image  are  relatively  short,  being  on  the 
order  of  a  kilometer  or  less. 

2.3.5  Wakes  With  Multiple  Components 


Two  wake  images  shown  in  Figure  2-7  show  curious  features 
going  beyond  those  described  above.  At  the  top  left,  the  wake  image 
shows  3  distinct  bright  components — two  on  the  left  forming  a  V  with 
a  -  5°  opening  angle  and  one  brighter  line  running  nearly  verti¬ 
cally  in  the  figure.  This  bright  line  forms  an  opening  angle  of 
~  20°  with  the  nearer  arm  of  the  narrow  V.  The  viewing  geometry  is 
virtually  the  same  as  with  Figure  2-6,  so  the  ship  making  this  wake 
is  apparently  travelling  nearly  parallel  to  the  SAR  velocity  vector 
V  ,  and  we  expect  little  displacement  of  ship  from  wake.  The  ship 
image  itself  appears  aligned  with  the  Isolated  and  very  bright  arm 
of  the  wake.  No  dark  turbulent  wake  is  visible.  Note  that  the 
background  against  which  this  wake  is  viewed  is  very  dark  indicating 
a  very  calm  sea  surface,  i.e.,  low  wind  speed. 

At  bottom  left  of  Figure  2-7,  we  find  a  ship  moving  almost 
directly  away  from  the  radar.  As  in  Section  2.3.1  above,  we 
calculate  the  ship's  speed  from  the  displacement  between  the  ship 
and  the  turbulent  wake.  The  displacement  d  "  250m;  hence  V  ■  6m/s 
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Note:  This  image  was  collected  on  orbit  782,  19  August  1978  and  digitally  Imaged  at  RAE  Famborough.  The  SAR 

platform  is  moving  from  bottom  to  top,  parallel  to  the  vertical  axis  of  the  figure,  and  views  the  scene  from  left 
of  the  image.  The  approximate  processing  resolution  Is  ay- 25m  in  the  azimuth  direction  and  Ax -23m  in  the 
range  direction.  The  image  measures  31.25  Km  along  both  the  abscissa  and  the  ordinate  The  image  was  kindly 
provided  by  Brian  Barber  at  the  Space  Department,  RAE  Famborough,  U.K. 


Figure  2-7.  SEASAT  SAR  image  showing  wakes  with  multiple  components  in  the 
Mediterranean  Sea  south  of  Monaco. 


or  12  kts.  The  wake  has  three  different  components:  an  upper  arm 
making  an  angle  of  ~  20°  with  the  central  turbulent  wake  component, 
and  finally  a  brighter  patchy  line  below  making  an  angle  of  about 
17°  with  the  central  turbulent  wake.  This  lower  arm  does  not  appear 
to  be  displaced  from  the  ship — this  arm  and  the  ship  both  appear  to 
be  displaced  relative  to  turbulent  wake  and  upper  arm. 

Finally,  a  long  turbulent  wake  is  seen  near  the  image  center. 
It  stretches  some  20  to  25  km  from  center  toward  upper  left  even¬ 
tually  dispersing  in  the  smoother  water  at  upper  left. 

2.3.6  Linear  Features  in  Wakes 

An  interesting  feature  of  the  bright  lines  in  wake  images  is 
the  very  fine  nature  of  the  lines.  In  most  cases,  the  lines  forming 
the  wake  are  some  several  resolution  cells  or  less  in  width.  This 
suggests  that  the  phenomena  causing  these  very  fine  lines  is  very 
narrowly  confined  spatially  along  the  water  surface.  This  suggests 
a  sharply  "resonant"  phenomena.  Hence  one's  attention  is  directed 
toward  items  in  wake  hydrodynamics,  radar  wave  scattering  from  the 
sea  surface  and  SAR  system  response  which  show  a  strong  functional 
dependence  on  some  relevant  parameter.  One  example  is  the  depen¬ 
dence  of  radar  backscatter  from  a  rough  surface  on  angle  of 


Radar  Backscatter  From  The  Ocean 


Scattering  Mechanisms 


Scattering  of  electromagnetic  waves  from  a  rough  surface  is  a 
difficult  problem  and  is  amenable  to  analytic  treatment  only  over  a 
limited  range  of  conditions.  Thus,  even  if  one  could  specify  the 
ocean  surface  exactly,  calculating  the  radar  echo  from  this  surface 
analytically  could  still  be  difficult.  We  point  this  out  because 
some  'simulated  SAR  Images'  have  simply  been  plots  of  some  ocean 
surface  property  and  did  not  include  the  radar  wave  scattering  part 
of  the  problem.  These  comments  notwithstanding,  analytical  models 
are  clearly  very  helpful  and  are  reasonably  accurate  over  a  wide 
range  of  conditions.  Below  we  briefly  discuss  some  analytical 
models  for  electromagnetic  wave  scattering  from  the  ocean  surface 
and  comment  on  their  relevance  to  the  imaging  of  ship  wakes.  An 
excellent  review  of  these  mechanisms  is  given  by  Hasselmann  et  al. 
(1984). 


2.4. 1.1  Bragg  Resonance  Scattering.  In  this  mechanism,  the 
incident  radar  wave  resonates  with  a  'lattice'  of  ocean  surface 
waves  which  satisfy  a  Bragg  resonance  condition,  namely 
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where  is  the  Bragg  resonant  wavenumber  and  k  the  incidence  radar 
wavenumber.  This  is  the  first  order  Bragg  resonance  which  generally 
dominates  this  scattering  process.  However,  in  some  circumstances 
higher  order  resonances  may  be  important  especially  if  the  directional 
distribution  is  broad  and/or  anisotropic.  In  this  Bragg  resonance 
process  the  radar  wave  acts  as  a  filter  on  the  surface  roughness 
selecting  out  waves  which  have  a  wavelength  -  1/(2  sin  6) 

(where  1  is  the  radar  wavelength)  and  are  traveling  either  directly 
toward  or  away  from  the  radar.  The  radar  backscatter  cross  section 
resulting  from  this  process  can  be  written  (Wright,  1968)  as 

o°  (0',  k)±J  -  16  ir2k4  cos40f  jg±j  ( ©  * ,  er)|2  ♦  (Kx-»  Ky*)  (2-3) 

where  +  is  the  directional  waveheight  spectrum  evaluated  at 

K  ,  «  2  k  sin  0,  K  „  -  0,  g  is  the  radar  reflectivity  of  a 
x  y  l  j 

planar  air-sea  water  interface,  and  ij  denotes  the  incident  and 
backscattered  polarizations.  For  the  SGASAT  SAR  ij  ■  HH,  i.e., 
horizontal  polarization.  The  most  important  feature  of  equation 
(2-3)  is  that  o°  is  proportional  to  the  waveheight  variance  spectrum 
evaluated  at  a  specific  K  and  is  thus  proportional  to  the  energy  in 


the  Bragg  resonant  waves 


This  scattering  process  is  thought  to  dominate  microwave 
backscatter  from  the  ocean  at  0'  ■  20°  to  70°  as  shown  in  Figure  2-8 
where  the  Bragg  resonance  model  is  labeled  as  the  'slightly  rough' 
curve.  Strictly  speaking  this  model  is  valid  only  for  a  slightly 
rough  surface  in  the  sense  that  |k?  cos  9  * |  «  1  where  C  is  the  rms 
height  deviation.  The  Bragg-resonance  model  is  supported  not  only 
by  the  argeement  with  observations  shown  in  Figure  2-8,  but  also  by 
measurements  of  the  radar  echo  frequency  spectrum  which  show  strong 
peaks  at  frequencies  corresponding  to  the  expected  Doppler  shift  of 
the  first  order  Bragg  waves  (Graf,  et  al.  1977).  Some  evidence  for 
higher  order  phenomena  were  presented  by  Gotwols  et  al.  (1984). 

2.4. 1.2  Quasl-Specular  Scattering.  This  scattering  mechanism 
applies  to  situations  where  C  nay  be  large,  but  the  surface  is 
gently  undulating  (radius  of  curvature  »  radar  wavelength).  The 
term  quasi-specular  refers  to  scattering  from  a  statistical  ensemble 
of  facets  (with  randomly  distributed  normals)  some  of  which  are 
oriented  so  as  to  backscatter  via  specular  reflection  (for  reference 
see  Bass  &  Fuks,  1979,  chapter  7).  With  reference  to  the  ocean 
surface  this  mechanism  is  thought  to  dominate  microwave  backscatter 
from  the  sea  at  6'  <  20°,  i.e.,  at  near  normal  incidence.  The 

radar  cross  section  for  the  quasi-specular  mechanism  is  given  in 
concise  form  by  Kodis  (1966)  as 


(2-4) 


0°  -  n  R1R2  <  N  > 

where  Rj  and  R2  are  the  principal  radii  of  curvature  of  the  surface 
and  <  N  >  is  the  number  of  facets  per  unit  area  properly  tilted  to 
yield  specular  reflection.  For  an  isotropic  surface  with  a  Gaussian 
distribution  of  height  fluctuations  Barrlck  (1968)  find  the  quasi- 
specular  a 0  to  be 

2 
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where  R(0)  is  the  Fresnel  reflection  coefficient  for  normal  inci- 
dence  and  s  is  the  total  slope  variance  for  ocean  waves  such 
that  A  >  A  .  To  obtain  the  'specular-point'  curve  in  Figure  2-8, 
R(0)  in  eq.  (2-5)  must  be  replaced  by  an  'effective'  value  which  is 
about  one-half  the  expected  value  of  R(0)  -  1.  This  "effective" 
value  of  |r(0)|  Is  probably  due  to  small  scale  roughness  on  the 
specular  point  facet  surface.  Such  roughness  would  reduce  | R( 0) | . 

O 

Further,  the  value  of  s  must  be  that  of  a  dirty  ocean  surface  where 
short  waves  are  more  suppressed  by  viscous  damping. 

2.4. 1.3  Wedge  Scattering  by  Cusped  Waves.  When  waves  are 
near  breaking  they  can  be  modeled  in  an  approximate  sense  as  a  wedge 
having  a  120°  included  angle  corresponding  to  the  limiting  slope  for 
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surface  gravity  waves.  Scattering  from  the  edges  of  such  wedges 
(local  radius  of  curvature  «  radar  wavelength)  can  be  treated 
analytically.  If  one  distributes  wedges  randomly  over  a  model  ocean 
surface  with  an  appropriate  density,  backscattering  from  cusped 
waves  may  be  estimated.  Lyzenga  et  al.  (1983)  have  performed  such  a 
calculation.  Their  results  show  that  backscattering  from  cusped 
waves  is  important  at  large  angles  of  incidence.  Comparison  with 
ocean  backscatter  measurements  shows  that  for  horizontal  polari¬ 
zation  at  X-band,  wedge  scattering  from  cusped  waves  dominates 
o°  at  0'  >  60°  and  agrees  much  better  with  measured  o°  for  high 
winds  (U  >  10  m/s).  At  L-band  for  horizontal  polarization  the 
cusped  wave  model  again  dominates  backscattering  for  8'  >  60°,  but 
yields  a  o°  somewhat  above  measured  values.  Wedge  backscattering 
from  cusped  waves  is  not  important  for  vertical  polarization  except 
for  8 '  very  near  grazing  incidence  where  shadowing  may  very  well 
reduce  its  importance. 


2. 4. 1.4  Two-scale  Scattering  Models.  Two-scale  scattering 
models  are  almost  always  used  in  calculating  the  radar  backscatter 
characteristics  of  the  ocean  surface.  In  connection  with  SAR  images 
of  ocean  surface  phenomena  two  distinct  two-scale  models  are 
relevant:  the  electromagnetic-hydrodynamic  (emh)  model  and  the  SAR 
model.  In  both  these  models  radar  backscatter  is  fundamentally  due 


to  radar  wave  interaction  with  small  scale  roughness  (<  10  X). 
However,  the  large  scale  behavior  of  the  ocean  is  allowed  to 
influence  the  properties  of  this  small  scale  roughness  in  particular 
ways.  The  dichotomy  between  large  and  small  scale  ocean  surface 
behavior  is  illustrated  in  Figure  2-9. 

In  the  emh  two-scale  model  Bragg  resonance  scattering  arises 
from  the  small  scale  roughness  (K^  »  Krmh  where  Kemh  is  the 
separation  wavenumber  for  this  model).  The  behavior  of  these  small 
scale  ripples  is  governed  in  terms  of  a  spectral  transport  equation 
where  the  large  scale  waves  advect  the  small  scale  waves  by  their 
orbital  velocities  and  refract  them  due  to  horizontal  shear  in  the 
orbital  velocities  of  the  long  waves.  A  spectral  transport  equation 
could  incorporate  source  terms  corresponding  to  generation  by  wind, 
energy  loss  due  to  damping  and  energy  transfer  via  wave-wave  (non¬ 
linear)  interaction. 

The  SAR  two-scale  model  divides  the  wave  spectrum  at  the  SAR 
resolution  scale  size  (Koa_).  Large  waves  (K  <  Koa_)  are  treated 
deterministically  and  small  waves  (K  >  K_a_)  are  treated  statlsti- 
cally.  As  it  turns  out,  the  two  models  may  be  combined  by  statisti¬ 
cally  averaging  the  results  of  the  emh  model  over  the  intermediate 
scale  waves  shown  in  Figure  2-9. 
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Figure  2-9.  Dichotomy  of  ocean  surface  behavior  into  large  and  small  scales  divided  by  a 
separation  wavenumber  (Kemh  for  the  emh  model,  for  the  SAR  model).  The 
first-order  Bragg  resonant  wavenumber  Kb  =  2  k  sin  6.  Generally  K,ar  <  Ymh 
<  Kb  and  for  SEASAT  (2  ir/25)m'1,  Kemh®  (2ir/3)m''  and  Kb®  (2ir/0.3)m1. 
Figure  after  Hasseimann  et  al.  (1984). 


Gravity  Waves 


Ocean  waves  of  length  10' s  to  100's  of  meters  are  often 
visible  in  SAR  ocean  images.  Three  major  mechanisms  are  thought  to 
contribute  to  SAR  imaging  of  long  (K  <  Kgar)  surface  gravity  waves, 
namely  the  tilt ,  hydrodynamic  and  velocity  bunching  mechanisms. 
Considerable  research  has  been  done  to  construct  mathematical  models 
corresponding  to  these  mechanisms  (see  reviews  by  Beal,  1982  and 
Vesecky  &  Stewart,  1982).  The  tilt  mechanism  refers  to  the  two- 
scale  models  above  in  which  small-scale  Bragg  (or  quasi-specular) 
scattering  regions  some  few  meters  in  size  are  tilted  with  respect 
to  the  radar  ray  path  by  large-scale  ocean  waves.  Since  those  Bragg 
(or  quasi-specular)  scattering  patches  on  the  ocean  more  nearly 
perpendicular  to  the  radar  ray  path  backscatter  more  strongly  (see 
Figure  2-8),  they  are  brighter  in  the  image.  Thus,  for  ocean  waves 
approaching  the  radar  the  portion  of  the  wave  somewhat  in  advance  of 
the  crest  would  appear  brightest  and  the  portion  somewhat  following 
the  crest  would  appear  darkest.  The  hydrodynamic  mechanism  refers 
to  the  straining  of  the  water  surface  by  a  passing  large-scale 
wave.  The  region  somewhat  in  advance  of  the  wave  crest  experiences 
a  velocity  convergence  which  compresses  ambient  small-scale  waves, 
thus  raising  their  amplitude  and  hence  their  radar  cross  section 


(see  eq.  2-3).  This  portion  of  the  wave  should  then  appear  brighter 
in  the  SAR  image.  The  velocity  bunching  mechanism  is  a  phenomenon 
unique  to  SAR  since  it  relies  on  Doppler  shift.  The  Doppler  shift 
introduced  by  the  relative  motion  between  the  SAR  platform  and  a  SAR 
pixel's  surface  location  allows  the  SAR  to  correctly  associate  back- 
scattered  energy  with  a  given  pixel  location.  The  orbital  velocity 
associated  with  a  large-scale  wave  introduces  a  periodic  perturba¬ 
tion  of  the  basic  Doppler  shift,  and  thus  periodically  shifts,  i.e., 
misplaces,  SAR  image  brightness  along  the  azimuth  direction 
(y'~  direction  in  Figure  2-2).  This  periodic  misplacement  could 
then  produce  a  wave  image  in  the  SAR  output.  There  are  other 
mechanisms  which  could  conceivably  affect  SAR  imaging  of  ocean 
waves,  e.g.,  cusped  waves,  breaking  waves,  foam,  etc.  However,  no 
quantitative  investigations  of  these  mechanisms  have  yet  been  done. 

The  modulation  transfer  function  R  relates  the  change  in  radar 
cross  section  (Ao°)  associated  with  a  large-scale  ocean  wave  to 
the  wave '8  physical  characteristics  by 

Aa° 
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where  y  -  arctan  [Im(R)/Re(R) ]  and  A,  K  and  fl  are  the  long  wave's 
amplitude,  vector  wavenumber  and  frequency;  r  is  distance  in  the 


(x',y')  plane.  It  is  interesting  to  note  that  when  R  is  calcu¬ 
lated  as  a  function  of  0*  for  the  tilt  mechanism  alone,  there  are 
two  maxima  and  one  maximum  is  near  8'  *  20®,  i.e.,  near  the  angle 

of  incidence  at  which  SEASAT  SAR  operated  (Bahar,  1983c  and  Vesecky 
et  al.,  1983).  This  implies  that  SEASAT  SAR  was  particularly  sensi¬ 
tive  to  small  tilts  of  the  ocean  surface.  Comparing  the  general 
visibility  of  ocean  wave  fields  between  the  first  shuttle  Imaging 
radar  (SIR-A)  which  operated  at  8’  ~  45°  and  SEASAT  we  find  waves 
much  more  prevalent  for  SEASAT  in  agreement  with  the  above  implica¬ 
tions.  Further,  the  aforementioned  calculations  of  R  indicated  that 
quasi-specular  as  well  as  Bragg-resonance  scattering  were  important 
for  SEASAT  SAR  observations  of  ocean  waves. 

2.4.3  Spatial  Modulation  of  Ocean  Roughness  by  Internal  Waves 

One  of  the  surprises  of  the  SEASAT  mission  was  the  almost 
ubiquitous  nature  of  Internal  wave  surface  effects  in  SAR  ocean 
Images.  Particularly  along  coasts  (out  to  the  continental  shelf) 
and  over  sea  mounts  strong  SAR  images  of  internal  wave  packets  were 
very  frequently  observed.  An  example  is  shown  in  Figure  2-10.  In 
many  cases,  such  as  Figure  2-10,  internal  waves  are  presumably 
associated  with  strong  tidal  currents  as  they  encounter  shallow 


water  near  coasts 


Note:  This  image  was  collected  17  September  1978.  at  about  1720  UT  near  79“  45'N  latitude.  114“  40W  longitude.  The  image 
is  approximately  100  km  top  to  bottom.  The  wavelengths  in  the  internal  wave  packets  range  from  about  0.5  to  5  km. 
The  white  specks  are  defects  in  the  archive  negative.  Spacecraft  ground  track  was  parallel  to  the  bottom  of  the  image 
going  from  right  to  left  along  340".  This  digital  image  was  kindly  provided  by  Jet  Propulsion  Laboratory. 


Figure  2-10.  SEASAT  SAR  image  (orbit  1183)  of  internal  waves  off  the  island  of  Angel  de 
la  Guarda  in  the  Gulf  of  California. 


Modulation  of  surface  waves  by  internal  waves  and  surface 

currents  has  been  examined,  Graf  et  al.  (1979),  in  laboratory 

experiments  and  theoretically  by  Zachariasen  (1973),  Gargett  & 

Hughes  (1981)  and  Phillips  (1973,  1981).  Phillips  (1981)  notes  that 

small-scale  Bragg  resonant  surface  waves  can  interact  strongly  with 

Internal  wave  surface  currents  provided  the  internal  and  surface 

waves  are  traveling  in  the  same  direction  and  the  group  velocity 

(c„)  of  the  surface  waves  is  near  resonance,  i.e.,  equal  to,  the 
8s 

phase  velocity  (Cp^)  of  the  internal  wave.  For  c^s  *  c  ^ 

internal-wave  Induced  surface  currents  can  modulate  the  SAR  resonant 

surface  waves,  i.e.,  ripples  with  K  -  K^.  In  shallow  water,  such  as 

along  the  continental  shelf,  near  islands  and  over  sea  mounts, 

internal  wave  amplitudes  can  be  large  and  induced  surface  current 

velocities  (u)  also  quite  large,  u  *  (c  c  ./h)  where  C  is 

max  pi  max 

the  maximum  thermocllne  displacement  and  the  thermocline  depth  is  h. 
Large  u  means  significant  modulation  of  small-scale  ripples  can 
occur  over  a  wider  range  of  (cpi  "  cg8)  ®nd  SAR  detection  of  the 
resulting  o°  modulation  is  more  likely.  Nevertheless,  Phillips 
(1981)  finds  it  remarkable  that  internal  wave  Images  are  so 
frequently  observed  by  SEASAT  SAR  in  continental  shelf  waters. 


One  clear  implication  of  the  mechanism  discussed  above  is  that 
internal  waves  traveling  approximately  in  the  range  direction  (y“- 


direction  in  Figure  2-2)  are  more  likely  to  be  imaged  than  those 
traveling  along  the  azimuth  (or  x'-  direction  in  Figure  2-2). 

Indeed  this  phenomena  occurs  in  Figure  2-10  where  we  observe  very 
bright  wave  crests  for  internal  wave  packets  traveling  approximately 
in  the  range  direction  at  center  right  and  upper  right.  Why  should 
this  be  so?  The  group  velocity  of  Bragg  resonant  waves  for  SEASAT 
SAR  is  -  35  cm/s.  Typical  internal  wave  phase  velocities  are  also 
~  10's  of  cm/s.  So  if  the  internal  wave  is  traveling  in  nearly 

the  same  direction  as  the  Bragg  resonant  ripples ,  then  relatively 
strong  Interaction  between  the  two  is  likely  and  relatively  strong 
modulation  of  surface  roughness  at  A  -  A^  is  likely  to  occur.  In 
Figure  2-10  we  expect  this  to  occur  for  internal  waves  traveling 
approximately  along  the  ordinate  of  the  figure.  We  associate  the 
bright  wave  packet  features  mentioned  above  with  this  mechanism  for 
modulation  of  small-scale  surface  roughness  particular  because  the 
bright  portions  of  the  wave  packet  rise  significantly  above  the 
background  <  o°  >.  Internal  waves  packets  traveling  in  other 
directions  in  Figure  2- 10  tend  to  be  more  a  darkening  relative  to 
the  background  <  o°  >  . 

The  mechanism  described  above  would  tend  to  enhance  SAR  image 
brightness  over  regions  of  internal-wave  induced  downwelllng 
(converging  surface  currents)  and  visa-versa  for  upwelling  regions. 
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Similar  enhancements  would  be  produced  if  some  type  of  radar  baclc- 
scatterlng  objects;  such  as  flotsam,  foam,  debris,  etc;  were 
concentrated  by  the  downwelling. 

An  alternative  mechanism  would  have  surface  films  being 
concentrated  in  current-convergence  regions ,  damping  the  local  SAR 
resonant  ripples  and  hence  reducing  SAR  image  brightness.  In 
current  divergent  regions  the  dispersal  and  breaking  of  surface 
films  would  allow  higher  energy  density  levels  of  SAR  resonant  waves 
and  hence  a  SAR  image  brightness  at  the  background  level  or  above. 

We  suggest  that  this  process  might  dominate  for  the  internal  waves 
traveling  more  nearly  along  the  abscissa  in  Figure  2-10  since  the 
SAR  resonant  ripples  would  be  traveling  nearly  perpendicular  to  the 
Internal  wave  surface  currents  and  thus  would  not  interact  strongly 
with  them. 

At  present  the  question  of  a  physical  mechanism  for  the 
visibility  of  internal  waves  in  SAR  images  must  remain  open  for  two 
reasons.  First,  to  the  author's  knowledge,  no  theoretical  attack 
has  yet  been  mounted  Incorporating  both  the  necessary  hydrodynamics 
and  radar  scattering  theory.  Second,  no  experimental  data  set 
containing  both  SAR  and  appropriate  in-situ  hydrodynamic  measure¬ 
ments  has  yet  been  analyzed.  Data  from  Boundary  Passage  collected 


during  Che  Georgia  Strait  experiment  provide  the  opportunity  to 
attack  the  latter  point  of  the  two  mentioned  above. 

Some  Interesting  empirical  data  from  a  comprehensive  survey  of 
SEASAT  SAR  Images  of  natural  internal  waves  in  the  northeast 
Atlantic  was  presented  by  Bagg  (1983).  After  examining  some 
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5  x  10  km  of  SAR  Images  in  this  region  Bagg  concluded  that: 

a.  Internal  wave  surface  features  can  be  imaged  by  SAR  in 
water  depths  ranging  from  very  shallow  to  very  deep. 

b.  Internal  wave  surface  features  can  be  Imaged  with 
winds  ranging  from  light  airs  to  gale  force. 

c.  The  likelihood  of  imaging  internal  waves  with  SAR 
decreases  with  increasing  water  depth  and/or  wind 
speed. 

d.  The  effect  of  increasing  wind  speed  causes  an 
increasing  limitation  on  internal  wave  imaging  as  the 
water  depth  Increases. 

Further  Information  on  Internal  wave  images  in  this  region  is 
supplied  by  Kaslschke  et  al.  (1983). 


2.4.4  Other  Surface  Roughness  Modulation  Observed  by  SAR 

2.4.4. 1  Winds .  Since  surface  wind  stress  is  the  primary 
source  of  radar  resonant  ripples,  radar  o®  is  clearly  a  function  of 
local  wind  speed.  Since  radar  measurement  of  ocean  winds  is 
normally  done  at  band  (X  ~  2  cm),  little  work  with  L-band 
(X  -  20  cm)  SAR  has  been  done.  Vesecky  and  Stewart  (1982)  review 
research  in  this  area.  Two  principal  research  results  standout: 

a.  SAR  a0  at  L-band  is  only  weakly  dependent  on  wind 
speed,  o°  being  approximately  proportional  to 
where  U  is  the  neutral  stability  wind  at  height  b 
above  the  surface  and 

Y  *  0.2  for  b  -  10  m  and  y  "  0.4  for  b  -  19  m. 

b.  SAR  o°  at  L-band  is  insensitive  to  wind  direction. 
(Although  roll  instabilities  may  cause  large  scale, 

-  few  km,  structures  aligned  with  the  wind 
direction.) 

However,  this  relatively  weak  sensitivity  of  a 11  to  wind  speed 


does  not  prevent  regions  of  differing  wind  speed  from  being  very 


evident  in  SAR  images.  Examples  in  Beal  et  al.  (1981)  and  Lyden  et 
al.  (1983)  aptly  demonstrate  this  point. 

2.4. 4. 2  Surfactants.  It  is  well  known  that  surfactant  films 
strongly  damp  surface  waves.  For  example,  Phillips'  (1977)  calcula¬ 
tions  indicate  that  Bragg  resonant  waves  for  SEASAT  SAR 
(A  ~  30  cm)  have  an  energy  density  decay  (e-folding)  time  of 

D 

t  "*  450  s  for  a  perfectly  clean  ocean  surface,  but  t  -  16  s  for 
a  surface  covered  by  an  inextensible  surface  film.  Features  in  SAR 
images  which  are  apparently  related  to  slicks  are  commonly  seen.  A 
particularly  striking  example  is  shown  in  Figure  2-11  where  a  ship 
Is  apparently  discharging  oil.  Little  experimental  research 
utilizing  SAR  images  and  surfactants  with  measured  characteristics 
has  been  done.  Huhnerfuss  et  al.  (1983)  report  experiments  with 
mono-molecular  films  observed  by  L  and  X-band  aircraft  SAR. 
Surfactant  films  reduced  X-band  radar  cross  section  by  a  factor  of 
-  15  to  200Z  ;  while  L-band  cross  section  was  reduced  by  about 
20Z.  Differing  chemical  structure  in  the  surfactant  does  alter  the 
effect  on  radar  cross  section.  Further,  Huhnerfuss  et  al.  (1981) 
have  investigated  the  fundamental  issue  of  surface  wave  attenuation 
by  artificial  surface  films  of  different  chemical  structure. 
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Figure  2-11.  Ship  in  the  Atlantic  ocean  apparently  discharging  oil  or  other  surfactant  which 
damps  the  surface  ripples  to  which  SAR  is  sensitive.  This  image  was  collected 
19  Sept.  1978  at  1442  GMT  about  500  km  out  from  Delaware  Bay  by  SEASAT 
SAR  on  orbit  1210.  After  Fu  &  Holt  (1982). 


2. 4. 4. 3  Spray,  Foam  and  White  Caps.  Spray,  foam  and  white 
caps  can  clearly  contribute  to  microwave  scattering  from  the  ocean 
surface.  However,  not  much  is  known  quantitatively  about  these 
processes.  SAR  images  of  breaking  waves  indicate  that  the  bright 
wave  crests  are  observed  by  the  SAR  at  a  Doppler  frequency  corre¬ 
sponding  approximately  to  the  wave  phase  velocity  (Lyzenga,  et  al. 


SAR  System  and  Wave  Scattering  Effects  Upon  Ship  Wake  Images 


2.5.1  Effects  Upon  Wakes  in  General 

There  are  several  SAR  system  and  wave  scattering  effects  which 
are  relevant  to  SAR  ship  wake  images  irrespective  of  the  physical 
mechanisms  which  cause  the  wakes  to  be  imaged.  These  are  discussed 
first  in  this  section  while  other  items  associated  more  specifically 
with  surface  or  internal  waves  are  covered  later  in  Sections  2.5.2 
and  2.5.3,  respectively. 

2 . 5 . 1 . 1  Contrast  Stretching  and  Fine  Wake  Features. 


Curiously  the  linear  features  which  form  the  bright  lines  in  SAR 
ship  wake  images  are  very  fine,  i.e.,  very  narrow.  Estimates  of  the 
width  of  these  lines  run  to  several  SAR  resolution  cells  or  less. 

If  indeed  these  lines  are  very  fine  features  contrasting  strongly 
with  the  background,  then  this  is  an  important  feature  of  SAR  ship 
Images  that  must  be  explained  by  any  candidate  model  for  the  wake 
Imaging  mechanism.  However,  many  of  the  estimates  of  wake  feature 
characteristics  are  based  on  photographic  prints  which  have  under¬ 
gone  a  "contrast  stretching"  process  described  in  Section  2.0  above. 
This  process  could  give  the  impression  that  wake  lines  are  very 
fine,  high  contrast  features  when  in  fact  the  backscattered  energy 


returned  by  the  features  could  be  only  tens  of  percent  above  the 
background.  Therefore,  It  is  necessary  to  know  what  contrast 
stretching  has  been  applied  to  an  image  before  making  any  conclu¬ 
sions  regarding  the  fineness  of  wake  features  in  the  image.  To  make 
quantitative  measurements  of  wakes  images  digital  format  data  and 
image  processing  apparatus  are  required  as  discussed  in  Section  2.0 
above . 


2 . 5 . 1 . 2  Spatial  and  Temporal  Limits  to  Coherent  Scatter. 
Fresnel  zone  size  (appropriately  calculated)  determines  approxi¬ 
mately  the  maximum  size  scale  over  which  an  individual  scattering 
element  (facet)  on  the  ocean  surface  can  backscatter  coherently.  In 
the  temporal  domain  an  element's  backscattering  properties  must  not 
fluctuate  significantly  over  the  integration  time  (T)  of  the  SAR 
system  in  order  to  contribute  coherently  to  the  received  SAR  signal. 
While  integration  time  for  a  SAR  system  can  be  varied  (up  to  some 
maximum  limit  related  to  the  system  bandwidth)  after  signal  collec¬ 
tion,  the  Fresnel  Zone  size  is  fixed  by  the  SAR  wavelength  and  the 
observational  geometry.  Thus  Fresnel  zone  size  could  be  an  impor¬ 
tant  parameter  in  the  scattering  of  SAR  signals  from  the  ocean. 
Fresnel  zone  size  for  a  plane  surface  perpendicular  to  k  ,  the 
radar  wave  vector,  is  -  /  AR/4. 


SAR  System 
SEASAT 


R 


X 


Typical 

Resolution 


850  km  23  cm  (L-band)  221  m  25  m 

High  Altitude 

Aircraft  SAR  10  km  23  cm  (L-band)  24  m  10  m 

High  Altitude 

Aircraft  SAR  10  km  10  cm  (X-band)  16  m  3m 

Low  Altitude 

Aircraft  SAR  3  km  23  cm  (L-band)  13  m  5m 

Low  Altitude 

Aircraft  SAR  3  km  10  cm  (X-band)  9m  2m 


We  see  from  the  table  that,  in  fact,  SAR  resolution  cell  size  is 
usually  much  smaller  than  the  Fresnel  zone  size.  So,  in  practice, 
it  is  the  SAR  resolution  cell  size  (discussed  below)  which  limits 
the  spatial  scale  of  coherent  backscatter. 

2.5. 1.3  Effects  of  Resolution  Cell  Size.  Referring  to 
Figure  2-9  we  note  that  in  the  SAR  two-scale  model  the  wavenumber 
KSar  (*  tt/SAR  resolution)  separates  the  small  scale  waves  K  > 


which  must  be  treated  statistically  and  the  large  scale  waves 


K  <  Kgar  which  are  resolved  by  the  SAR  and  can  be  treated  determi¬ 
nistically.  For  the  Bragg  resonance  scattering  mechanism  the 
individual  scattering  facets  (consisting  of  sets  of  ~  10  resonant 
waves  A  »  A^)  will  be  tilted  by  intermediate  scale  waves 

K  <  K  <  K  .  .  In  some  situations,  e.g.,  horizontal  polarization 
sar  emh 

and  angles  of  incidence  6'  <  30°  it  is  important  to  average  the 
contributions  of  the  Bragg  scattering  facets  over  the  slope  proba¬ 
bility  distribution  of  these  intermediate  scale  waves.  Clearly  the 
value  of  K  is  a  determining  factor  in  this  averaging  process. 

This  process  is  discussed  further  below. 

2. 5. 1.4  Wave-Wave  and  Wave-Current  Interaction  Hydrodynamics . 
The  action  spectral  density  N(&,  x,  t)  (-  E(K)/£1',  where  E  is 
wave  energy  spectral  density  and  (2*  is  the  intrinsic  wave  frequency) 
is  commonly  used  to  describe  the  state  of  small  ripples  which  inter¬ 
act  with  radar  waves  via  Bragg  resonance  scattering.  The  radiation 
or  action  balance  equation  (a  WKB  approximation)  describes  how  N  is 
transported  or  radiated  in  space,  time  and  wavenumber  (x,  t,  K) 
upon  Interaction  with  a  surface  current  u(x,  t)  or  other  perturba¬ 


tion,  to  wit 


[For  references  see  Keller  &  Wright  (1975),  Alpers  and  Hasselman 
(1978)  or  Wright  (1978).]  The  waves  propagate  in  four-dimensional 
phase  space  along  trajectories  described  by  ray  equations  as  follows 


_3n 
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9K 
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- 
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where 


(2-8) 


a  ( X ,  K,  t)  -  ft'(K)  +-K  •  u(x,  t) 

and  u  is  the  current  with  which  the  waves  of  action  density  N  are 
interacting.  S  is  the  source  term  describing  energy  input  to  and 
loss  from  the  waves  described  by  N. 

In  certain  special  cases,  S  ■  0  and  (2-7)  becomes  a  more 
simple  equation  expressing  the  conservation  of  action  spectral 
density.  In  this  case,  the  waves  propagate  without  energy  gain  or 
loss  through  the  slowly  varying  medium.  Phillips  (1981)  uses  this 
case  to  describe  the  interaction  of  small  A  -  A^  ripples  with  long 
ocean  waves.  The  interacting  current  u  is  the  tangential  surface 
velocity  of  the  long  wave  at  some  given  point.  Further  gravita¬ 
tional  acceleration  g  ♦  g'  where  g'  is  the  component  of  g  normal  to 
the  long  wave  surface  at  the  given  point  modified  by  the  accelera¬ 
tion  due  to  the  long  wave  orbital  motion. 


If  a  perturbation  approach  (based  on  the  ad  hoc  assumption 
that  the  action  tends  to  return  to  an  equilibrium  state)  is  used 
with  regard  to  N,  u,  and  S  (Alpers  and  Hassleman,  1978),  then  the 
radiation  balance  equation  (2-7)  can  be  approximated  as 

~  -  -w«N  (2-9) 

where  y  is  the  relaxation  rate  associated  with  a  relaxation  time 

■  u  1  and  AN  is  a  perturbation  in  action  spectral  density  with 
respect  to  a  constant  spectrum  Nq(K).  The  relaxation  time  is 
determined  by  the  totality  of  energy  gain  and  loss  processes,  such 
as  wave  generation  by  wind,  energy  exchange  by  resonant  wave-wave 
Interaction  and  energy  loss  by  viscous  damping  and  wave  breaking. 
Notwithstanding  a  lack  of  knowledge  of  the  physics  of  these 
processes,  is  thought  to  be  of  the  order  of  10  to  100  wave 
periods.  For  Bragg  resonant  waves  in  the  case  of  SEASAT  SAR 
-  30  cm  would  be  -  5  to  50  seconds.  No  experimental 
measurements  of  have  yet  been  made  in  the  open  ocean  to  the 
author's  knowledge.  However,  recent  work  by  Alpers  (1984)  suggests 
that  t  probably  lies  toward  the  upper  end  of  the  range. 

The  important  point  here  is  that  for  the  interaction  of  Bragg 


resonant  ripples  (A  -  A^)  with  currents  due  to  gravity  waves  of 


lengths  like  that  of  waves  in  a  Kelvin  wake  or  in  ocean  swell  the 
source  term  S  In  (2-7)  can  be  set  to  zero  with  some  confidence  and 
the  spectral  density  conservation  equation  used.  This  is  because 
the  current  associated  with  the  long  waves  varies  over  a  time  scale 
<  (P/4)  where  P  is  the  wave  period,  and  this  time  is  typically 
about  an  order  of  magnitude  less  than  t^.  Thus,  the  hydrodynamic 
Interaction  of  Bragg  resonant  waves  with  long  gravity  waves  is 
relatively  weak  since  the  currents  associated  with  the  long  gravity 
waves  do  not  act  in  a  given  direction  for  a  long  enough  time. 

Other  surface  currents  may  vary  relatively  slowly,  i.e.,  on 

time  scales  much  longer  than  T  .  Hence,  use  of  the  S  ♦  o  approxi- 
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mation  is  probably  not  justified  and  the  form  (2-9)  is  appropriate. 
Work  by  Alpers  (1984)  on  currents  related  to  underwater  typography 
in  shallow  depths  suggests  that  analysis  based  on  (2-9)  (see,  for 
example,  Alpers  and  Hasselmann,  1978)  would  predict  relatively 
strong  modulation  of  Bragg  resonant  ripples  and  hence  relatively 
strong  radar  cross  section  modulations  due  to  equally  strong 
currents  associated,  for  example,  with  shallow  water  internal 
waves.  Relatively  here  refers  to  the  case  of  gravity  waves  in  the 
paragraph  above. 


The  Important  point  here  la  that  for  gravity  wave  modulation 
of  Bragg  resonant  ripples  the  surface  currents  can  not  act  for  a 
time  comparable  to  Tr  and  so  S  can  be  set  to  zero  in  (2-7). 

However,  for  slowly  varying  surface  currents,  the  interaction  time 
is  probably  much  longer  and  the  use  of  a  source  term  is  probably 
appropriate,  e.g.,  equation  (2-9).  Experimental  work  suggested  in 
Section  2.6  below  could  help  sort  out  the  role  of  S  in  describing 
the  SAR  imaging  of  Internal  waves.  The  approach  to  S  described  here 
is  probably  the  simplist  aside  from  letting  S  go  to  zero.  Other 
treatments  exist  addressing  particular  source  mechanisms,  e.g.  wind 
energy  input  to  waves  (Hughes,  1978). 


2.5.2  Effects  Associated  with  Kelvin  Wakes 

2.5.2. 1  Elementary  Properties  of  Kelvin  Wakes.  Dashen 
et  al.  (1982)  give  a  concise  review  of  the  properties  of  Kelvin 
wakes.  Further  information  is  given  in  Section  4.0  of  this  report 
and  by  Lighthill  (1978)  and  references  therein.  To  make  matters 
clear,  we  reproduce  a  convenient  illustration  from  Dashen  et  al.  as 
Figure  2-12.  First,  we  remark  that  the  classical  solution  to  the 
Kelvin  wake  problem,  used  here  to  illustrate  several  points,  is  an 
idealized  model  of  a  real  ship  wake.  Hence,  we  expect  that 
elementary  calculations  made  here  will  have  to  be  modified 
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Figure  2-12.  First,  we  remark  that  the  classical  solution  to  the  Kelvin  wake  problem,  used 
here  to  illustrate  several  points,  is  an  idealized  model  of  a  real  ship  wake. 
Hence,  we  expect  that  elementary  calculations  made  here  will  have  to  be 
modified  significantly  in  order  to  apply  to  a  real  ship  wake.  Further,  we  Figure 
2-12:  SEASAT  SAR  Observes  a  ship’s  Kelvin  wake.  The  ship  is  moving  at  an 
angle  0  relative  to  a  line  parallel  to  the  SEASAT  velocity  vector.  The  SEASAT 
angle  of  incidence  is  8  '~23° .  Adapted  from  Dashen  et  al.  (1982). 


significantly  in  order  to  apply  to  a  real  ship  wake.  Further,  we 
note  that  many  ships  utilize  modified  bows  which  reduce  the  wave 
drag  on  the  ship,  and  hence  modify  the  Kelvin  wake  pattern;  for 
example,  reducing  the  wave  height  of  the  long  wavelength  waves  near 
the  cusps. 

Perhaps  the  most  Important  feature  of  the  Kelvin  wake  pattern 
is  the  similarity  of  the  wake  form  along  a  line  of  constant  6  (see 
Figure  2-12).  First,  the  orientation  (in  the  x,y  plane)  of  the 
divergent  or  transverse  wave  crests  (or  any  locus  of  constant  wave 
phase)  is  constant  along  constant  8  .  Second,  the  wavenumber  of 

the  waves  in  either  the  divergent  or  transverse  wakes  is  constant 
along  a  line  of  constant  8.  As  Dashen  et  al.  (1982)  point  out, 
these  properties  would  lead  to  a  V-shaped  wake  if  there  is  some 
selection  process  which  leads  to  a  local  maximum  of  oq  at  a  given 
wavenumber  K  or  crest  orientation. 


2 . 5 . 2 . 2  Qua8l-Specular  and  Bragg  Resonance  Backscatter. 

Bahar  et  al.  ( 1983a, b,c)  discuss  radio  wave  scattering  from  rough 
surfaces  which  is  relevant  to  the  ocean.  This  work  is  a  compre¬ 
hensive  treatment  of  rough  surface  scattering  in  that  the  formalism 
includes  both  the  Bragg  resonance  and  quasi-specular  mechanisms  is  a 
self-consistent  manner  as  well  as  shadowing  effects  important  at 


large  angles  of  Incidence .  Further ,  this  theory  claims  to  handle  a 
more  general  class  of  composite  rough  surfaces;  namely,  those  that 
cannot  be  decomposed  into  small-scale  surface  perturbations 
(amenable  to  Bragg  resonance  theory)  riding  on  a  filtered  large- 
scale  surface  (amenable  to  quasi-specular  theory).  A  draw  back  of 
this  formalism  is  that  the  mathematics  are  relatively  complicated 
and  hence  not  too  useful  in  analytic  work  as  they  stand.  However, 
suitable  approximations  can  be  made  and  reduce  some  expressions  to 
tractable  forms  useful  in  analytic  work.  For  example,  under  the 
appropriate  restrictions,  Bahar  et  al.'s  formalism  reduces  to 
equations  (2-2,  2-3  and  2-4)  in  Section  2.4  above. 

Bahar  et  al.  (1983c)  calculate  the  backscatter  cross  section 
for  arbitrarily  oriented  composite  rough  surface  regions  of  limited 
extent.  This  is  just  the  sort  of  information  one  needs  to  deal  with 
the  resolution  cell  sized  patches  observed  by  SAR  within  a  ship's 
Kelvin  wake. 

Bahar  et  al.  (1983a, b,c)  view  ocean  surface  height  fluctua¬ 
tions  c(x,y)  as  being  composed  of  two  components 

?(x,y)  *  ?1(x,y)  +  Cg(x,y)  (2-10) 
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where  1  and  s  refer  to  large  and  snail-scale  spectral  components 
The  mean  radar  cross  section  of  the  surface  Is  thus  given  by 


<  a 


PQ 


<  o 


PQ 


>!  +  < 


PQ 


(2-11) 


where  PQ  refers  to  the  received  and  transmitted  polarizations,  e.g., 
PQ  «  HH  for  SEASAT.  The  cross  section  from  the  large-scale 
components  is  given  as 

<  oPQ  >x  -  |x8  (v  •  ng) | 2  <  oPQ  >  .  (2-12) 

g 

X  is  the  characteristic  function  of  the  small  scale  surface 
structure  (where  for  backscatter  v  »  2  it  and  ng  is  the  unit 
vector  normal  to  the  surface  at  a  specular  point).  <  oP<^  >  is  the 
specular  point  cross  section  for  the  large  scale  surface.  The 
expression  for  <  oP<^  >  (Bahar  et  al.,  1983c,  p.  677)  is  somewhat 

O 

involved,  but  analogous  to  eq.  (2-4).  The  x  term  reduces  the 
PQ 

cross  section  <  o  >  to  account  for  the  fact  that  the  small-scale 

roughness  (C  in  2-10)  makes  the  specular  points  less  effective 
s 

scatters  (c.f.  section  2. 4. 1,2  and  Figure  2-8).  For  a  Gaussian 
distributed  rough  surface  ( 


corresponds  Co  Che  first- 


PO 

The  flrsC  Cerm  in  Chis  series  <  a  >  , , 

s  1 

order  Bragg  scaCCerlng  (c.f.  eq.  2-3). 

Alchough  Bahar  eC  al.  do  noC  calculaCe  which  Cerm  dominaCes  in 
(2-10)  under  various  condiCions,  Vesecky  eC  al.  (1983)  have  shown 
ChaC  for  surface  CilC  modulation  of  0°,  quasi-specular  scattering 
is  an  important  factor  at  Che  SEASAT  angle  of  incidence.  Alchough 
there  is  still  a  little  controversy  it  is  pretty  clear  that  scatter¬ 
ing  from  Bragg  resonant  ripples  will  be  perceived  by  the  radar  as 
having  a  doppler  shift  corresponding  to  the  small  (~  0.7  m/s) 
phase  velocity  of  the  Bragg  waves  themselves  plus  the  comparable 
orbital  velocity  of  the  large  waves  which  advect  the  small  waves 
(Graf  et  al.,  1977).  However,  it  is  not  clear  what  doppler  shift 
should  be  associated  with  the  quasi-specular  scattering  process. 
Hasselmann  et  al.  (1984)  suggest  that  for  quasi-specular  scatter  the 
long  wave  phase  velocity  is  important.  If  this  is  so,  wake  waves 
scattering  by  the  quasi-specular  process  would  have  substantial 
doppler  shifts  and  the  SAR  process  would  misplace  scattered  energy 
by  at  least  several  resolution  cells  along  the  azimuth  direction. 

We  note  that  in  most  cases  (e.g..  Figure  2-4)  the  ship  and  its  wake 
are  not  displaced  together.  This  argues  that  the  doppler  shift  of 
the  ship  and  its  wake  are  distinctly  different.  Further  our  calcu¬ 
lations  in  Section  2. 5. 2.3  show  that  the  size  of  the  displacement  of 


the  ship  from  its  wake  is  very  close  to  what  one  would  expect  from 
the  ship's  speed  alone .  This  implies  that  the  doppler  shift  of  the 
wake  (in  Figure  2-4)  is  virtually  stationary  with  respect  to  the 
mean  ocean  surface  and  that  the  phase  velocities  of  the  large 
(greater  than  twice  SAR  resolution)  waves  in  the  wake  are  not 
determining  the  doppler  shift  observed  by  the  SAR.  Hence  for  Figure 
2-4  we  conclude  that  quasi-specular  scattering  plays  at  most  a  minor 
role.  However,  in  Figure  2-7  there  are  very  bright  linear  features 
in  some  wakes  which  appear  to  end  on  the  ship  (rather  than  being 
displaced  as  is  often  the  case).  This  suggests  that  these  bright 
features  may  be  scattering  via  the  quasi-specular  mechanism,  and 
therefore  have  a  doppler  shift  comparable  to  the  ship  speed. 

While  large  scale  surface  height  fluctuations  (K  <  K  ) 

s&r 

would  be  associated  with  the  large  scale  waves  of  a  Kelvin  ship 
wake  (6  not  near  zero  in  Figure  2-12);  small-scale  roughness 

(K  >  K  )  must  be  generated  either  by  an  ambient  wind  or  by  the 
emn 

ship  itself  (see  Figure  2-9).  If  an  ambient  wind  is  involved,  then 
the  wake  would  have  small  scale  structure  approximately  uniformly 
distributed  over  the  wake,  but  subject  to  modification  by  large 
scale  roughness  via  tilting  and  wave-wave  interaction  as  discussed 
below.  Also,  both  wake  features  and  the  background  on  which  the 
wake  is  superimposed  are  similarly  enhanced  as  wind  speed  Increases. 


While  wind  direction  does  not  seem  very  important  at  the  L-band 
frequency  of  SEASAT  it  is  more  so  at  higher  radar  frequencies.  This 
comment  is  based  on  the  known  wind  direction  dependence  of  ocean 
radar  cross  section  at  various  frequencies  (Jones  &  Schroeder, 

1978). 

Small-scale  roughness  is  also  associated  with  the  Kelvin  wake 
itself  for  regions  close  to  the  velocity  axis  of  the  ship. 

Classical  Kelvin  wake  theory  predicts  small  scale  roughness  for 
small  0  (see  Figure  2-12)  in  that  the  waves  near  the  ship  velocity 
axis  must  travel  slowly,  and  therefore  via  their  dispersion  relation 
have  small  wavelengths.  Dashen  et  al.  (1982)  discuss  this 
phenomena.  In  addition  small  scale  roughness  may  be  generated  by 
wave  breaking  and  other  processes  which  occur  when  a  real  ship  of 
finite  size  generates  a  wake.  For  example  generation  of  small  waves 
by  a  real  ship  probably  occurs  non-uniformly  (in  both  space  and 
time)  over  the  length  of  the  ship  (and  possibly  behind  it). 

Finally,  the  quasi-specular  process  being  related  to  large 
scale  roughness  does  not  require  generation  of  small  scale  waves  as 
does  the  Bragg  mechanism.  Hence  for  large  scale  (>  20A)  waves 
generated  in  a  ship's  Kelvin  wake  quasi-specular  backscatter  can 
occur  whether  or  not  a  wind  is  blowing. 
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2. 5. 2. 3  Tilt  Mechanism.  SEASAT  SAR  observed  the  ocean  at  an 

0 

angle  of  incidence  0  ~  23  with  horizontal  polarization.  These 

parameters  are  particularly  well  suited  to  distinguishing  resolution 

cell  sized  patches  tilted  at  different  angles  relative  to  the  radar 

line  of  sight.  The  tilts  here  refer  to  the  tilting  introduced  by 

the  waves  of  the  Kelvin  wake.  Figure  2-13  shows  calculations  by 

Bahar  et  al.  (1983c)  of  fractional  changes  in  radar  cross  section 

do/<o>  with  respect  to  tilts  in  (dft)  and  perpendicular  to  (dr) 

the  plane  of  incidence.  Although  these  curves  are  for  X-band 

(X  -  3  cm)  ,  not  L-band  (X  ■  23  cm)  as  used  by  SEASAT  SAR,  and 

use  a  particular  sea  surface  model  (see  Bahar  et  al.,  1983c, 

p.  626),  we  can  note  general  features  using  a  similar  calculation  by 

Vesecky  et  al.  (1983)  at  L-band  for  reference.  First,  tilting  in 

the  plane  of  incidence  (with  respect  to  ft  )  is  much  more  effective 
HH 

in  modulating  <  a  >  than  is  tilting  perpendicular  to  the  plane 
of  incidence  (with  respect  to  t  ).  Thus  so  far  as  tilting  induced 
by  waves  in  the  Kelvin  wake  is  concerned,  waves  in  the  divergent 
portion  of  the  wake  should  be  more  prominent  when  the  ship  is 
traveling  nearly  along  the  azimuth  direction  ($  near  0  or  180°  in 
Figure  2-12).  Waves  in  the  transverse  portion  of  the  wake  should  be 
more  visible  when  the  ship  is  traveling  near  the  range  direction 
(<fi  near  90°  or  270°  in  Figure  2-12).  Referring  to  Section  2. 5. 2. 3 
above,  we  see  from  Figure  2-4  that  indeed  the  transverse  waves  are 
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Figure  2-13.  Rate  of  fractional  change  of  radar  backscatter  cross  section  (horizontal  transmit 
and  receive  polarization)  with  respect  to  tilts  of  the  resolution  cell  (a)  in  the  plane 
of  incidence  (d  <  crHH  >  /  dfl)  /  <  crHH  >  and  (b)  perpendicular  to  the  plane 
of  incidence  (d  <  cthh  >1  &t)  I  <  crHH  >).  Triangles,  octagons  and  squares 
refer  to  resolution  cell  sizes  of  3  m,  10  m  and  25  m  (SEASAT),  respectively. 
e'0  is  the  angle  of  incidence  in  the  absence  of  any  tilt.  Radar  wavelength  is  \ 
=  3  cm  (X-band).  After  Bahar  et  al.  (1983c). 


more  prominent  in  this  case  when  $  »  -112°.  Further  the  trans¬ 
verse  waves  in  the  upper  portion  (above  the  dark  turbulent  wake)  of 
Figure  2-4  are  more  visible  than  in  the  lower  portion  of  the  wake. 
And  indeed  $  is  closer  to  90°  in  the  upper  portion  of  the  wake. 

As  a  crude  estimate  of  the  expected  peak  to  trough  cross  section 
variation  for  the  transverse  wake  waves  in  Figure  2-4,  we  let  the 
transverse  waves  be  sinusoidal  with  a  wavelength  A  *  94  m  (see 
Section  2. 5. 2. 3  above)  and  amplitude  A  »  0.5  m  (1  m  peak  to  trough). 
The  slope  thus  varies  by  +_  AK  *  _+  0.034  radians  *  _+  1.9°  around 
0^  -  23.8°  ,  the  mean  angle  of  incidence  for  this  image.  Referring 

iTU  uu 

to  Figure  2-13(a)  [(d  <o  >/d£l)/<  o  >]  *  12  and 

6  <am>  -  +  AK  [(d<oIIH>/dO)/<oHH>l  -  +  (0.5)(2ir/93)(  12)  -  +  0.4 
which  implies  a  maximum  to  minimum  cross  section  variation  of 
"  2.5  dB.  One  would  have  to  do  a  quantitative  calculation  using 
the  digital  version  of  Figure  2-4  to  check  this  figure,  but  qualita¬ 
tively  the  expected  and  observed  cross  section  variations  are  small. 

Another  implication  of  Figure  2-13  is  that  since  the  angle  of 

incidence  for  SEASAT  was  9^  ~  23®,  SEASAT  operated  near  the 

HH 

combination  of  maximum  mean  <?  >  (see  Figure  2-8)  and  maximum 

HH 

sensitivity  of  6  <o  >  to  changes  in  surface  slope.  Note  that 
Figure  2-13  refers  to  a  tilted  resolution  cell;  hence  we  have 
implicitly  assumed  that  any  wave  slopes  referred  to  are  for  waves 


with  wavelength  at  least  two  resolution  cells  in  size,  i.e., 

A  >  50  n  for  SEASAT. 

In  some  cases,  e.g.,  Figure  2-4,  a  wake  feature  associated 
with  the  line  of  cusps  at  0  =■  +  19.5°  in  Figure  2-12  can  be 
identified.  The  wave  shape  for  a  Kelvin  wake  is  unique  along  the 
cusp  line  in  that  first,  the  wave  height  is  generally  largest  there 
and  second^,  the  wave  decays  in  nearly  every  direction  away  from  the 
crest.  In  terms  of  the  tilt  mechanism  one  would  expect  enhanced 
radar  backscatter  if  the  radius  of  curvature  of  the  "boss"  shaped 
features  at  the  cusp  were  large  enough  (much  larger  than  a  resolu¬ 
tion  cell  size)  and  properly  oriented.  Otherwise,  one  could  expect 
reduced  radar  cross  section  since  the  boss  would  either  have  only 
small  portion  tilted  more  nearly  (than  the  mean  surface)  toward  the 
radar  and  the  rest  tilted  less  nearly  toward  the  radar.  This  second 
case  appears  to  obtain  in  Figure  2-4,  since  the  cusp  line  is  dark  in 
the  SAR  image.  However,  there  are  in  other  cases  bright  cusp  lines, 
e.g.,  Figure  22  of  Shuchman  et  al.  (1983)  where  the  X-band  aircraft 
SAR  resolution  cell  is  much  smaller  than  SEASAT's  25  m. 

In  the  above  discussion  we  have  not  included  the  ambient  ocean 
waves  which  may  be  present  and  of  higher  amplitude  than  waves  in 
Kelvin  wakes.  In  general,  a  spectrum  of  ambient  ocean  waves  large 
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compared  to  the  Kelvin  wake  should  tend  to  obscure  the  Kelvin  wake 
since  they  increase  the  mean  ocean  radar  cross  section ,  and  hence 
this  cross  section's  variance  from  one  resolution  cell  to  the  next — 
cross  section  variance  is  approximately  equal  to  the  mean  for  one- 
look  SAR  processing. 

2. 5. 2. 4  Hydrodynamic  Wave-Wave  Interaction.  Both  experi¬ 
mental  and  theoretical  studies  indicate  that  hydrodynamic  wave-wave 
interaction  causes  short  waves  (wavelengths  ~  A^)  riding  on  longer 
waves  to  be  strained  by  the  long  waves  such  that  the  amplitude  of 
the  short  waves  is  enhanced  somewhat  in  advance  of  the  long  wave 
crest  at  the  position  where  a  surfer  would  ride.  For  waves 
traveling  toward  the  radar  this  position  coincides  approximately 
with  the  portion  of  the  long  wave  tilted  toward  the  radar  K 
vector.  Hence,  the  tilt  and  hydrodynamic  mechanisms  would  both 
enhance  the  radar  cross  section  at  nearly  the  same  portion  of  the 
long  wave  for  waves  traveling  toward  the  radar.  However,  for  waves 
traveling  away  from  the  radar  the  two  effects  would  tend  to  cancel. 
Overall  one  would  expect  to  find  wake  waves  traveling  toward  the 
radar  to  be  more  prominent  in  SAR  images  than  waves  traveling  away, 
provided  the  tilt  id  hydrodynamic  mechanisms  are  about  equally 
effective  in  modulating  the  radar  cross  section  wake  waves.  Because 
of  this  effect  one  would  usually  expect  one  arm  of  a  wave  image  to 


be  more  prominent  than  the  other.  There  are  many  examples  of  wakes 
with  only  one  arm  visible  in  the  SAR  image.  However,  there  are 
images  in  which  two  (or  more)  arms  of  a  V-shaped  wake  are  visible, 
e.g..  Figures  2-4,  2-5  and  2-6.  Shuchman  et  al.  (1983,  p.  33)  find 
that,  "Usually,  one  side  of  the  Kelvin  wake  envelope  was  more 
visible  than  the  other  side." 

2. 5. 2. 5.  Surface  Motion  Doppler  Effects — Velocity  Bunching 
and  Resolution  Smearing.  Motion  of  certain  local  portions  of  the 
ocean  surface  with  respect  to  the  motion  of  the  large  scale  mean 
surface  produces  two  major  effects.  First,  if  the  motion  in  a 
resolution  cell  sized  patch  of  surface  is  largely  coherent,  i.e., 
there  is  a  dominant  velocity  in  terms  of  both  space  and  time  on 
scales  comparable  to  the  SAR  resolution  cell  size  and  the  SAR 
integration  time,  then  the  radar  echo  energy  associated  with  this 
dominant  velocity  will  be  doppler  shifted  differently  than  the  mean 
large  scale  surface  and  the  echo  energy  will  be  shifted  along  the 
azimuth  direction  to  an  incorrect  location.  The  ship  in  Figure  2-4 
illustrates  this  point  well.  Second ,  if  the  motion  in  a  resolution 
cell  sized  patch  is  spread  over  a  range  of  velocities,  i.e.,  there 
is  no  strong  dominant  velocity  and  a  range  of  velocities  exists  over 
the  aforementioned  spatial  and  temporal  scales,  then  the  radar  echo 


energy  is  spread  in  the  azimuth  direction  over  a  number  of  SAR 


resolution  cells.  This  smearing  makes  long  ocean  waves  traveling 
nearly  along  the  azimuth  direction  difficult  to  image  (viz.  Beal  et 
al.,  1983  and  Vesecky  et  al.,  1984). 

As  Alpers  and  Rufenach  (1979)  and  Alpers  et  al.  (1981)  have 
pointed  out  coherent  and  systematic  spatial  patterns  of  surface 
motion  can  cause  ocean  surface  phenomena  to  be  imaged  by  a  SAR  (see 
Section  2.4.2).  For  example,  the  periodic  variations  of  surface 
motion  related  to  long  ocean  waves.  This  mechanism  is  most 
effective  for  gravity  waves  traveling  near  the  azimuth  direction 
(y"-  direction  in  Figures  2-2  and  2-12)  and  is  least  effective  for 
gravity  waves  traveling  in  the  range  direction  (x'-  direction  in 
Figures  2-2  and  2-12).  In  connection  with  natural  ocean  gravity 
waves  there  has  been  considerable  controversy  over  the  effectiveness 
of  the  velocity  bunching  mechanism.  Both  Harger  (1984)  and  Jain 
(1981)  argue  that  for  natural  ocean  gravity  waves  no  velocity 
bunching  modulation  should  occur. 

As  analyzed  by  Alpers  et  al.,  velocity  bunching  is  in  general 
a  nonlinear  process.  The  parameter  which  governs  the  nonlinearity 
is  C  where  the  change  from  linear  to  nonlinear  behavior  occurs  for 


C  >  0.3  where 


C  -  (R/V)  ftA  (sin  0'  sin  *  +  cos^  S’)'2  cos  *  (2-14) 

and  R  is  range  to  target,  V  is  SAR  platform  velocity,  A  and  ft  are 
wave  amplitude  and  frequency  and  6'  and  *  are  given  in  Figure  2-2. 
Since  R/V  is  relatively  large  for  SEASAT,  namely  128,  one  expects 
nonlinear  behavior  to  set  in  for  rather  small  values  of  wave  ampli¬ 
tude.  Thus  for  ship  wakes  where  wave  amplitude  may  be  low,  but 
ft  is  relatively  high  nonlinear  behavior  is  very  likely. 

Since  velocity  bunching  relies  on  spatial  modulation  of 
surface  velocity  which  is  on  a  scale  larger  than  the  SAR  resolution 
cell,  resolution  cells  which  contain  more  than  one  wavelength  of  the 
waves  making  up  the  Kelvin  wake  at  a  given  location  can  not  show 
much  velocity  bunching  modulation  (from  one  resolution  cell  to  the 
next).  That  is,  resolution  cell  sizes  greater  than  the  spatial 
scale  of  orbital  velocity  variations  effectively  average  out  any 
velocity  bunching  modulation  present.  We  discuss  how  this  trans¬ 
lates  into  wake  effects  in  Section  2. 5. 2. 6  below. 

An  effect  related  to  the  resolution  cell  averaging  above  is 


degradation  of  azimuth  resolution  by  surface  velocity  variations 
within  a  resolution  cell.  These  velocity  variations  effectively 
smear  out  (along  the  azimuth  direction).  The  backscattered  energy 


which  should  be  associated  with  a  single  resolution  over  two  or  more 
resolution  cells.  Since  any  resolution  cell  sized  patch  of  ocean 
surface  contains  a  statistical  distribution  of  ambient  ocean  waves 
in  addition  to  any  wake  related  waves  a  statistical  distribution  of 
surface  velocities  will  be  present  to  degrade  azimuth  resolution. 
Similarly  Beal  et  al.  (1983)  model  this  effect  as  a  gausian-shaped 
impulse  response  function 

hv(y)  -  exp  (-y2/2  £^  )  (2-15) 

o 

where  £  ■  (R/V)  £  and  £  is  the  variance  of  the  radial  (with 

yv  v  v 

respect  to  the  radar)  velocity  distribution  averaged  over  a  SAR 
resolution  cell.  Similarly  if  there  is  acceleration  in  a  coherent 
(over  the  resolution  cell)  surface  motion,  degradation  of  resolution 
along  the  azimuth  direction  results.  Alpers  et  al.  (1981)  give  this 
degradation  in  azimuth  resolution  as  a  function  of  undegraded  reso¬ 
lution,  SAR  Integration  time,  number  of  looks,  long  wave  amplitude 
and  frequency  and  other  observational  parameters  such  as  (R/V).  Beal 
et  al.  (1983)  and  Alpers  et  al.  (1981)  give  different  expressions 
combining  both  effects.  The  principal  point  here  is  that  resolution 
cell  degradation  can  effectively  average  out  any  velocity  bunching 
(or  other)  modulation  mechanisms  by  enlarging  the  resolution  cell  to 
a  size  where  waves  traveling  along  the  azimuth  direction  can  not  be 


resolved,  l.e.,  the  Nyquist  sampling  criterion  cannot  be  satisfied. 


Clearly  waves  observed  on  a  calm  day  [low  E^  in  equation  (2-15)] 


would  be  most  likely  to  be  imaged  via  the  velocity  bunching 
mechanism.  Indeed  ship  wakes  are  more  frequently  imaged  on  calmer 
days.  However,  ship  wakes  are  more  frequently  imaged  when  the  wake 
waves  are  traveling  near  the  range  direction  where  the  velocity 
bunching  mechanism  is  weakest. 


2. 5. 2. 6  Effects  of  Finite  Resolution  Cell  Size.  Although  one 
usually  thinks  of  SAR  as  creating  a  map  of  the  surface  distribution 
of  radar  cross  section,  the  process  of  making  the  map,  i.e., 
achieving  high  spatial  resolution,  actually  alters  the  cross  section 
observed.  Further  the  finite  resolution  of  a  SAR  causes  variations 
in  radar  cross  section  on  scales  smaller  than  the  resolution  cell 
size  to  be  averaged. 


With  regard  to  changes  in  cross  section  we  refer  to  material 
in  Section  2.4  and  to  Figure  2-9.  We  note  that  in  averaging  Bragg 
resonant  scattering  over  the  intermediate  scale  wavelengths  the 
resulting  cross  section  is  determined  by  an  integral  in  which  Kgar 
is  the  lower  limit.  Clearly  then  this  integral  is  dependent  on  its 
lower  limit  and  thus  on  resolution  cell  size.  Work  by  Vesecky 
et  al.  (1983)  illustrates  these  points. 


Most  of  the  discussion  above  concerns  expected  radar  cross 
section  variations  from  one  portion  (crest,  trough,  etc.)  to  another 
of  a  surface  gravity  wave.  Obviously,  it  is  of  interest  to  know  how 
many  wavelengths  of  Kelvin  wake  waves  there  are  in  a  SAR  resolution 
cell  for  different  locations  in  a  Kelvin  wake.  Because  if  there  is 
more  than  about  one  wavelength  within  a  SAR  resolution  cell  the 
variations  between  different  portions  of  wave  will  be  averaged  out 
and  little  variations  (from  the  mechanisms  discussed  above)  from  one 
resolution  cell  to  the  next.  This  does  not  mean  that  difference 
between  one  resolution  cell  and  the  next  might  not  occur  because  of 
changes  in  wave  or  observational  parameters  from  one  resolution  cell 
to  the  next — we  discuss  this  notion  below.  In  Figure  2-14,  we 
display  computations  of  the  number  of  wavelengths  (of  the  local 
Kelvin  wake  waves)  per  SAR  resolution  cell.  The  Kelvin  wake  wave¬ 
lengths  are  computed  using  the  formulas  given  by  Dashen  et  al. 
(1982).  The  Important  parameters  are  angle  9  in  figure  2-12,  ship 
speed  U  and  SAR  resolution.  We  consider  only  the  divergent  waves  in 
the  wake  since  most  of  the  ship  wakes  of  Interest  (narrow  Vees  as  in 
Figure  2-5)  apparently  Involve  the  divergent  waves  rather  than  the 
transverse  ones.  It  is  interesting  to  note  that  the  angle  6  at 


which  SAR  resolution  equals  Kelvin  wake  wavelength  decreases  with 
increasing  ship  speed  from  about  17°  for  Ug^^p  ■  10  m/s  to  about  12° 
at  15  m/s  to  9°  at  20  m/s.  There  is  a  suggestion  that  narrower 


Wavelengths/SAR  Cell  Wavelengths/SAR  Cell 


wakes  in  terms  of  6  are  associated  with  faster  ships,  for  example. 
Figure  2-5  is  thought  to  be  a  fast  bulk  carrier  while  Figure  2-6  is 
thought  to  be  a  fleet  of  relatively  slow  fishing  ships.  Neverthe¬ 
less,  from  these  calculations  one  would  not  expect  to  see  individual 
wake  waves  imaged  for  6  less  than  9°  and  there  are  clearly  wake 
images  associated  with  0  -  3°  or  4°  . 

2. 5. 2. 7  An  Hypothesis  for  SAR  Imaging  of  Kelvin  Wake  Waves 
for  9  ~  10°  .  To  explain  the  fine  linear  features  in  ship  wakes 
illustrated  in  Section  2.3  one  must  find  mechanisms  to  enhance  the 
radar  return  at  one  or  more  singular  values  of  0.  The  hypothesis 
described  here  makes  use  of  the  tilt,  hydrodynamic  and  velocity 
bunching  mechanisms  which  are  thought  to  make  natural  ocean  surface 
waves  visible  in  SAR  Images.  As  we  have  just  argued  in  Section 
2. 5. 2. 6  above,  this  hypothesis  can  not  apply  to  the  region  of  the 
wake  where  0  <  10°  because  the  wavelengths  of  the  Kelvin  wake 
waves  are  too  small,  i.e.,  smaller  than  a  SEASAT  resolution  cell. 
This  hypothesis  argues  that  for  a  given  wake  and  observation  situa¬ 
tion  a  particular  wavenumber  (K*)  in  the  Kelvin  wake  will  cause  a 
constructive  combination  of  one  or  more  of  the  above  mechanisms. 


This  particular  wavenumber  corresponds  to  a  particular  0  -  0  *  and 
hence  to  a  bright  line  in  the  wake's  SAR  image.  At  wavenumbers  near 
but  not  at  K*  the  tilt,  hydrodynamic  and  velocity  bunching 


mechanisms  do  not  operate  "in  phase"  and  hence  appear  much  like  the 
ocean  clutter  background.  The  constructive  combination  occurs  over 
only  a  small  range  of  6  because  ocean  wavenumber  in  the  wake 
changes  rapidly  with  9.  The  constructive  combination  occurs  at 
different  values  of  9  for  different  cases  because  the  wake  itself 
as  well  as  observational  parameters  change,  e.g.,  aspect  angle  ♦ 
and  incidence  angle  9'.  Figure  2-15  illustrates  that  the  strength 
of  the  tilt  mechanism  is  a  strong  function  of  angle  of  incidence 
9'  especially  near  the  SEASAT  observational  angles  of  incidence 
0i  m  190-250  . 

Some  wakes,  such  as  Figures  2-6  and  2-7,  may  be  explained  by 
this  hypothesis.  Even  some  wakes  having  bright  lines  meeting  at 
narrow  angles  may  be  explained  by  arguing  that  there  are  two 
singular  values  of  9  *,  namely  9^  *,  and  92*  on  the  same  side  of 
the  ship  velocity  vector  and  greater  than  -  10°  .  This  could 
produce  two  bright  linear  features  meeting  at  a  small  angle 

6°).  The  wake  at  upper  left  of  Figure  2-7  appears  to  corre¬ 
spond  to  this  situation. 

Very  bright  linear  features  could  be  produced  by  the  tilt 
mechanism  if  the  tilts  were  large  enough  to  cause  a  significant 
specular  or  quasi-specular  echo.  Except  for  short  wavelength  waves 


Figure  2-15.  Fractional  variation  of  radar  cross  section  with  tilt 

in  the  plane  of  incidence  —  (d  <  >  /  dft)  /  <  oHH  >  for 

different  angles  of  incidence  6' .  The  triangles,  circles  and 
squares  refer  to  SAR  resolution  cells  3  m,  10  m  and  25  m  in  size 
respectively-  The  radar  frequency  is  X-band,  A  =  3  cm.  After  Bahar 
et  al.  (1983c). 


(Section  2. 5. 2. 8)  it  is  unlikely  that  the  tilt  of  a  wake  wave  would 
present  a  specular  reflection  to  SEASAT — the  tilt  would  have  to  be 
at  least  some  19°.  However,  significant  quasi-specular  echo  can  be 
produced  at  tilts  of  much  less  than  19°  because  of  ambient  ocean 
waves  combining  with  the  wake  waves.  If  the  quasi-specular  mechan¬ 
ism  were  responsible  the  doppler  shift  would  be  near  the  phase 
velocity  of  the  wake  waves  and  hence  a  quasi-specular  wake  would 
tend  to  be  attached  to  the  ship  rather  than  shifted  as  in  Figure 
2-5.  Breaking  or  sharply  cu&ped  waves  could  produce  radar  back- 
scatter  similar  in  some  respects  to  the  quasi-specular  mechanism. 

If  one  examines  the  very  bright  wake  features  at  bottom  and  top  left 
of  Figure  2-7,  there  are  examples  of  which  may  be  explained  by  this 
hypothesis . 

Since  the  wake  waves  discussed  here  are  of  lengths  greater 
than  25  m,  they  can  persist  for  long  enough  times  that  a  long  wake 
feature  is  possible.  At  25  m  the  period  is  4  s  and  a  ship  at  10  m/s 
goes  40  m  in  4  s.  If  the  wave  persists  for  several  hundred  periods 
corresponding  to  several  e-folding  decay  times,  a  wake  feature 
dependent  on  these  waves  could  be  several  times  4  km  long.  Small 
Bragg  waves  needed  for  Bragg  resonance  scatter  would  be  generated 
continuously  by  the  wind. 
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This  hypothesis  is  presently  only  a  sketchy  notion.  It  needs 
further  development  using  empirical  and  theoretical  models  of  ship 
wakes.  At  present  the  scattering  theory  and  SAR  imaging  mechanism 
theory  is  enough  in  hand  to  make  appropriate  calculations  given  a 
wake  model. 

2. 5. 2. 8  Narrow  Wakes  at  9  ~  10°  .  For  Kelvin  wakes  at 
8  <  10°  the  waves  Involved  are  smaller  than  the  SAR  resolution 
cell  size.  To  explain  wake  features  here  one  must  rely  on  changes 
in  scattering  properties  (of  the  sea  surface)  which  are  averaged 
over  one  or  more  wavelengths  of  the  wake  waves  at  a  given  8. 

Hence,  the  wave  imaging  mechanisms  discussed  above  are  not  directly 
relevant — crests  and  troughs  lie  in  the  same  resolution  cell.  What 
is  required  here  is  a  local  (in  terms  of  8  )  enhancement  of  the 
radar  cross  section  averaged  over  the  numerous  wake  waves  in  a  SAR 
resolution  cell.  For  the  first  order  Bragg  mechanism  this  requires 
an  enhancement  of  wake  waves  at  A  ■  A^  and  propagating  along  the 
radar  range  direction.  And  Indeed  the  simple  Kelvin  wake  model 
provides  this  local  enhancement  since  surface  gravity  waves  are 
dispersive.  Waves  at  longer  wavelengths  (and  in  directions  other 
than  along  the  radar  range  direction)  might  be  involved  if  higher 
order  Bragg  scattering  were  Involved  (see  Johnstone,  1975  for 
reference).  For  the  quasl-specular  mechanism  a  statistical  ensemble 


of  wave  slopes  with  a  high  variance  is  needed  so  that  a  significant 

number  of  specular  points  will  appear  in  a  resolution  cell.  The 

short  wavelength  of  the  wake  waves  in  this  region  (6  <  10°) 

contributes  to  a  relatively  high  slope  variance  for  a  given  wave- 

height.  However,  the  short  wavelength  decreases  the  local  radius  of 

curvature  and  hence  lowers  the  radar  cross  section  (see  eq.  2-4). 

HH 

Whether  these  two  contending  effects  raise  or  lower  <  o  >  depends 
on  the  details  of  the  wake  in  this  region. 


In  order  to  make  a  first  order  estimate  of  the  radar  cross 
section  (using  eqs.  2-3  and  2-4)  in  this  portion  of  a  Kelvin  wake 
one  needs  to  know  the  basic  observational  parameters,  the 
directional  waveheight  spectrum  at  K  -  K.  along  the  range 
direction,  and  the  total  waveslope  variance — all  as  functions  of 
0.  To  make  a  better  estimate  the  complete  directional  waveheight 
spectrum  from  Ksar  upwards  to  should  be  known.  Breaking  waves  if 
they  exist  in  sufficient  quantity  could  be  important  for  narrow 
wakes.  Recall  that  if  there  is  sufficient  random  surface  motion  in 
a  resolution  cell  the  azimuth  resolution  of  the  SAR  will  be  degraded 
and  a  smearing  of  the  image  along  the  azimuth  direction  will  occur. 
Since  most  narrow  wakes  are  linear  features  along  the  azimuth  direc¬ 
tion,  this  smearing  would  tend  to  accentuate  features  along  this 
direction.  Azimuthal  smearing  would  be  enhanced  if  the  scattering 


features  were  exhibited  doppler  shifts  associated  with  wave  phase 
velocities  rather  than  wave  orbital  velocities.  Breaking  waves 
would  tend  to  travel  at  wave  phase  velocities  and  hence  would  be 
smeared  in  azimuth.  Evidence  of  this  phenomena  is  given  by  Shuchman 
et  al.  (1984).  Whether  one  might  expect  enhanced  wave  breaking  in 
the  interior  portions  of  a  Kelvin  wake  is  at  present  an  open 
question. 

Relevant  wake  hydrodynamics  are  discussed  elsewhere  in  this 

o 

report.  Several  factors  which  may  be  important  in  the  o  <  10 

region  are  mentioned  below.  The  turbulent  wake  smoothes  the  water 

and  reduces  <  o  >  for  very  small  values  of  0  .  Short  wavelength 
o  o 

waves  at  1  <  a  <  10  could  be  generated  by  mechanisms  other  than 

the  Idealized  Kelvin  wake  point  source;  for  example,  by  breaking 
waves  along  the  side  of  the  ship,  by  interaction  of  turbulent  wake 
related  currents  with  ambient  internal  waves  and  or  existing  short 
wavelength  waves,  and  by  wave-wave  interaction  between  the  short 
waves  of  the  divergent  wake  and  the  relatively  long  waves  of  the 
transverse  wake.  The  lifetime  of  Bragg  resonant  waves  can  be  short 
relative  to  the  length  of  the  observed  wakes  (  v  10  ki  ). 


2.5.3  Effects  Associated  with  Internal  Waves 


Although  the  hydrodynamics  of  surface  wave  and  Internal  wave 
wakes  are  obviously  very  different  (see  other  sections  of  this 
report),  there  are  many  scattering  and  radar  effects  which  are 
analogous  with  those  discussed  above  in  connection  with  surface  wave 
wakes . 


2. 5. 3.1  Relevant  Scattering  Mechanisms.  The  scattering 
mechanisms  Involved  In  Internal  wave  wakes  are  the  same  as  those 
discussed  In  Section  2.4.1  above.  One  relevant  difference  which  we 
discuss  below  is  that  Internal  waves  do  not  tilt  the  surface.  Hence 
quasi-specular  scattering  should  play  a  smaller  role  than  in  connec¬ 
tion  with  surface  wave  wakes. 

2. 5. 3. 2  Mechanisms  for  Spatial  Modulation  of  Radar  Cross 
Section.  There  are  several  means  by  which  internal  wave  surface 
currents  may  influence  radar  cross  section,  namely: 

a.  Modulation  of  local  waveheight  and  slope 
(Including  breaking  waves)  by  wave-current 


interaction 


b 


Current-induced  concentration  of  surfactants 


which  damp  existing  waves. 

c.  Current -induced  concentration  of  surface  debris. 

d.  Velocity  bunching  due  to  surface  currents. 

Each  of  these  items  is  discussed  briefly  below.  It  is  Important  to 
note  that  any  of  these  mechanisms  might  conceivably  dominate  under 
some  circumstances.  Item  (a)  is  the  mechanism  most  often  discussed 
in  connection  with  internal  wave  wakes.  However,  surfactants  can 
strongly  Influence  ripples  which  are  Involved  in  Bragg  scattering 
and  surface  currents  imply  velocity  bunching  (see  Section  2. 5. 2. 5). 

Increased  wave  height  and  slope.  As  discussed  in 
Section  2.4.3  above,  surface  currents  can  spatially  modulate  the 
ripples  and  longer  waves  involved  in  both  Bragg  and  quasl-specular 
scattering.  If  an  internal  wave  packet  could  sufficiently  enhance 
existing  surface  waves  at  current  convergence  zones,  then  clearly 
the  Internal  wave  wake  would  be  imaged.  The  degree  of  enhancement 
depends  on  the  Induced  currents  and  the  wave  current  interaction 
discussed  briefly  in  Section  2.4.3  and  at  length  in  Section  3.0 
below.  Two  points  of  interest  here  are  that  some  of  the  other 


mechanisms  discussed  immediately  below  may  be  involved  and  that  a 
source  term  may  need  to  be  involved  in  calculating  the  wave-current 
interaction  (see  Section  2.4.3). 

Radar  wave  scattering  from  breaking  waves  is  not  well  under¬ 
stood.  However,  as  discussed  in  Section  2. 5. 2. 8  above,  breaking 
waves  do  scatter  strongly  and  the  SAR  image  is  smeared  along  the 
azimuth  direction.  If  the  internal  wave  current  convergence  were 
strong  enough,  an  increase  in  the  area  density  of  breaking  waves 
would  presumably  occur.  The  increased  density  of  strongly  scatter¬ 
ing  breaking  waves  would  be  particularly  evident  for  a  current 
convergence  zone  oriented  nearly  along  the  azimuth  direction  because 
of  the  concurrent  azimuthal  smearing. 

Concentration  of  Surfactants.  Current  convergence 
zones  associated  with  internal  waves  would  presumably  concentrate 
surfactants  as  well  as  concentrating  the  spectral  density  of 
existing  ripples.  Clearly  these  two  effects  would  tend  to  counter¬ 
act  each  other  if  they  were  "in  phase,”  i.e.,  active  at  the  same 
location  of  the  internal  wave  as  one  expects  they  should  be.  Work 
by  Huhnerfuss  (1983)  shows  that  surfactants  do  indeed  show  up  as 
dark  patches  on  SAR  images.  Further  there  is  some  evidence  that  in 
the  tropical  Atlantic  ocean  surfactants  of  natural  origin  are  almost 


always  present  (Huhnerfuss,  et  al.,  1977).  Work  by  Bagg  (1983) 
indicates  that  there  is  no  pronounced  latitudinal  dependence  (from 
40°N  to  66°N)  in  the  frequency  of  observation  of  internal  waves  by 
SEASAT  SAR. 

Concentration  of  Surface  Debris.  In  some  waters  particularly 
along  coasts  there  is  a  significant  amount  of  surface  debris.  For 
example  debris,  such  as  foam,  flotsam  and  natural  floating  matter, 
collect  along  the  boundaries  of  the  Gulf  Stream.  If  such  matters 
were  concentrated  by  Internal  waves  in  a  ship  wake  if  could 
conceivably  cause  enhanced  backscatter.  At  this  point  such  a 
mechanism  is  speculative. 

Velocity  Bunching.  Since  surface  currents  associated  with 
Internal  waves  Impart  a  Doppler  shift  to  a  radar  echo,  the  different 
portions  of  the  surface  corresponding  to  the  varying  phase  of  an 
Internal  wave  will  have  their  backscattered  energy  misplaced  in  the 
SAR  image  as  discussed  in  Section  2. 5. 2. 5.  If  the  velocity  varies 
significantly  over  a  resolution  cell,  the  image  will  be  smeared 
along  the  azimuth  direction.  For  SEASAT  SAR  a  surface  current 
velocity  of  ux  m/s  directed  in  the  radar  wave  plane  of  incidence 
(along  x'  in  Fig.  2-1)  is  sufficient  to  shift  backscattered  energy 
by  Ax  where 


(2-16) 


Ax  -  U  (R/V)  sin  0' 
x 


and  for  SEASAT  SAR 


Ax  *  50  u  . 

x 

Hence  a  surface  current  along  x  of  50  cm/s  is  sufficient  to  shift 
backscattered  energy  by  a  resolution  cell  size  (25  m).  It  has  been 
suggested  that  wake-related  internal  waves  could  have  surface 
current  velocities  comparable  to  50  ca/s.  The  phase  velocity  of 
Bragg  scattering  waves  is  -  30  cm/s  and  is  also  seen  as  a  Doppler 
shift,  and  hence  a  position  shift  by  the  SAR.  The  two  effects  are 
additive;  but  whereas  the  Bragg  phase  velocity  would  be  homogeneous 
over  a  large  number  of  resolution  cells,  the  surface  current  would 
vary  substantially  over  several  resolution  cells  for  a  ~  few 
hundred  meter  wavelength  Internal  wave  associated  with  a  ship. 

Combination  Effects.  Above  there  are  listed  several  Internal 
wave  effects  which  could  spatially  modulate  the  radar  cross  section 
mapped  into  a  SAR  image.  Depending  on  the  nature  of  the  surface 
current  field  associated  with  the  internal  wave  wake,  the  ocean 
surface  conditions  and  observantional  geometry  the  relative  strength 
and  phase  of  these  effects  could  enhance  or  degrade  the  backscatter 


cross  section  making  a  wake  pattern  visible  or  undetectable  in  the 
SAR  image  speckle  background.  An  important  point  here  is  that 
consideration  of  the  combination  of  effects  could  produce  a  result 
significantly  different  (in  terms  of  visibility)  from  the  considera¬ 
tion  of  one  effect  alone. 

Background .  The  visibility  of  Internal  wave  effects  in  SAR 
images  clearly  depends  strongly  on  the  background  upon  which  the 
wake  is  imposed.  Some  notion  of  the  background  issue  can  be 
obtained  from  the  work  of  Bagg  (1983).  The  average  ocean  back- 
scatter  strength  increases  with  wind  speed  for  SEASAT  angles  of 
Incidence.  Bagg  noticed  in  studying  128  observations  of  natural 
Internal  waves  by  SEASAT  that  over  2/3  of  the  observations  were  in 
winds  of  10  kts  or  less  while  only  1/10  of  the  observations  were  in 
wind  speeds  of  20  kts  or  more.  Thus,  as  the  mean  background  back- 
scatter  level  was  raised  by  increased  wind  speed  the  number  of 
internal  waves  observed  decreased  dramatically.  Similarly  SAR  image 
features  associated  with  wake  related  internal  waves  should  become 
less  detectable  as  wind  speed  Increases. 


2.6  Conclusion  and  Recommendations  Associated  with  Radar  Effects 


2.6.1  Conclusions 


All  important  radar  scattering  and  SAR  system  effects  must  be 
combined  with  analytical  (or  emperical)  models  of  the  surface 
morphology  and  currents  in  a  wake  to  order  to  obtain  an  accurate 
estimate  of  what  a  SAR  would  image. 

Purturbations  of  ocean  surface  roughness  and  currents  by 
Kelvin  wakes  can  produce  different  images  when  viewed  by  a  SAR  under 
different  observational  and  environmental  conditions.  This  is  also 
true  for  internal  wave  wakes. 

Photographs  corresponding  to  SAR  Images  of  ship  wakes  may 
produce  a  distorted  idea  of  ship  wake  features  because  contrast 
stretching  is  used  to  obtain  a  "Interesting"  photograph. 

First  order  Bragg  resonant  scattering  has  almost  always  been 
assumed  to  cause  ocean  surface  backseat ter  observed  by  SAR.  It  now 
appears  that  quasi-specular  and  perhaps  higher  order  Bragg  resonant 
scatter  or  mechanisms  related  to  breaking  waves  may  also  be 
Important. 
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2.6.2  Recommendations 

SAR  image  data  on  ship  wakes  should  be  examined  in  digital 
form  using  an  image  processing  system  because  photographic  prints 
can  be  misleading  if  they  are  contrast  stretched.  Some  issues  of 
interest  regarding  ship  wake  Images  are:  line  width  of  wake 
features,  multiple  or  striated  wakes,  smearing  of  wake  features 
along  the  azimuth  direction  and  length  of  wake  features. 

It  would  be  useful  to  do  a  comprehensive  computation  of  the 
SAR  image  produced  by  even  a  crude  model  of  wake  surface  roughness 
and  currents.  This  idea  would  be  to  get  a  feel  for  how  a  SAR 
interprets  the  surface  morphology  and  currents  expects  from  Kelvin 
and  Internal  wave  wakes. 

Data  on  natural  internal  waves  gathered  both  by  SAR  and  by  in 
situ  measurements  should  be  examined  to  better  understand  SAR  images 
of  internal  waves.  In  particular  is  a  relaxation  time  approach  to 
the  modulation  of  ocean  surface  roughness  viable  -  what  is  the 
relaxation  time? 

Further  research  is  needed  to  discover  the  role  of  quasi- 
specular,  higher-order  Bragg  resonant  backscatter  and  other 
mechanises  in  SAR  imaging  of  the  ocean  surface. 
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3.0  INTERNAL  WAVE  WAKES 


3.1  Internal  Wave  Wakes 


Physically,  a  ship  moving  on  the  surface  of  the  sea  generates 
an  Internal  wave  wake  by  disturbing  the  ambient  density  profile. 

The  internal  waves  are  thus  the  oscillations  of  this  profile  around 
its  equilibrium  value.  The  disturbance  can  be  produced  by  several 
effects,  including  the  potential  flow  around  the  ship,  and  the 
turbulent  wake  it  leaves  behind,  among  others.  For  simplicity,  we 
will  concentrate  here  on  the  first  of  these;  that  should  be  suffi¬ 
cient  for  a  rough  estimate,  at  least,  of  the  surface  currents  the 
Internal  wave  wake  will  produce.  As  will  be  discussed  in  the 
following  section,  it  is  these  surface  currents  which  generate  the 
radar  return  from  the  internal  wave  wake. 


The  vertical  displacement  c(x,z)  at  depth  z  (z  is,  of 
course,  positive  upwards)  and  horizontal  position  x  satisfies  the 
Internal  wave  equation 


[~_2'(~— 2  +  ?*)  +  N2(z)V2]  ?  (x,z,t) 
3t  3z 


S(x,z,t) 


(3-1 


where  S  is  a  source  and  N2(z)  5  -  ^  is  the  buoyancy 

p  a  z 

frequency.  We  specialize  to  a  moving  source  S  -  S(x  -  Ut,z)  (we 
take  U  in  the  positive  x  direction)  and  look  in  the  rest  frame 
of  the  source.  In  this  frame  S  and  C  both  become  independent  of 
time,  and  the  solution  to  (3-1)  may  then  be  written  in  the  form 


C(x,z)  -  /  ~~~~2 
(2*)Z 


/  dz’elk  *  * 


</(k)W(z,k)W  (z'.k) 
r  n  n _ n _ 

n  k2(k2U2  -  w2(k)) 


S(k,z')  . 


(3-2) 


Here  S(k,z)  is  the  horizontal  Fourier  transform  of  the  source,  and 
and  u)^  are  the  eigenfunctions  and  eigenvalues  of  the  operator 
in  (3-1),  so  that 


[c/(k)(i-s-  -  k*)  +  k2N2(z)]  W  (z,k)  -  0 
3z*  n 


(3-3) 


subject  to  the  appropriate  boundary  conditions  that  W^  »  0  at 
z  ■  0  and  z  ■  -  *  . 


The  integral  over  k^  can  now  be  evaluated  by  contour 

integration,  locating  the  poles  at  k  *  +  w  (k)/U  in  the  lower 

x  -  n 

half  plane  to  ensure  causality.  We  note  that,  for  Internal  waves, 


l/u(3<i>n(k)/3kx)  «  1  and  k^  »  k^  ;  with  these  simplifications 
(3-2)  becomes 

*  dk  ik  y 

C(x,z)  -  0(-x)  I  2^- e  y  /  d*' 

o»  (k)  m  (k) 

l  A—  W (z,k)W  (z',k)sin  x  •  S(k,z»)  (3-4 

n  k  u  n  n  u 

where  now  k  ■  jk^|  to  a  very  good  approxiaation. 

Thus  C  is  a  sum  of  contributions,  one  from  each  Internal 

wave  mode.  We  may  write  C  ■  5"  C  ,  with  ?  identified  in 

n  n 

n 

Eq.  (3-4). 

In  the  far  field,  the  Integral  over  ky  in  Eq.  (3-4)  can  be 
evaluated  by  stationary  phase.  This  yields  two  wave  systems,  one 
traveling  in  the  +y  direction,  with  ky  -  k  >  0  and  k  determined 
from  the  stationary  phase  condition 

1  3“n<k> 

tan0  s  y/(-x)  -  ^  — -  (3-5 

u  3k 

and  another  traveling  in  the  -y  direction  with  ky  -  -k  <  0. 
Evidently,  8  in  Eq.  (3-3)  is  the  angle  off  the  wake  axis.  It  is 
sufficient  to  fix  our  attention  on  only  one  of  these,  say  that  in 


the  +y  direction,  call  it  '  .  Evaluating  the  integral  over 

k  then  yields 

y 


(+)>  x  e(-x) 

C  (x,z)  - -  - 

n 


“n(k) 


/-2itx 


2  2 
.  a  0)  k  u 
1  n 


U  3k2 


(3-6) 


u>_(k) 


Wn(z,k)  /  dz*Wn(z'  ,k)S(k,z')8in(k  +  — x  +  x/4) 


where  k  is  determined  from  Eq.  (3-5). 


Now  let  us  take  as  a  source  the  potential  flow  around  the  ship 
hull  in  an  unbounded  ocean  where  the  ship  is  made  symmetrical  about 
the  plane  z  -  0  ;  call  its  vertical  displacement  Cp0T(x,z)  .  Then 

S(x,z)  -  -  N2(z)V2cp0T(x,z)  , 


or 


S(k,z)  -  k2N2(z)?  (k,z) 


Let  us  also  specialize  to  a  two  fluid  model  of  the  thermocline.  Then 


there  is  only  one  mode,  and  we  have 


N2(z)  -  g'«(z  +  d) 


with  g'  -  gAp/p  ,  where  Ap  is  the  density  jump  at  depth 
z  ■  -  d  .  We  have,  also, 


w(k) 


-  cm 

\J 1  +  co 


k _ 

cosh  kd 


and 


W(z,k) 


Thus  we  get 


1  sinh  kz 
/g  sinh  kd 
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/g 
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CpOT(k*~d>8in(ky  +  ^  X  +  w/4) 


sinh  kz 
sinh  kd 


-  d  <  z  <  0 


Let  us  simplify  further  by  limiting  ourselves  to  the  case 
led  »  1  .  Under  these  circumstances,  Eq.  (3-7)  simplifies  to 


(+) 


(x,Z) 


9(-x)  (gf)3/2 

/-line  8U3(tan6)5/2 


CP0T*k,^8ln(ky  +  ^  *  +  */4) 


sinh  kz 
slnh  kd  ’ 


-  d  <  z  <  0 

-  »  <  z  <  -  d 


(3-8) 


where 

k  -  g’/8U2tan20 


and 


u  ■  g'/4Utan0 


Finally,  to  obtain  a  completely  explicit  answer,  let  us 
suppose  the  ship  is  spherical.  Then 


Cp0T(x,z)  -  8w  V 


r*  2 


(*-> 


T37I 


(3-9) 


3-6 


and  hence 


(3-10) 


?POT^k,_d^  2  V  6 

where  V  Is  Che  ship  volume.  As  will  be  made  clear  In  the 
following  section,  the  quantity  relevant  for  radar  backscatter  Is 
the  surface  strain,  defined  by 

e  .  i£ 

3z  z-0 

From  Eqs.  (3-8)  and  (3-10),  we  see  that  the  surface  strain  Is 

e  -  -  0(-x)sin(ky  +  ^  x  +  ir/4) 

1  3  V(g*)5/2  _-g'd/8U2tan2e  ,,  ,, 

- - «T  q/o  e  •  (3-11 

128  ir(tan0)  /z 

As  a  numerical  Illustration,  suppose  we  take  the  following 
values : 

V  -  10A  m3 

U  ■  10  m/sec 

-2  2 
g'  ■  10  m/sec 

x  -  103  m 

d  ■  102  m 

0  -  2°  -  1/30 
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Then  we  find  the  strain  amplitude  to  be 


e  «  4  x  10 

Since  our  approximation  assumed  kd  »  1  ,  while  with  the 
above  parameters  kd  is  only  9/8,  this  result  is  probably  a  slight 
overestimate.  For  more  accurate  values,  one  should  use  Eq.  (3-7) 
with  the  correct  value  of  o»(k),  not  assuming  kd  »  1  . 

3.2  Radar  Observations  of  Internal  Waves 


3.2.1  Abstract 


Radar  images  of  weak  internal  waves  are  obscured  by  the 
speckle  associated  with  any  form  of  coherent  illumination.  It  can 
only  be  removed  by  averaging.  A  precise  method  is  given  for 
identifying  and  counting  the  necessary  degrees  of  freedom. 

3.2.2  Introduction 


It  is  well  known  that  large  amplitude  internal  waves  can  be 
detected  by  their  effect  on  surface  roughness.  The  basic  hydro- 
dynamic  mechanism^ ^  is  straining  of  the  surface  by  the  internal 


waves  which  changes  the  slopes  of  small,  ambient  ripples.  The  eye 
then  sees  a  change  in  the  reflected  sky  brightness.  Recently,  there 
has  been  considerable  interest  in  the  use  of  radar  to  monitor  the 
sea  surface.  The  radar  cross  section  is  a  direct  measure  of  surface 
roughness . 

Radar  images  from  SEASAT^’^  clearly  show  some  internal 
waves.  However,  these  are  believed  to  be  rather  special  waves  that 
are  both  long  crested  and  have  large  amplitudes .  It  is  interesting 
to  ask  if  a  space-based  SAR  could  detect  weaker  (more  typical) 
internal  waves.  The  problem  is  speckle.  Speckle  is  a  graininess  or 
noise  in  the  image  which  is  always  present  when  a  rough  surface  is 
illuminated  coherently.  It  can  be  removed  only  by  incoherent 
averaging.  Low  contrast  features  (e.g.,  weak  internal  waves)  will 
be  lost  in  a  random  pattern  of  speckle  unless  enough  averaging  can 
be  done. 

Detection  of  weak  internal  waves  from  a  satellite  is  not 
a-priori  hopeless.  Satellites  can  image  large  areas  of  the  ocean 
and  considerable  averaging  is  possible.  It  is  a  quantitative  issue; 
one  has  to  compute  the  signal  to  noise.  Simple  numerical  estimates 
can  be  obtained  with  the  following  two  rules: 


(i)  The  variance  in  the  radar  cross-section  is  equal  to 

its  average.  This  is  true  independent  of  the  size  of 
the  resolution  cell. 

(ii)  An  average  with  N  degrees  of  freedom  will  reduce  the 
variance  of  the  cross  section  by  a  factor  of  1//N- . 

Thus  an  average  over  N  degrees  of  freedom  will  yield  the  (true) 

average  of  the  cross  section  to  within  a  fractional  error  of 

1//N  .  Roughly  speaking,  if  an  internal  wave  makes  a  fractional 

-2 

change  e  in  the  radar  cross-section,  then  an  average  with  e 

degrees  will  be  required  in  order  to  observe  it.  Very  large  ambient 
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waves  have  c  ~  .1  ;  a  more  typical  number  is  e  ~  10  . 

The  subtle  issues  arise  when  one  tries  to  rigorously  count  the 
available  degrees  of  freedom.  It  seems  intuitively  clear  that 
looking  at  several  pieces  of  the  ocean  will  produce  a  better 
average,  but  how  about  looking  at  the  same  patch  at  different  angles 
or  (radar)  frequencies?  The  purpose  of  this  report  is  to  give  a 
definitive  answer  to  these  questions.  The  calculations  are  a  bit 
longer  than  one  would  like,  but  the  answer  is  simple  and  satisfying. 
In  particular,  the  way  in  which  the  fluctuations  depend  on  the  size 
of  a  resolution  cell  is  worked  out  in  detail.  It  turns  out  that  if 


one  uses  all  available  Information,  then  the  size  of  a  cell  is 
irrelevant.  However,  this  is  generally  not  true  in  practice. 

It  is  simplest  to  work  things  out  for  an  ordinary  monostatic 
radar  system.  The  bulk  of  the  text  is  devoted  to  this  analysis. 
SGASAT  is  actually  a  synthetic  aperture  radar  (SAR).  This  intro¬ 
duces  some  further  complications,  but  once  the  simple  monostatic 
case  is  understood,  it  is  not  hard  to  make  the  transition  to  SARs. 

All  the  calculations  are  done  under  the  assumption  that  the 
Born  approximation  is  valid.  The  results  are  actually  more  general, 
but  the  proof  will  not  be  given  here.  In  practice,  the  Born 
approximation  is  usually  adequate.  Also,  the  problem  will  be  set  up 
in  terms  of  a  scalar  wave  equation.  The  results  apply  separately  to 
either  vertical  or  horizontal  polarization. 

The  report  is  organized  as  follows.  Section  3.2.3  is  a  sum¬ 
mary  of  the  statistical  properties  of  ripples  and  the  way  in  which 
they  are  affected  by  an  internal  wave.  The  information  contained  in 
the  ripples  is  computed  in  Section  3.2.4.  Section  3.2.5  discusses 
simple  radar  scattering  and  the  transition  to  a  SAR  is  carried  out 


in  Section  3.2.6 


3.2.3  Wave  Statistics 


Let  h(r)  be  the  height  of  the  ocean  surface  at  the  point 
r  •  (x,y)  and  let  a(k)  be  its  Fourier  transform 
h(r)  -  j  d\  a(k)  e1-*-  .  It  will  be  assumed  that,  in  the  absence 
of  a  perturbing  internal  wave,  the  a(k)  are  independent  random 
variables  with  a  Gaussian  distribution.  The  unperturbed  statistics 
are  then  specified  by  the  correlation  function 

<a(k)  a*(k')>  -  d2  (k-k*)  F(k)  (3-12) 

where  F(k)  is  adequately  described  (Ref.  1)  by 

F(k)  -  |  k"4  ,  B  -  0.005  (3-13) 

for  the  wavelengths  (tens  of  centimeters)  which  will  be  of  Interest. 

Physically,  one  can  think  of  the  ocean  surface  as  a  sum  of  plane 
ik#r 

waves  a(k)  e  ~  ~  whose  phases  are  random  and  whose  amplitudes  have 

1/2 

the  average  value  [F(k)j 

Now  supppose  that  this  random  collection  of  plane  waves  is 
allowed  to  interact  with  a  weak  internal  wave  whose  wavelength  and 
period  are  long  compared  to  those  of  the  relevant  ripples.  What 
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will  happen  is  that  a  plane  wave  of  amplitude  |a|  and  wave  number 
k  will  be  converted  into  a  new  wave  of  amplitude  |a'|  and  wave 
number  k*  .  The  phases  remain  random.  Thus,  the  effect  of  an 
internal  wave  is  to  take  a  distribution  of  random  plane  waves  and 
replace  it  by  a  different  one.  For  a  weak  internal  wave  the  changes 
6k  »  k'  -  k  and  6 | a |  =  |a'|  -  |a|  will  be  small.  Also  for 
small  | a j  ,  which  is  the  situation  in  practice,  6k  will  be 
independent  of  ja|  and  6 | a |  will  be  proportional  to  |a|  .  Thus 
6k  and  6  j  a | / 1 a |  ,  in  contrast  to  5 1 a |  itself,  are  not  random 
variables  but  are  completely  determined  by  k  and  the  internal 
wave.  This  allows  one  to  write 


2fi|a(k)|  a 

F'(k)  -  F(k)  =  6F(k)  -  “ppy  FCk)  -  6k-  F(k) 
6 |a(k) |  4k  •  6k 

-  (2  paTg)J  +  ■  -~2"  )  F(k)  -  2Y(k)  F(k) 


(3-14) 


where  F'  is  the  new  distribution  function,  Eq.  (3-13)  has  been 
used,  and  the  last  line  defines  Y(k)  .  The  factor  Y(k)  is  a 
dimensionless  number  which  can  be  computed  from  hydrodynamics.  It 
is  a  numerical  factor  of  order  one  times  the  internal  wave  strain 
and  is  small  for  weak  (linear)  internal  waves.  The  precise  form  of 
Y  will  not  be  needed  here. 
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♦ 


The  total  effect  of  an  Internal  wave  la,  then,  to  change  the 


probability  distributions  of  the  independent  random  quantities 
a(k)  .  No  coherent  effects  or  new  correlations  are  introduced.  If 

v 

the  unperturbed  a's  were  Gaussian-distributed,  so  are  the  perturbed 
ones  a  +  6a  .  The  function  6F(k)  in  Eq.  (3-14)  therefore 
contains  all  available  information.  It  is  important  to  keep  this  in 
mind. 


3.2.4.  Information  Contained  in  the  Ripples 

The  wavelength  of  interesting  internal  waves  is  always  much 

longer  than  that  of  the  surface  ripples.  One  can  take  advantage  of 

these  quite  different  scales  of  length  in  the  following  way.  Divide 

the  area  to  be  observed  into  square  patches  whose  linear  dimension 

L  is  small  compared  to  the  internal  wavelength.  Label  these 

patches  by  a  -  1,  2,  ...  ,  and  let  Ro  be  the  center  of  the 

a-th  patch.  Next  choose  some  length  l  which  is  long  compared  to 

the  ripple  wavelengths  of  interest  but  small  compared  to  L  . 
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Subdivide  each  L-by-L  patch  into  (L/i)  smaller  squares  of 
linear  dimension  l  .  Let  R  +  r.  (1*1,  2,  .  .  .  , )  be  the 

vqj 

center  of  the  small  patches  in  the  a-th  large  patch  and  define  the 


Fourier  coefficient 


A4(r  ,  k  ) 

i'-'-cr  'Tun-' 


1 


-ik  •  x 

j  e  ''nm  h(R  4.  j-  +  x)d2x  (3-15) 

<2w)2 

where  k  *  (2irn/4,  2xm/4)  and  n,  1  ■  0,  1,  2 . The 

~nm 

quantity  A^R^,  k^)  is  evidently  a  local  Fourier  transform  which 
is  capable  of  resolving  wave  numbers  which  differ  by  more  than 
2ir/4  .  By  the  definition  of  4  ,  this  is  adequate  to  resolve  the 
ripples  of  interest.  Furthermore,  by  the  definition  of  L  ,  the 
statistical  properties  of  A^(Ra,  k^)  must  be  (nearly)  independent 
of  i  .  In  fact,  no  information  is  lost  if  one  considers  only  the 
average  quantity 

g(r  ,  k  )  =  (:J-)2  y  Ia.(r  ,  k  )|2  .  (3-16) 

'"-a*  ''nm'  ^L'  “  1  i'"-a*  'nm;| 

Any  method  of  observing  the  ocean  surface  that  determines 

g(r  ,  k  1  provides  all  the  information  that  is  relevant  for 

detecting  the  Internal  waves.  Any  more  detailed  information  about 

the  sea  surface  is  superfluous.  The  phases  of  the  A^  are  purely 

random  and  therefore  uninteresting.  The  individual  amplitudes 

|ai|  are  also  uninteresting.  Since  the  internal  wave  is  supposed 

to  change  little  over  a  length  L  ,  the  individual  A^'s  differ 

from  G  only  in  that  they  fluctuate  more.  According  to  what  was 

2 

said  in  the  last  section  these  fluctuations  of  A^  around  its 


average  are  the  same  in  the  presence  or  absence  of  internal  waves 


average  value  of  G  In  the  absence  of  an  internal 


wave  is 


<G(R  „  k  )>  -  F(lc  )  ,  Ik  »  2w/£ 

Vvnm'  1  ^nm 


(3-17) 


with  the  stated  proviso  that  Ik  I  is  large  compared  to  2w/£  . 

» ^nm  * 

With  a  perturbing  internal  wave,  the  shift  in  the  average  is 


p^k  ]  "  2T(Sa’  iSom) 

'•'Tim-' 


(3-18) 


where  Y  depends  on  the  patch  location  R  . 


Having  decided  what  there  is  to  be  measured,  the  next  question 
is  signal  to  noise.  For  Ik^J  »  2w/ A  it  can  be  shown  from 
Eqs.  (3-13)  and  (3-16)  that  the  G(Ra,  k^  )  are,  to  a  good  approxi¬ 
mation,  independent  Gaussian-distributed  random  variables  satisfying 


<G(R  ,  k  )G(R  ,,  k  ,  -  6  ,6  ,6  ,  [l  +  2(r)2]F2(k  )  , 

v-o*  '*nm-'  a ’  •wn,m,y  aa1  nn'  mm' L  ML'  1  'mm' 


k  »  2tt/£ 
-nm 


(3-19) 


A  reasonable  definition  of  signal-to-noise  £  is  then 


r2  r  L2  1  !<G(5a.  )>' 

1  ’  „L  *2  2  ^ 


l  i_  y(r  ,  k  )  . 

2  ~a  ~nm 

nma  £ 


(3-20) 


The  quantity  e  L  has  a  simple  interpretation.  It  is  the  proba¬ 
bility  of  "false  alarm"  when  an  ideal  matched  filter  is  used. 


To  compute  £  it  is  convenient  to  replace  the  sum  \  by  the 
2  2* 

integral  (i/2ir)  /  d\  .  Similarly,  the  sum  \  can  be  replaced  by 

_2  a 

L  /  dR  if  6<G>  varies  slowly  with  Rq  .  The  definitions  of 
l  and  L  were  just  such  as  to  guarantee  these  conditions.  Thus, 
regarding  G  as  a  function  of  the  continuous  variables 
R  and  k  ,  one  has 

v  v 


l  -  (4 


2  (2w)2 


,  2  «<G(R,  k)>2  1/2 

/  d  xd  k  — to - ) 


F(k) 


^  [J  d2kd2Ry(R,k)]1/2  . 

(3-21) 


Notice  that  the  lengths  £  and  L  have  cancelled  out  of  the 
final  answer. 


All  of  the  above  discussion  refers  to  a  single  observation  of 
the  ocean  surface.  One  might  make  several  observations  at  different 
times.  If  the  observations  are  done  over  a  time  T  ,  the  number 
of  statistically  Independent  looks  will  be  T  divided  by  the 
relaxation  time  of  ripples  tg  .  Then  by  the  usual  rules  of 
probability,  one  can  define  an  effective  7 


3.2.5  Radar  Scattering 


For  scattering  of  microwave  radar  off  small  ripples,  the  Born 
approximation  is  adequate.  If  the  radar  illuminates  an  area  S  the 
scattered  power  is  proportional  to 


P  « 


ik  •  x  .  ,  .  .2 
e  ~  ~  h(x)  d  x 


(3-23) 


where  k  is  the  change  in  the  radar  propagation  vector.  Evidently, 

2 

if  the  illuminated  area  S  is  a  square  of  area  t  then  the 

scattered  power  is  just  | A^Ck) |  defined  in  the  previous  section. 

2  2 

The  power  averaged  over  the  L  /l  radar  resolution  cells  in  a 
patch  is  proportional  to  G(R,k)  .  To  get  the  normalization,  it  is 
easiest  to  take  the  ratios.  If  a  is  the  unperturbed  radar  cross 
section  averaged  over  a  patch  and  6a  is  the  variation  due  to  an 
Internal  wave  then 


6o(R,K) 

-w  ■  2T(S'S> 

V 


(3-24) 


Detailed  formulas  for  y  can  now  be  obtained  by  comparing  with 
Ref.  (2). 


Taking  the  average  over  a  patch  assumes,  of  course,  that  l 
is  of  a  reasonable  size.  What  is  a  reasonable  size?  Remember  that 
the  ratio  1/ L  canceled  out  of  the  final  expression  for  signal  to 
noise.  In  fact,  one  can  easily  see  that  l  could  have  been  any 
length  larger  than  a  ripple  wavelength  and  smaller  than  or  equal 
to  L  .  Thus  it  appears  that  the  signal  to  noise  is  independent  of 
the  radar  spot  size.  In  a  sense  it  is. 

Although  the  formal  expression  for  £  is  independent  of 

l  ,  one  suspects  that  somehow  a  price  must  be  paid  if  l  is 

large.  This  is  Indeed  true.  In  Eq.  (3-20)  for  £  one  has  to  sum 

over  all  values  of  k  »  (2irn/i,  2ma/i)  up  to  the  point  where 

'-nm 

2x/|k  I  is  smaller  than  the  smallest  ripple  wavelength  of 

» ^nm  • 

2 

interest  •  There  are  roughly  (*/\d_n)  terms  in  the  sum. 

Therefore  when  X.  is  large,  one  has  to  take  measurements  for  many 
values  of  k  in  order  to  even  approximate  the  ideal  signal-to-noise 
ratio.  This  would  have  to  be  done  either  by  changing  the  radar 
frequency  or  by  looking  at  the  same  spot  from  many  different  angles. 

The  above  discussion  about  radar  is  a  good  illustration  of  an 
interesting  general  principle.  The  length  i  can  be  thought  of  as 
the  resolution  length  of  the  system  in  coordinate  space.  The  ideal 
signal-to-noise  ratio  is  thus  obtained  only  if  the  resolution  in 


wave  number  space  is  the  best  possible,  namely  2ir/£  .  If  the 

resolution  in  k-space  is  changed  to  £  2k /i  with  £  >  1  ,  then 

the  number  of  available  terms  in  the  sum  over  n  and  m  in 

-2  r 

Eq.  (3-20)  is  clearly  reduced  by  a  factor  of  £  ,  and  l  is 

correspondingly  reduced  by  a  factor  of  £  .  Evidently,  the 
following  general  statement  can  be  made.  Given  a  system  with 
resolution  dr  in  coordinate  space  and  a  resolution  Ak  in  wave 
number  space,  subject  only  to  the  conditions  that 


Ar  <  a  typical  internal  wavelength 


(3-25) 


2ir/Ak  >  relevant  ripple  wavelengths 


then  the  effective  £  is 


2tt 


Jeff  Ak  Ar 


(3-26) 


where,  of  course,  the  usual  inequality  Ak  Ar  >  2ir  must  be 
satisfied.  Furthermore,  if  either  of  the  Inequalities  in  Eq.  (3-25) 
is  violated  so  that  the  characteristic  patterns  are  not  resolved, 
then  is  less  than  that  given  in  Eq.  (3-26).  Note  that  the 

conditions  of  Eq.  (3-25)  allow  2ir/Ak  Ar  to  be  quite  small.  The 
extra  "resolving  power"  needed  to  achieve  the  ideal  has  nothing  to 


do  with  the  characteristic  sizes  of  the  internal  waves  or  ripples. 
Its  role  is  simply  to  reduce  the  effect  of  fluctuations  in  the 
ambient  background. 

Evidently,  it  is  possible  to  make  tradeoffs  between  resolution 
in  k  and  in  r  .  Radar  generally  has  a  poor  Ar  and  therefore 
requires  a  high  resolution  in  k  to  achieve  good  statistics.  High 
resolution  in  k  means,  of  course,  the  ability  to  accurately  map 
out  G  as  a  function  of  k  ,  not  just  making  precise  masurements  at 
one  value  of  k  . 

3.2.6  Synthetic  Aperture  Radars  (SARS) 

A  SAR  like  that  on  SEASAT  gets  it  azimuthal  resolution  from 
doppler  information. 

It  is  convenient  to  do  the  analysis  in  three  steps.  In  step 
one  the  sea  surface  is  assumed  to  be  Independent  of  time.  The 
motion  of  ripples  and  internal  waves  is  included  in  step  two  and  in 
step  three  the  remaining  effects  are  considered. 

3.2.6. 1  Frozen  Sea  Surface .  When  there  is  no  motion  of  the 
sea  surface  all  doppler  effects  are  due  to  the  velocity  of  the  SAR 


platform.  In  this  case  the  analysis  of  sections  (3)-(5)  requires 
only  trivial  changes. 


From  the  linearity  of  Maxwell's  equations  and  the  Born  approx¬ 
imation,  it  follows  that  after  the  synthetic  aperture  has  been 
formed,  the  scattered  power  can  be  written  as 


P  * 


*  h(x)  W(x)  d2x 


(3-27) 


where  W  is  a  weighting  function  which  depends  on  the  precise 
algorithm  used  to  form  the  aperture.  If  W  were  equal  to  one 
inside  of  a  region  S  and  zero  outside,  then  Eq.  (3-27)  would  be 
identical  to  Eq.  (3-23).  It  is  not,  but  the  point  is  that  a  SAR  is 
designed  so  the  W  is  almost  equal  to  unity  in  a  resolution  cell  and 
is  nearly  zero  outside.  Thus  for  a  well  designed  SAR,  Eqs.  (3-23) 
and  (3-27)  are  almost  indentical  and  consequently  the  analysis  of 
sections  (3)  and  (5)  remains  valid. 


3. 2. 6. 2  Ripples  and  Internal  Waves.  For  backscatterlng  from 
a  ripple  convected  by  an  internal  wave  there  is  a  doppler  shift 

Im  -  +  «k  +  k  •  U  (3-28) 
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where  *  /  gk  (for  gravity  waves)  is  the  intrinsic  frequency  of 
the  ripple,  k  is  the  scattering  vector  (twice  the  horizontal  wave 
number  of  the  radar,  and  U  is  the  internal  wave  water  velocity. 

v 

This  extra  doppler  shift  causes  a  SAR  to  slightly  displace  the  image 
of  the  ripple.  It  is  sufficient  to  work  to  first  order  in  Au>  : 
the  two  terms  in  Eq.  (3-28)  then  represent  independent  additive 
effects.  The  first  term,  iii^  ,  is  a  purely  kinematic  effect  which  is 
independent  of  the  internal  wave.  It  can  be  ignored.  The  second 
term,  k  •  U  ,  produces  a  local  distortion  of  the  image  which  is 
proportional  to  the  internal  wave  strength.  In  the  language  of 
section  3  this  is  a  deterministic  effect  which  can  be  absorbed  into 
a  redefinition  of  y  .  In  radar  parlance,  it  produces  the  "velocity 
bunching"  term  in  the  SAR  cross  section  .  The  detailed  way  in  which 
velocity  bunching  enters  into  y  can  be  obtained  by  comparing 
Eq.  (3-24)  to  the  formulas  given  in  Ref.  (2). 

3. 2. 6. 3  Other  Sea  Surface  Motions.  Currents  and  other  flows 
with  a  horizontal  velocity  at  surface  act  just  like  internal  waves. 

The  only  other  phenomenon  of  any  significance  is  the  tilting 
of  ripples  by  long  surface  waves.  This  has  the  effect  of  adding  a 
non-deterministic  perturbation  to  the  oj^  term  in  Eq.  (3-28).  The 
net  result  is  a  small  random  distortion  of  the  image.  When  compared 


to  the  velocity  bunching  term  this  whole  effect  is  down  by  a  factor 

2 

of  m  ~  0.1  ,  where  m  is  the  mean  square  slope  of  the  sea 
surface.  It  is  not  likely  to  be  of  much  consequence  in  practice. 

In  summary,  to  a  good  approximation  the  only  difference 
between  an  ordinary  radar  and  a  SAR  is  that  in  the  latter  case  the 
velocity  bunching  term  should  be  Included  in  y  . 

3.2.7  Conclusions 


The  main  results  of  this  report  are: 

(i)  Eqs.  (3-20)  and  (3-21)  which  give  the  signal  to  noise 
ratio  for  any  method  of  detecting  an  internal  wave  via 
surface  roughening. 

(ii)  Eq.  (3-24)  which  relates  the  quantities  needed  in  (i) 
to  radar  cross  sections. 

(Hi)  The  fact  that  to  a  good  approximation  (i)  and  (ii) 

apply  to  SAR  provided  only  that  the  velocity  bunching 
term^  is  included  in  — 


(or  y  ) 


As  an  example,  consider  SEASAT.  An  analysis  of  the  SEASAT 

o 

geometery  and  processing  algorithm  shows  that:  (i)  the  function 

2 

W  in  Eq.  (3-23)  is  close  to  unity  in  a  region  of  area  ~  100  m  and 
nearly  zero  elsewhere  and,  (il),  SEASAT  measures  the  cross  section 
only  at  one  value  of  k  .  Then  according  to  Eq.  (3-20) 

l2  -1^2  y  (Rj  -  4  /  d2Ry(R)  (3-29 

a  l  l 

2 

where  the  k  dependence  of  y  has  been  suppressed  and  l  is  to 

2 

be  interpreted  as  100  m  .  From  Eq.  (3-29)  it  is  clear  that  the 

2 

number  of  degrees  of  freedom  available  to  SEASAT  is  one  per  100  m 
of  ocean  surface. 

Finally,  it  is  possible  to  extend  the  analysis  to  the  more 
general  problem  of  finding  a  special  Internal  wave  in  a  background 
which  contains  both  speckle  and  random  ambient  internal  waves.  At 
formal  level,  the  only  difference  is  that  there  is  now  a  random 
component  in  y  .  One  is  free  to  assume  that  the  random  part  of  y 
averages  to  zero  so  that  the  expected  values  of  G  and  a  are 
unchanged.  What  does  change  is  the  variance  of  G  (Eq.  (3-19) J. 

In  computing  this  variance  one  can  assume  that  the  fluctuations  in 
y  and  h  are  independent  (the  coupling  between  surface  and 
Internal  waves  is  weak).  The  result  is  an  additional  term  on  the 


right  hand  side  of  Eq.  (3-19)  which  is  diagonal  in  the  wave  number 
indices  n  and  m  but  is  not  diagonal  in  the  patch  indices 
a  and  a'  .  To  proceed  further  it  is  necessary  to  make  a  Fourier 
transform  with  respect  to  the  a  indices.  The  variance  is  then 
diagonal,  and  the  rest  of  the  calculation  goes  through  as  before. 
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4.0  KELVIN  WAKES 


4.1  Kinematics 


The  geometry  of  the  Kelvin  wake  behind  a  ship  is  determined  by 
the  fact  that  the  wave  pattern  is  steady  in  the  rest  frame  of  the 
ship,  and  by  the  consequent  dispersion  relation 


(4-1) 


A 


2  2 

where  k  *  /k  +  k  is  the  wave  number  of  the  wave,  and  k  and 
'  x  y  ’x 

ky  are  its  x  and  y  components  in  a  coordinate  system  in  which 
the  ship  is  moving  with  velocity  U  along  the  negative  x  axis 
and  y  is  the  crosstrack  direction. 


In  the  far  field,  the  phase  k  x  +  k  y  of  the  wave  pattern 

x  y 

will  be  stationary,  so  that 


__1 _ 

tan  8 


(4-2) 


where  8  is  the  angle  in  the  wake  measured  from  the  x  axis: 
x  -  r  cos6  and  y  ■  r  sin8  where  r  is  the  distance  to  a  point 
in  the  wake  from  the  source.  The  geometry  is  shown  in  Figure  4-1. 


Numerically,  if  U  ■  10  m/sec  ,  a  wake  angle  of  6°  corresponds  to  a 
wavelength  of  3  m;  a  wake  angle  of  3°  corresponds  to  a  wavelength  of 
less  than  1  m. 

The  angle  the  wave  crests  make  with  the  axis  is  called  $  , 
and  is  defined  by  tan*  ■  k  /k  .  Then 

*  y 

tan*  -  tan*  (4-7) 

2  +  tanZ* 

and  for  small  6  ,  *  »  28  . 

To  derive  the  geometry  of  the  Kelvin  wake  is  easy;  to  derive 
the  amplitudes  and  slopes  in  the  wake  is  hard,  especially  at  small 
wake  angles.  As  a  point  of  reference  for  more  sophisticated 
calculations,  we  describe  below  a  standard  approximation  known  as 
"thin  ship  theory." 

The  linearized  equations  of  motion  obeyed  by  the  velocity 
potential  *  ,  written  in  the  rest  frame  of  the  ship  and  with  z 
vertically  upwards,  are: 


(i)  V2*  -  0 


(ii)  3$/3z  +  U  /g  3  $/3x  »  0  at  z  -  0 


(iii)  n»$$*n«U*nU  on  hull. 

x 

The  first  of  these  says  the  flow  Is  Incompressible  and  lrrotatlonal ; 
the  second  Is  the  (linearized)  free  surface  boundary  condition;  the 
third  says  that  the  normal  (to  the  hull)  water  velocity  equals  the 

A 

(normal)  ship  speed  on  the  hull  (n  is  the  unit  normal  to  the  hull). 
The  surface  displacement  ;  can  be  calculated  from  the  velocity 
potential  through  the  (linearized)  relation 

r  u  ii 

g  3x  z-0 

Eqs .  (i)  through  (iii)  define  the  "Kelvin-Neumann"  problem, 
which  has  never  been  solved.  (A  further  discussion  of  it  appears  in 
Section  A. 3.)  "Thin  ship  theory”  consists  of  replacing  Eq.  (iii)  by 
the  much  simpler  condition 

(iii)'  n  1^*  n  U  at  y  *  +  0  ; 

y  3y  x  J  ’ 


that  is,  the  actual  hull  is  replaced  by  an  infinitely  thin  (in  the 
transverse  direction)  plate.  The  solution  to  (i),  (11),  and  (iii)' 


can  be  written  down  easily,  and  one  obtains  for  the  amplitude  at  a 
point  x,y  in  the  wake 


C(x,y) 


U2  /a 


♦<V 


(4-8) 


where 


F(8) 


1  fA+i'1/4  /A(A-l) 

2  1  2J  J  A  -  2 


(4-9) 


and 


•  o  -ik  x 
/  dx  /  dz  e  x 


("  '* K z  3Y(*,0 


dx 


(4-10) 


In  Eq.  (4-10),  kx  is  given  by  Eq.  (4-3a),  and  y  -  Y(x,z)  is  the 
equation  of  the  ship's  hull.  Thus  the  integrals  over  dx  and  dz 
actually  extend  only  over  the  ship's  side,  where  Y  *  0  . 


The  thin  ship  solution,  for  divergent  waves,  gives  wave  ampli- 
3/2 

tudes  behaving  like  8  for  small  6  ,  near  the  wake  axis,  and 

-1/2 

wave  slopes  behaving  like  8  .  In  this  approximation,  then, 

larger  wave  slopes  occur  at  smaller  angles.  Obviously  this  solution 


fails  at  small  9  ,  presumably  where  6  is  of  the  order  ship  beam 
divided  by  ship  length.  If  the  correct  theory  cuts  off  the  slope 
for  small  9  ,  we  would  expect  a  maximum  in  slope  at  some  small  wake 
angle;  possibly  this  maximum  could  explain  the  SEASAT  observations. 
We  must  therefore  find  an  improvement  on  the  thin  ship  solution, 
valid  as  9  -*■  0  ,  to  see  if  such  a  maximum  actually  occurs. 

Attempts  to  do  this  occupy  Sections  4.2  and  4.3. 

4.2  Kelvin-Neumann  Problem  (A) 

4.2.1  Introduction 

Wakes  produced  by  ships  have  been  subject  to  investigation  for 
many  years. There  has  been  particular  Interest  recently  because 
of  the  pictures  obtained  with  SEASAT. 

Here  we  present  an  approximation  technique  to  calculate  wakes. 
As  a  basic  assumption  we  assume  linearization  is  acceptable. 

Further,  we  will  assume  low  Froude  number.  For  some  cases  this  may 
be  of  Interest.  In  any  event,  the  structure  of  the  solution  will  be 
instructive.  In  particular,  it  will  be  found  that  in  this  limit  the 
displacements  are  indeed  small.  This  lends  credibility  to  the 
linearization  procedure. 
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It  has  been  stated that  the  solutions  have  essential 
singularities  at  zero  Froude  number.  In  our  treatment  it  will  be 
seen  how  these  can  come  about.  Such  a  singularity  does  not  mean  we 
cannot  approximate  for  small  Froude  number.  It  merely  means  that  a 
convergent  power  series  cannot  be  obtained. 

As  an  example  consider  the  function 

xsin  — 

f(x)  -  e  X  -  1 

This  indeed  has  an  essential  singularity  at  x  *  0  .  However, 
for  small  (real)  x  it  is  well  approximated  as 

f(x)  -  x  sin( 1/x) 

Our  method  consists  of  using  Green's  theorem  to  form  an 
integral  representation  for  the  desired  potential  function.  Passing 
to  the  limit  where  the  field  point  is  on  the  body  of  the  ship  yields 
an  Integral  equation  for  the  potential  there.  For  small  Froude 
numbers  the  kernel  is  essentially  that  for  the  potential  of  the  flow 
past  the  ship  plus  its  reflection  in  the  ocean  surface  but  without  a 
free  surface.  The  approximation  consists  in  using  the  solution  of 
the  simplified  problem  as  the  zeroth  order  solution  of  the  integral 
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equation.  The  solution  elsewhere  is  then  obtained  from  the  integral 
representation.  Higher  order  corrections  (including  neglected  non¬ 
linear  terms)  can  then  be  found  by  iteration — at  the  cost  of 
increasing  complexity. 

This  note  is  organized  as  follows:  In  Section  4.2.2  the 
general  form  of  the  linear  problem  is  given.  This  immediately 
suggests  the  form  for  a  perturbation  theory.  The  possible  occur¬ 
rence  of  an  essential  singularity  is  made  clear  in  Section  4.2.3. 
Further,  an  ambiguity  which  is  to  be  resolved  by  physical  considera¬ 
tions  is  found. 

To  see  with  clarity  how  the  approach  is  Implemented, 

Section  4.2.4  is  devoted  to  the  two  dimensional  case.  Finally,  in 
Section  4.2.5  the  result  is  obtained  for  a  two  dimensional  example. 
This  is  immediately  generalized  to  a  large  class  of  bodies. 

4.2.2  Formulation 


Consider  a  ship  moving  with  velocity  U  in  the  direction  x 
on  a  sea  with  surface  z  -  0  .  Transforming  to  a  coordinate  system 
at  rest  with  respect  to  the  ship  we  can  write  the  velocity  potential 
as 
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Notice  that  Eq.  (4-18)  plays  two  roles: 


(1)  If  r  Is  on  Sj  this  gives  an  inhomogeneous  integral 


equation  for  $  on  S.  ,  (inhomogeneous  since  — — 

1  3n* 


is  given). 


(2)  Given  $  on  we  then  have  $(r)  everywhere. 


What  is  needed  is  an  approximation  to  the  solution  of  the 


integral  equation.  For  this  we  first  construct  the  Green's 


function. 


4.2.3  The  Green's  Function 


This  is  readily  constructed  by  Fourier  transforming  in  x 


and  y.  The  result  can  be  written  as 


G(r',r)  -  G  (r ' ,r)  +  G.(r',r) 


(4-19) 


where 


g0<S'»E>  “  “  —  (7-^-7  +  — ~ 

4tt  r'-r  rf-r 

I  V  V  I  V  V  I 


(4-20) 


„*  >  .*•  7.  .7 


.A. 


i 


r  -  (x,y ,-z) 


and 


•  V 


(2ir)‘ 


If 


d2k 


k  [*|-k2. 
U 


k(z+z ' )  i[kx(x’-x)+k  (y'-y)] 
e  e  y 


(4-21) 


k-  + 

A  number  of  points  should  be  noted: 

(1)  Gq  Is  precisely  the  Green’s  function  which  would  be 
used  to  calculate  the  potential  flow  in  an  Infinite 
medium  caused  by  a  body  consisting  of  the  ship  and  its 
image  in  the  z  ■  0  plane  moving  with  velocity  U  . 

In  particular,  then,  the  solution  of  the  integral 
equation  with  G  replaced  by  GQ  can  be  determined 
by  any  of  the  classical  methods. 


For  small  Froude  number  Is  In  some  sense  small. 

Thus  if  we  Introduce  dimensionless  quantities  in  terms 

of  a  typical  length  of  the  boat,  the  Eq.  (4-21) 

2 

involves  g  /U  in  the  denominator.  For  this  number 
large  should  be  small.  The  obvious  perturbation 

theory  is  then  to  take  as  the  lowest  order  solution  of 
the  integral  equation  the  solution  q  of  the 
equation  with  G  replaced  by  GQ  .  In  particular, 
the  lowest  order  wave  field  far  from  the  ship  can  be 
computed  from  Eq.  (4-18)  with  replaced  by  . 
Iterating,  we  can  compute  corrections  to  on  the 
body  and  then  corrections  to  the  wave  field. 

In  this  way  we  do  not  obtain  a  power  series  in 
2 

U  /g£  .  Such  a  power  series  would  be  obtained  if  in 
Eq.  (4-21)  we  introduced  the  expansion 


However,  the  terms  in  this  series  are  successively 

more  and  more  divergent.  This  is  the  origin  of  the 

"non-analytic”  behavior  at  zero  Froude  number. 

Successive  terms  in  our  iterative  series  will  indeed 

be  smaller  and  smaller — but  it  will  not  be  a  power 
2 

series  in  U  /£g! 

(4)  The  G  of  Eq.  (4-21)  looks  peculiarly  symmetric  in 

x  and  -x  .  But  we  expect  waves  to  appear  far  behind 
the  ship  and  not  far  in  front. 

(5)  There  is  an  ambiguity  in  Eq.  (4-21).  Indeed  singular¬ 
ities  occur  where  the  denominator  vanishes.  How  are 
these  to  be  treated?  Physically  the  reason  for  the 
singularities  is  clear.  They  occur  at  wave  numbers 
such  that  the  corresponding  free  surface  waves  travel 
with  velocity  U  . 

The  points  (4)  and  (5)  are  intimately  connected.  Indeed,  as 
we  will  see,  the  treatment  of  the  singularities  is  completely 
determined  by  the  requirement  that  there  be  waves  far  behind  the 
ship  but  not  in  far  in  front. 


It  may  be  noted  that  with  our  choice  of  Green's  function  the 
integral  around  the  water  line  in  Eq.  (4-18)  is  not  small.  Indeed 


it  is  of  order 


However,  we  will  see  that  part  of  the  integral  over  cancels 

this  identically  and  so  we  will  be  left  with  an  expression  of  the 
form 


♦  -  ««%! 


These  points  become  especially  clear  if  we  look  at  the  two 
dimensional  problem  to  which  we  now  turn. 


4.2.4  The  Two-Dimensional  Problem 


The  formulas  become  much  more  transparent  when  we  consider  the 
"ship"  to  be  infinite  in  the  y  direction  and  therefore  omit  the 
dependence  in  this  direction.  With  obvious  modifications  we  can 
take  over  the  previous  work.  Thus  the  integral  representation  (and 
the  Integral  equation)  become 


/  {G(r\r)  -  *<r’)  -  t(r')  —  G(r',r)}  dS 

SL  3n '  3n'  1 

+  —{♦(.')  —  G(r',r)  -ii-o(r'.r))  , 

g  ax'  3x' 

2 

~  ~  {♦(£')  - —  G(r’,r)  -  G(r *  ,r) }  , 


Thus  the  Integral  over  the  surface  of  the  boat  Is  replaced  by  one 
over  the  arc  Sj  ,  and  the  integral  over  the  waterline  is  replaced 
by  contributions  from  the  front  and  back. 


G  becomes : 


G  -  G  +  G, 
o  1 


0„<r',E)  -  1/2.  {in  |r'-r|  +  In  |r'-r|) 


»  dk  k 

Gl(-'  "  2*  $  TFT  Tfk 


[- 


|k  |(z+z')  ikx(x'-x) 
e  e 
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-k2] 
x J 


which  Is  conveniently  written  without  absolute  values  as 


,  o  ..  k  (z+z')  lk  (x'-x) 

_  J_  r  <fcy  e  *  * 

2x  J-  k  -  a- 


x  u2 


dk 


2v  J  .  2 

°  kx  "  g/U 


k  (z-z')  ik  (x'-x) 
x  x  x 

e  e 


(4-25 


Now  It  Is  Intuitively  clear  that  the  dominant  contribution 


to  Gj  for  | x ' —x |  large  will  come  from  the  poles,  at 


k  *  +  g/U  .  For  simplicity  consider  z  ■  z'  *  0  .  Then  Gj  ca 


be  written  as 


,  •  dk  lk  (x'-x) 

1  j  x  _  x' 


2x  -•  k  + 

x  ^2 


.  *  dk  ik  (x’-x) 

+  —  /  x  -  x 
9ir  J 


2ir  -  k  -  -4 
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X  y2 


.  *  dk  ik  (x'-x)  1  *  dk  lk  (x'-x) 

.  1  (  x  x  r  x  x 

+  “  /  - 7  e  - / 


2*  o  k  +  -f 

x  y* 


2x  -  [kx  -  -f] 

X 


Now  integrating  by  parts  we  readily  see  that  for  large 
x-x'|  the  sum  of  the  last  two  terms  is  o(— i— ] 


1 


4.2.5  An  Example 


To  Illustrate,  we  consider  the  case  where  the  "ship"  Is  an 
Infinite  semi -elliptic  cylinder,  i.e..  In  the  x,z  plane  the  bottom 
of  the  ship  is  described  by  the  lower  half  of  the  ellipse 

2  2 

^  +  ^2  -  1  .  (4-30) 

a  b 

The  potential  flow  around  the  full  ellipse  described  by 
Eq.  (4-30)  is  readily  found. It  can  be  written  so: 

If 

x  *  2.2  cosh  u  cos  6 

a  -b 


z  *■  2.2  sinh  u  sin  0 

a  -b 


then 


*  -  Ub  ~  cos  0  e^  .  (4-31) 

O  O 


With  K  ■  cosh  u  ,  the  surface  S,  is  given  by  £  *  £  where 

a  o 
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For  a  point  on  Eq.  (4-31)  gives 


$  -  bU  cos  6 

o 


(4-32a) 


34 

To 

3n 


_ bU  cos  6 _ 

(a2-b2)(Co2-cos2e) 


(4-32b) 


To  calculate  the  surface  waves  produced  we  use  Eq.  (4-22) 
with  G  replaced  by  the  asymptotic  form  of  Eq.  (4-28).  (The  terms 
in  Eq.  (4-22)  involving  GQ  merely  reproduce  $  or  give  other 
terms  which  fall  off  algebraically  behind  the  ship.)  Let  F  denote 
the  contribution  to  $  due  to  the  integral  over  ,  i.e., 


fi  +  f2 


where 


F  -  /  G5a8)(r\r)  —  *  (r')  dS 
1  «  1  3n '  °  ~  1 


and 


*  -  -J  ♦  <*’>  — — -  G^SS^(r ' ,r)  dS, 


(4-33) 


(4-34) 


We  obtain: 


£-{z-ix) 

.2  2x 


iUb  /  cos 9  d0 


[bsin  9  +  iacos  9] 
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+  complex  conjugate 


(4-35) 


and 

^2<z-ix)  2if  ^2[b8in9+iacos0] 

F2  *  1  ^  /  cos9  d0  [asin9  +  ibcosQ]  e^ 

+  complex  conjugate  (4-36) 

These  integrals  are  readily  evaluated  in  the  limit 

2 

gb/U  »  1  (which  we  have  been  assuming  all  along).  Thus  only  the 
immediate  vicinity  of  0  ■  x  and  9  *  2x  will  give  contributions. 

Thus  the  first  approximation  to  the  integral  in  Eq.  (4-35)  is 

_  [bsin0+laco80] 

zx  z 

/  cos9  e  d0 


(4-37) 


-iga  iga  -g^ 

.[-e^  +eu2]/’eu2  de-^isinfif 

o  8  U 


£-<z-ix) 


1-  •  &  •  y 


F1  - 


3  2 

-2U  U 

- e 

g 


sin  +  c.c. 
IT 


Similarly,  in  lowest  order  we  find 


F2 - 2Ub  e 


^(z-ix) 

u2 


cos  +  complex  conjugate 
U 


(4-38) 


One  readily  checks  that  the  contributions  of  Eqs.  (4-37)  and 
(4-38)  exactly  cancel  the  "end  point"  terms  in  Eq.  (4-22).  Remark: 
This  is  a  general  result.  The  lowest  order  terms  of  the  surface 
Integral  cancel  the  line  integral  in  our  representation. 


To  get  a  non-zero  radiated  field  we  must  calculate  to  one 
2 

higher  power  in  U  /gb  .  This  only  arises  from  F2  and  is  due  to 
the  8in0  in  Eq.  (4-36).  Our  result  for  the  radiated  field  far 
behind  the  ship  then  becomes 


£Z 

♦  *  4aU  sin  sin  ^  eU 

s  U  U 


(4-39) 


u 

Remark:  This  is  indeed  small  if  r—  is  small. 

bg 
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A. 2. 6  Generalization 


Looking  at  what  has  been  done,  it  is  seen  that  all  that  has 
been  used  were  the  values  of  at  the  endpoints  and  the  form  of 

the  surface  in  the  immediate  neighborhood.  This  suggests  a  general 
class  of  problems  for  which  the  solution  is  immediate. 

Thus  let  the  backend  be  at  -a^  and  the  front  at  &2  * 

Assume  that  the  surface  in  the  vicinity  of  these  points  can  be 
described  as  portions  of  ellipses  with  axis  ai>b}  ,  and 
respectively.  Following  our  argument  exactly  we  obtain 


A. 2. 7  Conclusion 


(1)  The  method  seems  to  produce  reasonable  results  for  small 
Froude  number. 


(2)  What  needs  to  be  done? 


(a)  Realistic  examples  in  3  dimensions. 

(b)  Estimates  of  the  corrections  obtained  by  iteration. 

(c)  Estimates  of  the  effect  of  neglected  non-linear 
terms . 

4.3  Kelvin-Neumann  Problem  (B) 

4.3.1  Introduction 

In  Section  4.2(A)  we  considered  the  wake  problem  in  the  low 
Froude  number  limit.  This  was  an  attempt  to  see  if  there  is  some 
limit  in  which  linear  theory  is  adequate. 

Using  a  coordinate  system  at  rest  with  respect  to  the  boat,  we 
obtained  an  integral  representation  for  the  velocity  potential 
<Kr)  •  When  r  is  on  the  boat  hull  (S)  this  gives  an  integral 
equation  for  4>  on  S  .  The  inhomogeneous  part  of  this  equation  is 
given  in  terms  of  the  known  normal  component  of  velocity.  The  exact 
solution  of  the  linear  problem  would  involve  solving  the  integral 
equation  and  then  using  the  integral  representation  to  obtain 
$  when  r  is  far  from  S  . 


An  approximate  solution  will  be  good  if  we  can  construct  an 
adequate  description  of  $  on  S  .  For  low  Froude  numbers  this 
seems  possible.  The  kernel  of  the  Integral  equation  decomposes  into 
two  parts.  The  first  corresponds  to  the  flow  past  the  boat  plus  its 
image  in  the  plane  z  *  0  in  an  infinite  medium.  The  second  part 
should  be  small  for  low  Froude  numbers. 

The  approximation  suggested  is  to  use  the  solution  of  the 
simple  potential  flow  problem  in  the  full  integral  representation. 

In  Section  4.2(A)  the  general  formulation  was  given.  For  illustra¬ 
tive  purposes  the  two-dimensional  case  was  treated.  It  was 
indicated  that  three-dimensional  examples  should  be  calculated. 

Here  we  consider  a  simple  but  not  completely  ridiculous  case. 

4.3.2  Boat  Model 

We  model  the  hull  as  half  of  a  prolate  ellipsoid.  To  be 
specific,  we  introduce  spheroidal  coordinates  such  that 

x  ■  c  coshn  cos  9 
y  *  c  sinhn  sin0  cosu> 
z  ■  c  sinhn  sin9  sinu 


(4-41) 

(4-42) 

(4-43) 


Then  the  hull  is  the  surface 


n  *  n  ,  0  <  9  <  it  ,  n<u<2ir 


I 

Thus  the  major  and  minor  axes  are  a  =  c  coshn  and  b  *  c  sinhn  , 

o  o 

respectively. 


4.3.3  The  Potential  Flow  Problem 


The  problem  to  be  solved  is  that  of  finding  the  solution  of 
Laplace's  equation  with  the  boundary  conditions  0  as  |r|  * 

and 


-r—  ♦  *  -n  U  (4-44 

3n  o  x 


on  the  surface 


n  *  n  ,  0  <  9  <  it  ,  0  <  w  <  2tt 

o 


In  general  orthogonal  coordinates,  Laplace's  equation  becomes 


In  our  present  coordinates 


Note: 


2  2  2  2 
h.  ■  h-  ■  c  sinh  n  cos  9  +  cosh  n  sin 
1.  £* 


m  c  sinhn  sinQ 


It  is  convenient  to  introduce 


t,  -  coshn  ,  y  *  cos0 


6 


Then 


2  2 

h!  -  h2  -  c  5  -  y 


h3  -  C  (c2-i)(i-P2) 


and  Laplace's  equation  becomes 


(i-p2)(c2-i 


a? 


3C 


3y 


3y 


The  close  relation  of  this  problem  to  the  two  dimensional 
problem  of  an  ellipse  described  in  Section  4.2(A)  suggests  looking 
for  a  solution  of  the  form 

♦0*  Ug(?)  .  (4-47) 

(One  readily  shows  that  the  boundary  condition  of  Eq.  (4-44)  can 
then  be  satisfied).  Inserting  this  ansatz  in  Eq.  (4-46)  yields 

—  (U2-  l)  -  2g(C)  .  (4-48) 

We  need  the  solution  of  this  which  vanishes  at  •  .  [An  overall 
constant  is  then  determined  by  Eq.  (4-44)1. 

One  solution  of  Eq.  (4-48)  is  found  by  inspection.  It  is 

gl(?)  *  ?  * 

This,  of  course,  does  not  satisfy  the  condition  at  infinity. 

However,  it  is  well  known  that  given  one  solution  of  a  linear  second 
order  differential  equation  a  second  can  be  obtained  in  terms  of 
quadrature.  In  this  case,  the  requisite  integrals  can  be  explicitly 


evaluated.  We  find 


One  readily  checks  that  Eq.  (4-48)  is  satisfied  and  82^)  +  0  as 
C  ♦  0  . 


Using  Eq.  (4-44)  we  then  obtain 


-  pUc 


{2  +  ?ln 

(C+l) 


2C 


+  In 


(V‘) 


=  -  pUcf(C) 


(4-50) 


'„2- 1 


<«.♦»> 


4.3.4  Asymptotics  of  the  Green's  Function 

The  radiated  wave  field  is  to  be  calculated  from  Eq.  (4-15)  of 
4.2(A)  by  inserting  4>q  for  $  ,  and  Gj  for  G  .  In  particular, 
for  the  far  field  we  need  Gj  for  large  distances  back  of  the  boat. 
Let  us  evaluate  the  asymptotic  form  in  the  low  Froude  number  limit. 


In  4.2(A)  Eq.  (4-21)  we  had  found  that 


^ _  Hz-*')  1(lt,<x'^)  +  Vy'-y)] 


/«.  Y 


•  A 


1  ff  A 


Guided  by  the  two-dimensional  case  we  expect  the  major  contri' 
bution  to  the  asymptotics  to  come  from  the  zeros  of  the  denominator 
These  are  where 


k 


2 


x 


.,2 


k 


and  kx  on  the  real  axis.  (We  remember  that  the  causality  require¬ 
ment  means  we  are  to  go  below  both  poles).  These  points  are  at 

kx  -  +  K(ky)  (4-52! 


where 


K(fcy) 


(4-53 


The  contribution  from  these  points  is  then 


.2  k(z+z ' ) 
k  e 


i[K(x'~x)+k  (y’-y)] 


We  further  approximate  this  when  |x'-x|  is  large,  by  using 
stationary  phase  to  do  the  ky  Integration.  For  simplicity  we  look 
at  the  not  uninteresting  region  where  (y— y ' |  «  (x-x'|  .  The 
stationary  phase  point 

Hr ' 0  (4'55) 

y 


yields 


k  -  0 

y 

The  stationary  phase  integration  is  then  readily  performed 
with  the  result 


-  A(x,y ,z)  B(z ' ,x' )  +  c.c. 


(4-56) 


where 


/K 


8  igy2/2U2x  g(z-ix)/u2 

V  C 

(-x)U2 


(4-57) 


and 


n  .  g(*,+ix,)/U2 


(  A-Sfl'i 


Here  we  have 

also  assumed  j x | 

»  kl  • 

|y|  »  |y* | 

4.3.5 

The  Radiated  Wave 

Field 

The  velocity  potential  far  behind  the  boat  is  now  to  be 

obtained  using  Eq.  (4-18)  of  4.2(A)  with  $  replaced  by  the  $  of 

o 

Eq.  (4-20)  and  G  replaced  by  Gj  of  Eq.  (4-26),  i.e., 

<Kr)  -  A[lx  +  I2  +  I3  +  I4  +  I5  +  I6]  +  c.c.  (4-59) 


Here, 


„2  3* 

h  mr~f  dy' 

L  g  c+  3x' 

2  3d> 

1 2--rJ  — Bdyt 

L  g  c-  3x' 


I  -  /  B  -  *  dS' 

J  S  3n»  ° 

I,  -  -  —  /  $  —  B  dy' 

4  g  i+  °  3x’ 

^“r/  *0— Bdy' 

3  8  c-  °  3x» 


I ,  -  -  /  ♦  - B  dS' 

6  S  °  3n' 


w-.’.w  ■. 


>,%  ,»  /.  .. 


sw 


Ue  now  calculate  these  integrals  to  the  first  non-vanishing  order  in 
Froude  number. 


Let  us  illustrate:  For 
line.  This  can  be  written  as 


at 

I.  we  need  — —  on  the  water 
3x’ 


(4-60) 


It  will  be  seen  that  the  dominant  contribution  to  our  Integrals  come 

2 

from  the  regions  p  *  1  .  Then  to  calculate  to  zeroth  and  first 


order  we  can  replace 


3x' 


3<f> 

o 

3x' 


(4-61) 


Then 


3 

-U  3f 
g  3? 


r  igacos6/u  2fl  . 

J  e  °  cos  0  dy 


(4-62) 


?o  c+ 


Here  c+  is  the  part  of  the  waterline  of  the  hull  for  x  >  0  .  The 
spheroidal  coordinates  of  this  curve  are  so: 


c+  is  0  <  0  < 


Then 

y  *  b  sin0  cosio  , 

z  ■  b  sind  sinui 

The  waterline  is  clearly 

a)  *  2ir  and  u)  *  n 

As  y  goes  from  0  to  b  ,  0  goes  from  0  to  tt/2  ,  and 

dy  -  b  cos0  d0  .  As  y  goes  from  -b  to  0  , 

0  goes  from  ^  to  0  ,  and  dy  -  -b  cos0  d0  .  Therefore,  the 

total  integral  over  c+  is: 

V-  n  J. 

g  0CO  °  O  ir/2 

i.e. , 

T  .  -2U3  3f  ~2~~  igacos0/U2  3a  .a 

I,  -  — - - c  C  -  lie®  cos  0  d0 

A  8  a ,  o  ' 

3Co  o 


(4-63) 


Similar  considerations  lead  to  the  result 


!  12!  c  {2_  !  /'  elSa=o,9/U2  ^3,  de 

2  g  350  °  «/2 

Putting  in  the  appropriate  surface  element  area  and  normal 
derivative,  we  have  for  the  formula 

I3  -  Uc2(Cp-  l)  —  /  d0  sin6  cos0  eigacO80/U  |  J 
3C  o 

°  (4-6 

with  (  }  -  f  . 

w 

The  principal  contributions  here  come  from  the  end  points.  To  the 
first  two  orders  we  get 


Then 


{  I  -  -2lj2 

*  '  gbsinO 


(4-6l 


21T 


?2-  1  — 

0  3{  o 


/"  eigaco80/lj2cos0  d0 


(4-6; 


Adding  Eqs.  (4-62),  (4-64),  and  (4-67)  we  obtain 


The  integral  in  Gq.  (4-72)  is  recognizable  as  a  Bessel 
function.  We  can  look  up  the  asymptotic  form  or  directly  use  the 


stationary  phase  approximation.  Thus 

/'  e1*’00867"2  cosfl  d8  -  i  iisi  -  !)  .  (4.73) 

o  *8  „2  4 

Then  Eq.  (4-72)  becomes: 

2  3/2 

I  -  -2i  /H*  UacfUJ  ±  (~)  sin(£|  -  J  )  .  (4-74) 


Using  this  result  and  Eq.  (4-57),  we  see  that  Eq.  (4-58)  says  that 


♦  -  -4Ucf (4  )  ^  !L  e^/0  8ln(Si  .  l)cog  [_*(*„  .  1] 

'  '  g  2  4  y2  2x  4 


(4-75) 


4.3.6  Conclusions 


(1)  For  smooth  hulls  of  the  type  considered  here,  we  do  seem 
to  have  a  consistent  expansion  in  Froude  number. 


(2)  All  the  contributions  to  the  far  wake  arise  from  the 
immediate  vicinity  of  the  bow  and  stern.  Presumably,  the  shape  of 
the  boat  there  will  largely  determine  a  specific  wake. 


(3)  Protuberances  which  emphasize  other  places  on  the  hull 


may  have  large  effects. 


4.4  Propeller  Wakes 


As  described  in  Section  4.1  of  this  report,  the  kinematics  of 
the  Kelvin  wake  of  the  ship's  hull  are  determined  by  the  dispersion 
relation 


/ik  =■  kxU 


(4-76) 


this  leads  directly  to  the  caustics  at  angles  of  +  19.5°  as  well  as 
to  the  wavelength  dependence  on  wake  angle  outlined  earlier.  When 
the  source  strength  is  not  steady,  but  oscillatory,  as  is  the  case 
for  the  ship's  propeller,  then  the  dispersion  relation  changes  to 


/gk  ■  k  U  +  <d 

X  o 


(4-77) 


where  is  the  angular  frequency  of  the  oscillating  source.  The 
kinematics  associated  with  the  Kelvin  wake  of  the  propeller  is  thus 
very  different  from  that  of  the  hull. 


As  before,  the  stationary  phase  point  is  determined  from 


tan6  3k 


(4-78) 


where  8  is  the  angle  in  the  wake  relative  to  the  wake  axis.  The 
caustic  occurs  when,  in  addition, 


0  . 


(4-79) 


From  Eqs.  (4-77),  (4-78),  and  (4-79),  it  is  straightforward  to 
show  that  for  an  oscillating  source,  the  caustic  angle  0c 
satisfies 


(y  +  l)(y  -  1/8)  -  —  (— -)y2  / y2  -  0  (4-80) 

/6  8  y  z  g 


where  y  »  tan  .  When  Ua)Q/g  •  0  ,  this  yields  y  -  1/8  ,  which 

corresponds  to  the  steady  source  caustics  at  8  ■  +  19.5°  . 

c  - 

However,  when  Uu>  /g  »  1  ,  the  solution  is 
o 


Tof  (“Vs  J  ' 


In  between,  y  becomes  infinite  at  Uu^/g 
y  versus  Uw  /g  is  shown  in  Figure  4-2. 


(4-81) 


-  flirt  . 


A  sketch  of 


Numerically,  typical  values  are  U  -  10  m/sec  and  ui  »  1  . 

o 

The  propeller  source  is  thus  in  the  Uu>  /g  »  1  regime.  For 

o 

to  ■  10  rad/sec  ,  we  have 
o 


The  propeller  caustics  are  thus  expected  to  be  at  quite  small 

angles.  Furthermore,  since  u>q  ,  the  propeller  blade  rate,  is 

roughly  proportional  to  ship  speed  U  ,  one  expects  that  6  is 

c 

roughly  proportional  to  U~*:  slower  ships  have  wider  propeller 
caustics . 

The  wavelengths  and  slopes  in  the  propeller  Kelvin  wake  are 

determined  by  Eqs.  (4-77)  and  (4-78).  For  small  values  of  the  wake 

angle  8(0  <  0  )  ,  the  relations  between  k  and  k  and  0  are 
c  x  y 

nearly  the  same  as  in  the  <i»q  *  0  case.  Thus  we  again  expect  large 
slopes  at  small  angles. 

The  difficulty  with  the  propeller  caustics  as  an  explanation 
of  the  SEASAT  images  lies  in  the  relatively  small  amplitudes 
probably  associated  with  them — though,  to  be  sure,  it  is  difficult 
to  make  reliable  quantitative  calculations.  As  a  model,  we  may 
describe  the  propeller  as  an  oscillatin'  pressure  distribution 


PQ(x)e  .  The  solution  of  the  linearized  problem  with  this 

source  results  in  wave  amplitudes ,  far  from  the  source  and  on  the 
caustics,  of 


|H  .  i  rq/3)  g1/6 _ (k  k  ) 

4npU  25/12  3l/3  yll/6  05/3  xl/3  W  yJ 

c 


(4-82) 


where  k  and  k„  are  to  be  evaluated  at  0*0 
x  y  c 

Eq.  (4-82),  we  have  specialized  to  the  case  0  «  1  ,  or 


In  deriving 

-1/3 


Uo>  /g  »  1  .  Note  the  characteristic  caustic  fall  off  of  x 
o 


with  distance  behind  the  ship,  in  contrast  to  the  x 
off  the  caustic. 


-1/2 


behavior 


The  source  strength  p  (k  ,k  )  is  proportional  to  exp  (-kd) , 

o  x  y 

where  d  is  the  effective  depth  of  the  propeller.  For  the  small 
angles,  and  short  wavelengths,  in  the  divergent  wave  system  of 
interest  to  us ,  k  is  very  large.  Therefore,  if  the  effective 
depth  d  is  anything  other  than  zero,  the  source  is  strongly  damped 
and  the  propeller  wake  amplitudes  and  slopes  are  far  too  small  to  be 
relevant  to  SEASAT.  Only  if  the  source  is  really  the  entire 
turbulent  flow  field  around  the  propeller,  and  if  this  flow  field 
extends  up  to  the  surface,  could  the  amplitudes  be  significant. 


APPENDIX  4A 


It  Is  Illuminating  to  treat  the  two  dimensional  case  in  some 
generality.  The  nature  of  the  cancellation  of  the  "endpoint"  terms 
and  the  implicit  assumptions  of  our  semi -ellipse  model  become  clear 

Let  the  hull  of  the  "boat”  be  described  by  the  curve 

x  »  ?<z)  .  (4A-i; 


The  endpoints  are  to  be  at  x  *  -a^  and  x  *  a2  » 


-aj  *  £(o) 
a2  -  C(o) 


-a. 


C(o) 


(4A-2; 


We  will  assume 


3z 


z-o 


0 


This  seems  to  be  a  necessary  assumption  if  we  are  to  argue  in  any 
generality.  Otherwise  the  boundary  value  problem  to  be  solved  for 
t  will  involve  a  surface  with  cusps  at  the  ends.  In  general  6 


vv* 


(or  some  of  its  derivatives)  will  be  singular  there.  Explicit 
asymptotic  results  will  depend  on  the  nature  of  the  singularity  and 
hence  be  very  model  dependent. 


The  endpoint  terms  in  Eq.  (4-22)  are 

2 

7"  W“al’°)  ~  G1(-ai»°;£) - -  (-a1.°)G1(-a1.0;r)}  (4A-3) 

g  °  1  3x'  1  3x' 

plus  similar  terms  evaluated  at  &■£  • 


The  integral  over  s^  in  Eq.  (4-22)  can  be  converted  to  an 
integral  over  z  of  the  form 


3? 


3*  (x\z') 

I-  /  to  e - =--=-4.U\z*) 

s,  3x'  3z '  8l' 


3z'  3z*  1 


(4A-4) 


For  Gj  we  use  the  asymptotic  form  given  by  Eq.  (4-29).  With 
the  assumptions  that  have  been  made  the  only  significant  contribu¬ 
tions  to  the  integral  come  from  the  vicinity  of  the  points 
x’- 


a^,z'«  o  and  x’-  a^iZ'  ■  0 


k? 


ft 


i* 


K 


tv 


k- 


a  r-.-rv  w  -  w  ;  v.  r  717^?  w .’  y — ^  -v-  vv 


'T^r^rrn^  ^  *.  vv.v.%  t  w  v 


Since 


3^ 


3z 1 


3? 


0  and  by  assumption  — — 


z  '*0 


3z' 


0  we 


z '  =*o 


readily  see  that  to  lowest  order  the  contribution  of  the  left 
endpoint  is 


I  m  $  (_a  0)  /  G.fr’.r)  dz' 

O  .  ,  O^  1  ‘  J  l'~ 

3x’ 


-  4>  (-a  ,0)  /  —  CjCr’.r)  dz» 
3x* 


(4A-5) 


Then  since  only  the  vicinity  of  z'  *  0  contribute  to  the 
integrals,  we  can  write  the  limits  as  zero  to  infinity.  But  using 
our  asymptotic  representation  we  then  obtain 


00  2 

/  Gi(j,»E)dzl  *  -7  Gi(^V0;s) 


(4A-6) 


and 


/  —  61(£,»e)<**'  “  —  —  G1(-a1,o)ir) 
o  3x 1  g  3x  ’  1  1 


(4A-7) 


Then  Eq.  (4A-5)  becomes 


II2  a 

I  -  —  { -  ♦o(-a1.0)G1(-a1,0;r)  -  *  (-a.,0)  —  G  (-a  ,0;r)  . 

0  g  3x'  0  1  11  01  3x'  1  1 

(4A-8) 
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I*  --  -  ^  *  g-  t— »  -A 


Comparing  with  the  expression  of  (4A-3)  we  see  that  these  just 
cancel. 


Similarly,  the  lowest  order  contributions  from  the  right  end 
point  cancel  the  other  terms  in  (4A-3). 

Thus  we  must  proceed  to  the  next  order  of  approximation  to 
Eq.  (4A-4).  It  is  readily  seen  that  this  arises  from  the  last  term 
in  (4A-4);  i.e.. 


I,  ■  /  Mr')  -i-  G  (r',r)  dz'  .  (4A-9) 

L  °  3z '  3z’  1 


Consider  the  contributions  from  the  left  end  point.  Since 
3c 

—2~  m  0  there  we  must  Taylor  expand.  This  yields 
3z' 


‘l  ■  ^2 


3z ' 


;i(e'»e) 


dz' 


x ’ *-a , 


(4A-10) 


Again  to  evaluate  the  integral  we  use  the  asymptotic  form  of 
and  extend  the  limit  to  Infinity.  The  result  is 


-fj-a,  »0) 


U 


G.(-a,,0;r) 


(4A-11) 
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5.0  TURBULENT  WAKES 


5.1  Abstract 

Turbulent  wakes  of  surface  ships  can  persist  for  considerable 
distances.  A  number  of  wind  tunnel  and  water  tank  studies  are 
available  for  wakes  of  bodies  when  cylindrical  symmetry  obtains. 
These  observations  are  supplemented  in  some  detail  by  theoretical 
calculations.  The  cases  of  cylindrical  symmetry  are  reviewed 
here.  The  SAI  FASTWAKE  code  can  model  a  surface  ship  wake  in  the 
rigid  lid  approximation.  A  review  of  some  calculations  done  with 
this  code  is  also  given. 

5.2  Introduction 

Turbulent  wakes  of  surface  vessels  may  often  be  observed  to 
extend  for  some  kilometers  on  a  relatively  calm  ocean  surface.  Such 
wakes  seem  also  to  appear  in  part  of  the  SEASAT  L-Band  radar  imagery 
of  surface  ship  wakes.  In  spite  of  the  ready  visibility  of  turbu¬ 
lent  wakes,  the  hydrodynamic  structure  and  the  interaction  with 
ambient  surface  waves  do  not  seem  to  have  received  much  study. 


In  constrast  to  this,  a  considerable  body  of  experimental  and 


theoretical  literature  is  devoted  to  studying  turbulent  wakes  of 
submerged  bodies.  This  may  be  associated  with  the  apparent  relative 
simplicity  of  a  wake  entirely  within  a  uniform,  or  slowly  varying, 
fixed  medium.  In  view  of  this,  we  are  driven  to  review  first  the 
modeling  of  submerged  wakes .  To  our  knowledge  there  is  only  one 
available  numerical  model  (the  Science  Applications  FASTWAKE  Model) 
which  can  treat  surface  ship  wakes — and  this  code  was  originally 
developed  to  study  submerged  body  wakes. 

5.3  Submerged  Body  Wake 

A  number  of  experimental  observations  in  wind  tunnels  and 
water  tanks  are  available  for  towed  and  self-propelled  bodies. 
Chevray^^  studied  the  mean  and  turbulent  flow  behind  a  spheroid. 
Naudascher^)  studied  the  flow  behind  a  disc  having  a  jet  at  the 
center  to  simulate  a  self-propelled  body.  Gran^^  studied  the  wake 
of  a  slender  body  in  a  wind  tunnel.  A  propeller  on  the  body  was 
driven  at  such  a  speed  to  give  a  momentumless  wake.  Lin  and  Pao^^ 
observed  the  wake  of  a  slender  body  in  a  towing  tank  (fresh  tap 
water  was  used).  Ten  hot-film  probes  and  20  constant  temperature 
hot-film  anemometers  were  mounted  on  a  carriage  moving  at  a  constant 
distance  behind  the  body. 


These  observations  have  provided  data  on  the  radial  and  down¬ 
stream  profile  of  mean  flow,  swirl,  turbulent  energy,  and  turbulent 
shear. 


Theoretical  analysis  of  these  observations  all  use  some  form 
of  turbulent  closure  model.  Lin  and  Pao^^  use  a  relatively  simple 
eddy  viscosity  modeling  (3  parameters),  Lewellen,  Teske,  and 
Donaldson^ employ  a  considerably  more  elaborate  model. This 
has  four  scale  parameters  which  were  taken  from  a  match  to  boundary 
I  layer  and  freejet  data  (not  re-adjusted  for  the  wake  calculations). 

The  FASTWAKE  model^^’^®^  uses  a  closure  scheme  of  intermediate 
complexity.  A  brief  discussion  of  these  models  is  given  in  the 
[  Appendix. 

In  spite  of  apparent  differences  in  the  turbulent  closure 
j  scheme,  these  models  all  show  very  similar  features  for  the  flow  and 

generally  agree  pretty  well  with  the  observed  data.  This  suggests 
that  momentum  and  energy  conservation,  combined  with  empirically 
I  scaled  transverse  mixing  rates,  may  dominat-e  the  physical  phenomena 


in  the  flow 
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To  describe  this,  we  choose  a  cylindrical  coordinate  system 


(r,  4,  z)  with  the  z-axis  being  an  axis  of  symmetry  along  the 
wake.  The  corresponding  velocity  components  are  (u,  v,  w),  with 


W  «  free  stream  velocity, 
o 

The  body  is  assumed  to  be  fixed  near  z  ■  0  .  We  write 

u  -  U  +  u'  , 
v  -  V  +  v'  , 

w  -  W  +  w'  ,  (5-1 

where  U  -  <u>  ,  etc.,  and  V  is  the  “swirl  velocity." 


Self-similar  scaling  appears  to  obtain  at  z/D  ~  10  to  10' s 
(D  -  body  diameter).  These  relations,  as  observed  in  refs.  (3)  and 
(A),  and  calculated,  are  for  a  self-propelled  body. 


w  —  r**°  *  1.0( — Z~\p$ )  *  (velocity  defect) 

°  00  D 

■  °*4( — (swirl  velocity) 
o  r-o  DC.,1' 0 


Here  CQ  is  Che  drag  coefficient  of  the  body  (nominally,  say, 

*  0.2),  z  is  the  downstream  distance  measured  from  a  point 
empirically  determined  (but  normally  within  two  or  so  body  diameters 
from  the  after  end  of  the  body). 


For  a  towed  body  the  mean  velocity  defect  and  turbulent 
velocity  decay  as  z”^^,  while  the  wake  radius  grows  as  z*^. 
In  contrast  to  the  above,  for  a  stratified  fluid 


1 


W-W 

\  w  J  /— 

o  on  axis  / z 


(5-3) 


Quantities  measured  and  calculated  as  functions  of  radius  and 
downstream  distance  are 


(w-WQ)  ,  V  ,  y< W,2>  ,  y<u»2>  ,  and  <u’w'>  . 


Detailed  comparisons  are  given  in  refs.  (4),  (5),  and  (8).  As  noted 
above,  the  agreement  of  calculations  and  observations  seems  to  be 
pretty  good.  Lin  and  Pao^  had  three  parameters  available  for 
their  model.  Lewellen  et  al.^^  had  four  parameters,  but  these  had 


previously  been  fixed  by  other  experiaents  with  different  flow 
geometry. 

5.4  Surface  Ship  Wakes 

Cylindrical  symmetry  obviously  does  not  obtain  for  wakes  of 
surface  ships.  The  portion  of  the  wake  behind  the  propeller  will 
have  a  momentum  excess.  The  turbulent  wake  from  the  hull  will  have 
a  momentum  deficit — sharing  with  the  Kelvin  wake  the  drag  force  on 
the  ship.  The  swirling  wake  of  the  propeller  will  grow,  entraining 
fluid  from  the  hull  boundary  layer.  Because  of  wind  resistance  and 
surface  wave  radiation,  the  total  wake  is  expected  to  have  a 
momentum  excess. 

The  Science  Applications,  Inc.,  FASTWAKE  is  a  three- 
dimensional,  time-independent  code  capable  of  handling  a  surface 
ship  wake  in  the  rigid  lid  approximation.  The  model  on  which  this 
code  is  based  is  described  briefly  in  the  Appendix. 

The  FASTWAKE  code  was  run  by  G.  E.  Innis*  *  to  model  the 
turbulent  wake  of  the  Quapaw  (60  m  LWL  and  D  ■  7  a  beam)  at  a  speed 
of  15  kts.  The  results  of  the  run  are  summarized  as  follows: 
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» 


(W  “  Wo ' I r«0  ~  z  for  D  <  400 


-  z"0,5  for  |  >  400 


(-  0.1  a/sec  at  3700  m) 


~  Z  l/2  for  1  <  100 


-  Z  ^  for  ^  >  100 


For  further  details  reference  (11)  should  be  consulted. 


It  Is  likely  that  the  turbulent  wake  is  sensitive  to  the 
radial  structure  of  the  hull  induced  drag  field.  This  sensitivity 
has  not  at  present  been  explored. 
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APPENDIX  TO  SECTION  5.0 


In  this  appendix  we  give  a  short  review  of  the  LP  {ref.  (4)), 
the  LTD  {ref.  (5)1,  and  the  FW  [ref.  (7)1  models.  For  simplicity  of 
presentation  we  assume  axial  symmetry  (only  the  FW  model  does  not). 


The  equations  for  the  mean  flow  are  given,  for  example,  on 
pg.  27  of  ref.  (9):  For  example,  incompressible  fluid  of  uniform 
density  the  equation  of  continuity  reads 


7  W  <*» 


0 


(5A-1) 


If  f  ■  (f^,  f^)  represents  a  body  force  acting  on  unit  mass  of 
fluid,  the  equations  of  motion  are 


i  II 

p  3z 


-4r<»,2>  (t  <»’*’>) + ■ 


dz 


T  3lf 


(5A-2) 


w  n  IX  +  w  .  _  J_  rY2<u'v’>) 

W  3z  U  3Y  Y  ^2  3y  j 


* 


(<v’w'>)  +  f 


w  3l+  U  3 Y  -  Y - pTy  -7  37  KU 

+  <v^!>  _  |_  (<u.w.>)  .  (5A-4) 

Here  p  is  the  fluid  density  and  P  is  the  mean  pressure. 

Equations  to  determine  <u'v'>  ,  <u'^>  ,  etc.,  are  reviewed  in 

Hinze.(^)  These  equations  contain  third  order  quantities 
3 

<u '  >  ,  ,  etc.,  for  which  there  is  another  set  of 

equations  involving  fourth  order  averages,  etc.  All  models  employ 
some  closure  scheme. 

Relations  expressing  conservation  of  momentum  and  angular 
momentum  may  be  obtained  directly  from  (5A-3)  and  (5A-4). 

First,  we  re-write  (5A-3)  as 


fc  (*> +  h  b 


-71? 


-  -5—  <v'w’>  +  f  . 


(5A-5) 


Let  the  propeller  be  at  z  -  0(+)  and  integrate  (5A-5)  as 


2irp  /  dz  /  y  dY  •  •  • 
o  o 


The  results 


2irp  /  VW  Y2  dY  +  2wp  /  <v'w'>  Y2  dY 
o  o 


n 

2rp  /  i  dY  /  dz  f .  *  torque  of  propeller 


Now.  In  the  far  wake,  W  “  W  and  <v'w’>  Is  a  second  order  small 

o 

quantity.  Also,  the  torque  Is 


n_P 


a) 


(5A-6) 


where  P  Is  the  power  delivered  to  the  propeller,  n  is  the 

P 

propeller  efficiency  (nominally,  in  the  0.6  -  0.8  range)  and  is 
the  angular  velocity  of  the  propeller. 


Thus,  in  the  far  wake  we  have 


2irpWQ  /  Vy  aY  ■  t  , 
o 


( 5A-7 ) 
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*  >  .  •  w  « 

-V- 
1 ... 


*.\V. 


v'/.V 


,  /, 

.V.'-* , 


a  rather  obvious  relation,  since 


2irp  /  (Y  V)  Y  dy  ■  angular  momentum/ unit  length 
o 

and  x  is  the  rate  of  generation  of  angular  momentum  in  the  wake. 
The  swirl  in  the  wake  is  then  constrained  by  the  power  delivered  to 
the  propeller. 


To  describe  momentum  conservation,  we  re-write  (5A-2)  as 


—  f  w2  1  +  —  —  (  y  uw]  +  —  — 

3z  V"  J  y  3y  ^  UWJ  p  3z 


1  3P 


(5A-8 


3z 


<W  >  **  -  (r  <u'w'>)  +  fa 


He  now  perform  the  Integral 


2irp  /  dt  /  Y  dy  .  .  . 


z  o 

o 


where  zQ  is  a  point  ahead  of  the  body.  This  gives  the  relation 


2irp  /  (W2  -  W2)  y  dy  +  2ir  /  (p  -  PQ)y  dy 


+  2*  p  J  <w'2>  Y  dy 
o 

■  2sp  /  f  dz  y  dy  •  total  force  acting  on  wake, 
z 

Here  Pq  is  the  ambient  pressure  ahead  of  the  body.  We  also  have 

2sp  /  fz  dz  Y  dy  -  FT  -  |  p  CD  Ag  W2  .  ( 5A-9 ) 

Here  is  the  drag  coefficient  of  the  body  and  Ag  is  its  area 

(»  it  (D/2)  ,  D  ■  diameter).  The  thrust  force  of  the  propeller  is 

written  here. 

T 

In  the  far  wake  we  take  P  -  P  ,  W2  -  W2  *  2W  fw  -  wj 

o  o  o  o 

and  have 

m 

Asp  Wo  /  (w  -  Wo)  Y  dy  -  FT  -  i  p  CD  Ag  w2  .  (5A-10) 

o 

For  a  se If -propel led  submerged  body  the  right-hand  side  of  (5A-10) 

vanishes.  For  a  towed  body  F^  ■  0  ,  while  for  a  surface  ship  we 
1  9  ' 

expect  FT  >  y  p  Cg  Ag  Wq  .  The  value  of  the  right  side  of 
(5A-10)  has  a  very  important  implication  for  the  development  and 
persistence  of  the  wake. 


A  further  integral  constraint  on  the  wake  is  given  by  Birkhoff 
and  Zarantonello^®^ : 


2ir  W  /  (w  -  W  )  Y3  dY  “  constant 
o 


The  evaluation  of  the  turbulent  energy 


1  >  ,2. 

q  -  "2  <v'  > 


(5A-11) 


is  described  by  [see  Hinze,  Eqs.  (1-111)]. 


fW  TRT  «  +  U  T7 
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+  [<*'2>%+«'2>%+<''2>Z 


+  <u'w'>  +  ff)  +  <u’v’>  Y  (^)l 


3  ,V, 


3z 


3u'  3u'  3u' 

V  <3xJ  ^3x7“ 


+  3^)>  5  "e  » 

j  i 


(5A-12) 


where  we  have  used  rectangular  coordinates  to  express  the  turbulent 
dissipation  rate  e  . 


The  three  bracketed  terms  on  the  right  describe,  respectively, 
advection,  production,  and  diffusion  of  turbulent  energy. 

The  LP  model^)  uses  the  simplest  closure  prescription.  They 

write 


2  2  2 
<u'  >  -  <v'  >  -  <w'  > 

«*'»•>  -  -  vT  -JJ  (w  -  Wo)  , 

<"'v'>  ■  -  u.  T  Tf  ’ 

■  <u'(|  V'2  +  )>  .  (5A-13) 

The  three  quantities  *  v8»  and  ve  were  determined  by  a  fit  to 
observed  data.  Generally,  the  final  fit  seemed  good. 


A  much  more  elaborate  model  is  used  by  LTD^^’^^.  They  write 


<“i 


p’>  + 


3 

3x, 


<UJ 


P’> 


(5A-14) 


3Xk  i  -J" 


[In  FW'  '•  It  is  assumed  that  <uj!  p’>  ■  0  ,  but  a  more  symmetric 

expression  is  used  for  the  right-hand  side.] 


3u!  3u!  ov 

2u  <T-i  3-J->  -  It1  5„  +  <uj  u.’> 

S*k  S\  3A  13  A2  1  i 

,  3u!  3u!  rz —  0 

<p"  ^  +  =  ■  nF  (<ui  uj>  ■  f q  6ij 


(5A-15) 


)  • 


(5A-16) 


,7.  3U.  3U 

[FW3''  adds  a  term  Cjq(-g^—  +  ■g~)  t0  the  right  of  (5A-16)]. 


The  parameters  are  determined'®''  from  data  on  boundary  layers 


and  jets: 


2.5  , 

|  ,  (5A-17) 

0.3  , 

0.2  n,  . 


LTD  insert  the  relations  (5A-14),  (5A-15),  and  (5A-16)  into 
the  differential  equations  for  <uj  uj>  ,  which  then  provides  a 
closed  set  of  equations.  LTD  argue  that  it  is  consistent  with  the 


5 


data  to  set  C^-  0  and  that  the  more  symmetric  form  of  (5A-14)  is 
unnecessarily  complicated. 

A  significant  simplification  is  introduced  in  FW^).  It  is 
assumed  that  deviations  from  isotropy  are  "small"  and  a  set  of 
coupled  algebraic  equations  is  obtained  for  the  <uj  uj>  . 

The  similarity  in  the  results  obtained  from  these  three  models 
suggests  that  conservation  of  momentum,  angular  momentum,  and 
energy,  coupled  with  a  "reasonable”  model  for  radial  transport, 
determines  the  gross  features  of  wake  development. 
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6.0  A  FUTURE  EXPERIMENT 


After  three  years  of  theoretical  effort  without  firm  conclu¬ 
sions,  it  is  time  for  an  experimental  effort. 

1.  The  first  step  was  a  JSS  83  experiment. ^  We  photographed 
a  number  of  ship  wakes  off  Southern  California  from  the  SIO  air¬ 
craft,  at  elevations  between  S00  and  1500  feet.  The  principal 
result  is  that  different  ships  have  different  wakes.  Sometimes 
the  +  19.5°  cusp  stands  out  very  clearly;  it  may  be  the  outer  edge 
of  a  broad  Kelvin  pattern,  or  simply  two  narrow  cusp  bands  without 
much  wave  structure.2  At  times  the  +  19.5°  cusp  was  not  visible, 
and  the  interior  wake  was.  The  turbulent  wake  can  be  single  or 
double  (depending  on  the  number  of  propellers?),  and  can  be  visible 
as  a  slick  (calm)  region  in  a  rough  ambient  sea,  or  as  a  roughened 
band  in  a  calm  background.  Winds  clearly  make  a  difference. 

Sun  glitter  can  be  used  to  advantage  by  photographing  the  wake 
at  the  edge  of  the  glitter,  where  a  slight  favorable  slope 


*R.  Davis,  pilot;  F.  Zacharlasen,  navigator;  W.  Hunk,  photographer. 

2This  could  be  because  the  cusp  attenuates  downstream  as  x'1^ 
compared  to  x-1'2  for  the  interior  wake. 


enhancement  from  the  Kelvin  waves  can  produce  an  enormous  increase 
in  specular  return.  In  this  way  we  could  see  many  wave  lengths  of 
short  divergent  waves  at  very  long  distances  behind  the  vessel. 

2.  We  participated  as  observers  in  a  joint  U.S. -Canadian 
experiment  led  by  Blyth  Hughes  of  the  Canadian  Defense  Establishment 
on  Vancouver  Island.  The  wake  produced  by  the  NOSC  Vessel  QUAPAW 
was  measured  by  equipment  forward  of  the  bow  of  the  ENDEAVOUR  during 
wake  crossings.  The  key  device,  designed  by  Hughes,  is  an  upward- 
looking  beam  from  a  submerged  laser  mounted  on  the  bow  of  the 
ENDEAVOUR.  The  beam  is  deflected  by  the  tilted  surface  facets,  and 
the  deflection  is  recorded  in  two  dimensions  against  time.  The 
illuminated  spot  is  a  few  mm,  and  the  response  time  is  a  few  ms,  so 
that  wavelengths  down  to  a  cm  can  be  measured.  Open  sea  measure¬ 
ments  using  Hughes'  tilt  device  were  published  some  years  ago.  The 
slope  statistics  so  obtained  were  in  very  good  accord  with  slope 
statistics  derived  from  sun  glitter. 

A  very  crude  ship-board  analysis*  indicates  an  amplitude 
pattern  roughly  consistent  with  thin  ship  theory  for  a  traveling 

*Thl8  was  summarized  in  a  letter  of  5  August  1983,  from  Munk  to 
Zacharlasen,  Dashen,  and  Davis. 


source.  However,  the  transverse  slope  dc/ty  falls  to  grow  with 

decreasing  y  ,  In  contrast  to  the  theoretical  expectation 
-1/2 

3C/3 y  -  y  ;  in  fact,  significant  slopes  occur  only  just  Inside 
the  19.5°  cusp. 

It  was  also  learned  that  the  wake  slopes  were  swamped  by 
ambient  slopes  in  the  presence  of  even  light  winds.  For  future 
experiments  we  shall  need  to  suppress  the  larger  contribution  to 
slopes  from  the  ambient  waves  that  are  shorter  than  the  ship  wake 
waves . 

3.  It  is  now  fairly  clear  that  the  Interior  (to  the  +  19.5° 
cusps)  behavior  of  transverse  slopes  Is  sensitive  to  the  type  of 
vessel,  including  hull  shape  and  configuration  of  screws.  He 
suspect  that  the  narrow  V's  seen  on  the  SAR  images  were  associated 
with  particular  vessels  that  produce  high  Interior  transverse 
slopes.  Our  first  recommendation  is  to  perform  wake  calculations 
for  various  ships  and  to  select  some  that  seem  suitable  from  the 
foregoing  considerations.  It  would  probably  be  advantageous  if  some 
suitable  Naval  vessels  could  be  identified. 


4.  With  regard  to  a  future  experiment,  we  propose 


(a)  To  make  simultaneous  coordinated  measurements  of  a 
Kelvin  wake 

(i)  with  a  SEASAT-like  SAR  from  a  spacecraft. 

(11)  optically  from  an  aircraft. 

(ill)  with  a  Hughes  device  from  a  vessel  (suitably 
low-passed). 

It  might  be  worthwhile  to  Include  coordinated  photography  from  the 
spacecraft;  it  might  also  be  possible  to  fly  a  second  SAR  on  the 
aircraft. 

(b)  To  cover  a  variety  of  parameters: 

(i)  Ship  speed  and  direction. 

(li)  At  least  two  types  of  vessels,  which  produce 
significant  Interior  transverse  slopes. 

(ill)  Wind  speed  and  direction  (relative  to  vessel), 
(lv)  SAR  viewing  angles  of  15°,  20°,  25°,  with  the 
vessel  (centered  on  SEASAT's  20°). 


(c)  The  experiments  should  be  closely  tied  with  numerical 


simulation 


It  ought  to  be  possible  in  such  an  experiment  to  identify 
the  V  signature  on  SAR  images  with  other  observables.  Once  this 
has  been  done,  the  theoretical  effort  can  be  focused  and  will  no 
doubt  produce  early  affirmative  results. 

5.  An  opportunity  for  such  an  experiment  exists  in  connection 
with  the  proposed  shuttle  flight  in  August  1984.*  The  very  purpose 
of  the  experiment  is  to  gain  an  understanding  of  radar  imagery  of 
land  and  sea  (including  wind  waves,  internal  waves,  .  .  .  ),  so  that 
a  study  of  Kelvin  waves  should  be  a  natural  part  of  the  experiment 
(they  are  not  included  as  of  now). 

There  would  be  a  total  of  about  10  passes  over  a  given  spot; 
this  is  a  serious  limitation.  The  desirability  of  placing  the  ship 
wake  at  the  edge  of  sun  glitter  will  call  for  careful  orbit 
planning.  Because  of  the  possibility  that  shuttle  flights  are 
postponed  or  canceled,  it  is  proposed  that  the  experiment  be  planned 
to  provide  critical  information  even  without  the  spacecraft. 

We  recommend  two  dedicated  wake -making  ships,  operating  as 
close  as  possible  without  wake  interference.  They  should  be  on  a 


*The  SIR-B  Plan.  JPL  82-78,  December  1982 


straight  course  for  at  least  half  an  hour  prior  to  satellite 
passage.  A  third  vessel  is  required  for  wake  crossings.  Aircraft 
glitter  (and  SAR)  measurements  and  surface  tilt  measurements  should 
be  made  not  only  during  satellite  passage  but  at  other  times  as 
well. 


Here  is  a  possible  expedition  profile: 


Satellite 
viewing  angle 

Two  ships' 
speed 

Cover  relative 
to  satellite 

Wind 

speed 

20° 

fast 

parallel 

low 

20° 

slow 

parallel 

low 

20° 

fast 

parallel 

high 

20° 

slow 

parallel 

high 

20° 

fast 

perpendicular 

low 

20° 

fast 

perpendicular 

high 

15° 

fast 

parallel 

low 

25° 

fast 

parallel 

low 

This  just  about  exhausts  the  satellite  measurement  program.  But 
many  other  configurations  should  be  tested  by  aircraft  and  surface 
measurements . 
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We  recommend  a  site  near  Maui,  to  provide  a  high  probability 
of  clear  sky,  and  a  variety  of  sites  with  (normally)  different  wind 
speeds . 

There  is  some  urgency  if  this  program  is  to  be  carried  out  one 
year  from  now.  A  responsible  group  needs  to  be  identified.  One 
possibility  is  APL  Johns  Hopkins,  who  have  experience  in  this  field. 
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